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FOREWORD 


This report was prepared by 1/Lt Thomas 0. Dobbins of the Rocket 
Propellant Section, Fuels and Oils Branch, Propulsion Sub-Systems Division, 
Propulsion Laboratory. Work was performed under Project 3148, "Development 
of Liquid and Selid Rocket Propellante". This work was also submitted by 


the author 9s a thesis for the degree Mastei' of Science at the Ohio State 


University, Desember 1959. 
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This report presents a concise sumiary of the thermodynamics of rocket 
propulsion and the evaluation of 100 prototype propellant combinations, 
Emphasis is placed on the thermodynemics of chemical propulsion and on flow 
through the converging-diverging nozzle. Fran this basic theory, two Y2)*—— 
rigorous methods are developed for conducting theoretical serformance 
calculations. One. of these utilizes thermodynamic equilibriun 
constents while the second utiiizes a completely new concept cf minimizing 
the total free energy of the combustion products, Sample caleculatians are 
presented for each method. Thermodynamic date are presented and discussed 
for several reactants and ell important products of cambustion consisting 
of the elements carbon, hydrogen, oxygen, nitrogen, chlorine, fluorine, 
lithium, boron, beryllium, alumi.um, and magnesium. A rigorous theoretical 
evaluation is presentec for 100 prototype propellent aysteim. Syecific 
impulse, chamber temperature, exhaust temperature, characteristic exisust 
velocity, and chamber moleculer weight are plotted as e function of the 
propellant compositicn. Also, complete serformanse data are tabulated for 
selected fuel-to-oxidizer ratios. 4 
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PREFACE 


In this report, the author hee attempted to present a concisa, 
comprehensive, fuadamental approach to the theoretical evaluation of rockat 
propellant performance. Firat, those basic principles of thermodynamics end 


fluid dynemica directly related to rocket propulsion are reviewed, Although 





many detailed textbooks have been prepared on these two subjects, fewer 
attempts have beer made to combine the phasea applicable to rocketry into 
one summary. From this fundamental theory, two rigorous methods of 
calculating propellant performance are presented. 


One of these approaches, the minimizetion-of-free_energy at constant 





enthalpy or constant entropy, has not been published previously. Then, 
using the latter method of calculation, the author has evaluated the 
performance of a series of prototype fuels and oxidizers. The basic 
thermodynamic date used in this evaluation are also presented and discussed. 
Many scientists have contributed indirectly to this report by their 
work in the basic areag of thermodynemics, fluid dynamics and propulsion. 
Unfortunately, there is no adequate way to repay this debt except through 
& general acknowledgement. ‘When material hes been taken from the 
publications of otbers, this has been ind&éated by footnotes in the text. 
The author's first specific acknowledgement is to Professor P. R. Jones of 
West Virginia University and Dr. W. B. Kay of the Ohio State University who 
introduced him to the acience of thermodynemics and kindled his interest in 
the subject. Furthermore, en additional acknowledgement is made to Dr. Kay 
for his guidance and interest while serving as adviser to the author on 


this project. Also, Dr. R. S. Brodkey of the Ohio State University is 
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recognized for his inspiring teaching of fluid aynenton. 

Substantial assistance wes proviied ey several others in prepsring, 
specific sections of this report. The author wishes to acknowledge 
Mr. A. E. Roberts, Jr. of General Kinetics, Inc. who programmed the 
minimizeation-of-free energy concept far the Universal Scientific Canput=r, 
Model 11034, and who offered valuable constructive camments on the author's 
presentation of this theory. Thanks are also due Mr. L. J. Gordon of the 
Aerojet-General Corporation for his guidance in evaluating the relative 
importance of propellant density and specific impulse. The author is 
grateful to Dr. C. B, Henderson and Mr. R. S. Scheffee of the Atlantic 
Research Corporation for their valuable assistance in assembling 
thermodynamic data on reactants and products-of-cambustion. Additional 
recognition is given to Dr. D. Re Stull of Dow Chemical Company, Mr. /chn S. 
Gordon of Reaction Motors Division of Thiokol Chemical Company, and to 
Dr.. Charles W. Beckett, 7.» and his colleagues at the National Bureau of 


Standards for their assistance and advice in obtaining accurate thermodynamic 





deta. The unpublished multiphase flow theory of Dr. David Altman of United 
Aircraft Research also contributes to the comprehensiveness of this report. 
The interest of each of these gentlemen is sincerely appreciated. 

Finally, but certainly not least, thanks are due to Miss Dora M. Wolfe 
and Mrs.. Lavera Howard of Wright Air Development Center for typing the text 
of this report, and to Mr. Nicholas Engler and Mr. William Powers of the 


University of Dayton for assisting with the plotting of graphs. 
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CheFTER I 


INTROUUCTION 
A, Historical Background 


Although the concept of rocket propulsion is an old one, in practice 
it is a relatively new aeroneutical science. Professor Robert H. Goddard, 
hailed as the father of American rocketry, began his atudies in the year 
1935. For the next 25 years, he and his followers pioneered rocket 
research within the United States until the outbreak of World War II 


stimulated considerable interest in their activities. Since that time, 


| 
rocket research and development has been pursued at an accelerating rate 
until it now constitutes a significac. part of the activity cof America's 
aircraft industry. 
Fran the time of Goddard's early experiments, rocket scientists have 
needed to know the performance of a rocket propellant syatem before design- 
ing « missile to use it. Originally, thermodynemists thsoretically 
calculated the performance of a few carbon-hydrogen-oxygen propellants by | 
neni. However, the complexity of the calculation coupled with the lack of | 
high temperature thermodynamic date limited their acticity. Aa same ! 
rocket thexrmodynemists attempted to expand tne scope of their work, short-cut | 
methods of calculation were developed. However, these techniques invoked | 
simplifying assumptions which introduced considerable error e# the chemical | 
complexity of the propellants increased, particularly if condensed phases | 


wero formed as products of combustion. 
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In the cerly 1950's, investigutors at the Lowis Laboratories of the 
Nutionul Advisory Comaittee for Aeronautics developed a program for celcula- 
ting the performance of propellants composed of carbon, hydrogen, oxygen, 
nitrogen, chlorine, and fluorine on the IBM 650 digital computer. Similarly, 
investigators at California Institute of Technology's Jet Propulsion Labore- 
tory programmed the Datatron 201 computer. This constituted an advencement in 
theoretical analysis work, However, these and many similar programs subse- 
quently developed by others have been designed to evaluate specific ingredi- 
ent elements which form specific products of combustion, In these cases, 
neither new ingredient elementa nor new products of combustion may be con- 
sidered without re-writing the computer program. Furthermore, the shortage of 
memory capacity in amall and me*’ w-size digital computers hes limited the 
complexity of the propellant system which could be evaluated. 

4s an extension and expansion of this earlier performance calculation 
work, a completely general computer progrem based upon a new concept of the 
minimization of free energy has been established at Wright Air Development 
Center. This program, preparsd for the Universel Scientific Computer, Model 
1103A, is neither limited by the number or type of elemental ingredients nor 
by the number or state of the products of combustion. Vor all practical pur- 
poses, any propellant system for which thermodynamic data are available may 
be evaluated, This achisvement may be considered a logical outgrowth of the 


work of earlier investigators. 


Be. Object and Scope 
This report has been prepare as a contribution to the theoretical 


evaluation of potcntial fuels and oxidizers es propulsion media. The objest 
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of the report is to present the thermodynemics of rocket propulsion and 

the evuluation of several prototype propellant combinatioas., The scope 
may be divided into four logical areas: (1) a concise, canprehensive 
treatment of the fundamental theory of rocket propulsion; (2) a detailed 
description of two rigorous methods for calculating propellent performance; 
(3) « discussion of the thermodynamic data for a number of selected 
reactants end products of cambustion; and (4) the evaluation of the 
performance of several prototype propellant systems. 

The development of fundamental theory is designed to incorporate all 
important concepts into one logical, comprehensive summary, The literature 
has been carefully screened in order to utilize the contributions of 
previous workers. Emphasis is placed on the thermodynamics of chemical 
propulsion and on flow through the converging-diverging nozzle. From this 
pabietheeen, two rigorous methods are developed for conducting theoretical 
performance calculations. One of these utilizes thermodynamic equilibrium 
constants to express the concentrations of minor products of combustion in 
texms of the major species. The second method utilizes a completely new 
concept of minimizing the totel free_energy of the canbustion products. 
Although b-th methods are equally rigorous, the latter one is more general 
end wes aj:iied v5 the calculational work. 

Thermodynamic data are presented and diacussed for several reactants 
and ali important products of cambustion consiating of the elements carbon, 
hydrogen, oxygen, nitrogen, chlorine, fluorine, lithium, boron, beryllium, 
aluminum, and magnesium. Ten protot-pe fuela and ten prototype oxidizers 


were then chosen for rigorous theoretical evaluation. Although a number of 











these systems do not constitute workable propellant combinations, they were 
selected because of the basic combustion information wich could be obtained, 
The study was designed to reveal the relative merits of various elements 

and chemical bonds as ingredients for chemical propellants. In conducting 
the calculations, all known products of cambustion were assumed to exist. 
Since the calculational technique was completely rigorous, the accuracy of 
the calculations in this report is limited only by the accuracy of the 

basic thermodynamic date used, Therefore, this report is sutmitted as a 
comprehensive resume’ of the present state-of-the-art in theoretical 


propellant evaluation. 
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GHAPTER IT 
SUMMARY 

The evaluation of shemicel rocket propellanta may be based directly 
upon the laws of thermodynamics. The energy available from any chemicai 
propellant combination is directiy proportional to the square root of 
the change in enthalpy of the system as it flows through a converging- 
diverging nozzle. This analysis is equally applicable to both campress#- 
ible and multiphase flow. In the latter situation, thermal and manentwm 
equilibrium between the gas phase and the condensed phase are very 
closely approached when the diameter of the condensed particle is less 
than ten microns. If this size is exceeded, the momentum lag may 
cause a losa of uszcful enersy. Although the density of a propellant 
eaubination must be considered in evaluating its effectiveness in a 
given missile application, the specific impulse is the most significant 
parameter in screening propellants. 

Rigorous performance calculations mat be made by assuming adiabatic 
caabustion in the chamber and isentropic expansion through the e nozzle. 
This may be accomplishad by either the method of equilibrium constants on 
the. method of minimization-of-frea energy. The latter is recomended for 
conducting calculations on a high-speed digital computer because of its 
completely general approach. 

High temperature thermodynamic data are available for the products 
of combustion of most of the interesting propellant ingredients. However, 


in meny cases these data are extrapolated fran experimental studies at 
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lower temperatures. Therefore, additional experimental studies ara 
needed to refine the thermodynamic properties of many substances. 

This is perticularly true of the cowpoundsj of the light metals lithium, 
aluminum, beryllium, boron, and magnesium. 

Fluorine and oxygen difluoride were found to be excellent oxidizers 
for most of the prototyyge fusla studied, Hydrogen and beryllium hydride 
were universally energetic fuels. For oxidizers, the O-7 bond wos 
found more energetic than the N-F bond which, in turn, surpassed the 
Gl-F bond. For fuels, hydrogen wes found to be an essential ingredient 
for high energy. Nitrogen appears more desirable than carbon as « 
carrier element. Eeryliium appears to tbe ths most prawisine ef the 
light motels, although boron amd aluminum offer considerable energy. 
Lithium is considerably more energetic if oxidizad with fluorine than - 
2f oxidized with oxygen. 











CHAPTER IIT 


FUNDAMENTAL THECRY OF ROCKET PROPULSION 


&  Bocket Propulsion Principles 
1. Genera) Principles 
All known methods for propelling a body through a fluid medium 
are based upon the reaction principle. One of these methods is rocket 
propulsion, which may be defined as a means of locomotion through which 
reaction is imparted to a vehicle by the sseuati bias of ejected matter. 
The matter is ejected as a high-speed gas stream, which is generated by 


tween a fuel and ai: oxidizer stored within the 
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vehicle, By accelerating the gas stream through a converging-diverging 
nozzle, supersonic velocities sre attained. Thus, rocket propulsion is 
achieved by the conversion of the thermal energy of a chemical reaction 
into the translational energy of the products of combustion, 

The theory of propulsion was first presented by Sir Isaac Newton 
in the year 1686 with the publication of his three fundamental laws of 
motion in Phi Naturalis Principia Mati tica. These laws 
state that: (1) every body tenés to continue either in a state of reat 
or in a state of uniform mit — unless acted upon by external forces; 
(2) a change of motion is proportional to the force introducing the 
change and occurs in a straight line along which the force acts; and 
(3) the action of a force is a)ways oppoged by an equal and opposite 
reaction. Newton’a second lew of motion constitutes the fundamental 
principle of fluid dynamics. CGombining the second and third laws, the 


magnitude and direction of all momentum changes are defined. 





thermodynamics, and the second law of thermodynamics, These laws are 
completely independent of the nature of the propelling medium. For a 
specific rocket propulsion process, additional subsidiary lawa may be 
introduced which are related to the particular nature of the propelling 
medium, These include the equation of state and equations expressing 


the thermodynamic properties of the fluid. 


2. Tke Momentum Principle 





In addition to Newton's laws, the behavior of compressible fluids 
is also governed by the law of conservation of mass, the first law of 
‘From the laws of motion, momentum is shown to be a vector quantity 


possessing the direction of the velecity vector. The time rate of | 
' change of the total momentum of a bounded body of fluid is equal to the 


resultant of the external forces acting on the boundaries, i.e. 


} LF! = $~ (av) (1) 
| 

| 

t Tu. net external force is also a vector quantity having the same magni- 
if 


tude and direction as the velocity vector representing the change in 


the rate of momentun, 


| 

In applying the principle of momentum to rockst propulsion, | 
another term must be added to the momentum equation when a difference | 
in pressure exists between the emergent fluid jet and the surrounding | 
fluid medium, Thus, the total thrust imparted to a vehicle moving | 
| 


through a homogeneous external fluid is composed of the momentum thrust 


rate of propellant and the exhaust velocity of the smergent fluid strean 


and the pressure thrust. The former is the product of the mass flow 
\ 
relative to the wehicle. The latter is the product of the cross-sectional 


8 





area of the emergent stream and the difference in pressure between this 


stream and the surrounding medium. Mathematically, 


where fF! 
m 


& 

Ve 
Pe 
Pa 
he 


condition defines 


through a converging-diverging nozzle, and will be discussed in more 
detail in a subsequent section. When the nozzle exit pressure exceeds 
the external pressure, the nozzle area is not sufficiently large to 
permit complete expansion of the working fluid. On the other hand, if 
the exhaust gases are overexpanded to a pressure less than the external 
pressure, @ shock wave will occur in the nozzle. Thus, either over- 
expansion or underexpansion induces a loss of energy and a corresponding 
reduction in the exhaust velocity. 

Actually, only the axial component of velocity will produce a 
net thrust in a converging-diverging nozzle. Therefore, in compering 


experimental and theoretical thrust values, a geometric correction must 


be applied, i.e. 


Fi= 5 (1 + coset) avg + (De-Pa) Ae (3) 


F'=2Bve + (po-Pa) Ae (2) 


m 


thrust 

mass flow rate 

acceleration of gravity 

velocity of exhaust gases at nozzle exit 


nozzle exit pressure 
external surrounding pressure 


exit area of nozzle 


the optimum thrust obtainable from isentropic flow 











where ©4 is the helf angles of nozzle divergence. Thus, the effective 
exhaust velocity for the general case is defined as 


(Pe-Pa Ags 


= Fig il 
via Shas (1+ cose) ve + = 


(4) 


This, in turn, leads to the definition of specific impulse, which is 


one of the cost important parameters for evaluating rocket performance, 


3- Useful Rocket Parameters 
Specific impulse is defined as 


Isp= = (5) 


or the pounds of thrust obtained per pound of propellant flow per second. 


[ 
| 
| 
| 
| 
| 


In practice, the dimensions of specific impulse are reducad to seconds. 
It will be shown later that specific impulse is a property of the 
chemical composition of the propellant and can be theoretically 
evaluated by thermodynamic analysis. 

In ordinary engineering practice, variations in gravity are un- 
important, and mass aid weight quantities are used interchangeably. 


Although the thrust of 2 rocket is actually a function of the mass flow 


rate (not of the weight flow rete), weight terminology is still preferred 


for use sy the rocket industry. Actually, specific impulse is independ-~ 
ent of gravity and other external effects and, in mass terminology, 
would have the dimentions of ft/sec and be numerically equal to the 


effective exhaust velocity. 
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In most sxperimentai rocket meter tests, the following quantities 

are either measured or known: 

(1) Thrust, F 

(2) Combustion pressure, pe 

(3) Exit pressure, De 

(4) Nozzle throat area, A* 

(5) Nozzle exit area, Ag 

(6) Mass rate of propellant flow, m 

(7) Propellant composition 
The overall performance, as measured by specific impulse or exhaust 
velocity, is evaluated experimentally from items (1) and (6) or computed 
theoretically from items (2), (3), and (7}. However, in some test work 
it is desirable to individually analyze the combustion and the expansion 
processes. The combustion process should determine the ability of a 
propellant combination to generate useful energy while the expansion 
process should determine the effectiveness of utilizing that energy. 


For this reason, two experimental quantities, the characteristic velocity 


and the nozzle thrust coefficient, are defined respectively as 


oe a PAte (6) 


x 
and Gp = ne (7) 
The characteristic velocity is actually the velocity of sound at the 
nozzle throat. In essence, this is a measure of the effectiveness with 
which the combustion products are generated in the rocket chamber. 
Closely related to the characteristic velocity is the weight flow 


coefficient, or discharge coefficient, which is defined as 


ey * oe (8) 








This parameter gives a direct measure of the rate of consumption of 
the propellant. The product of the thrust coefficient and the character- 
istic velocity equals the exhaust velocity. 

The total impulse may be expressed either as a function of thrust 


or of specific impulse: 


T= Ft = Tep + W, (9) 
where I = total impuise 
t = operating time 
Wy = total weight of propellant consumed 
Isp = specific impulse 
F' = thrust 


This is a measure of the total energy possessed by a given missile, and 


che 


may consist of a large thrust lasting 


z 
» 


‘oz a short duration or a smaller 


t 


thrust acting for a longer time.. 


B. Analysis of the Ideal Rocket 


1. Fundamental Technical assumptions 
The quantitative treatment of fluid dynamics may be based 


directly upon the laws of thermodynamics .4 4& strictly thermodynamic 
approach Goes not require any consideration of the mechanism of flow. 
Rather, it constitutes an energy analysis and establishes the limiting, 
ideal conditions of flow. A frictionless, one-dimensional flow is 
assumed which is independent of fluid viscosity, the transport property 


of momentum. This approach is quite justified from an engineering view- 


point since the viscous forces are negligible compared to the inertial forces, 
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Frequently, the measured performance of rocket engines exceeds 95% 

of the ideal, theoretical valves. Furthermore, the accepted practice 
for designing rocket engines is to utilize ideal rocket perameters and 
modify these by empirical corrections. 

Before introducing the equations of thermodynamics which govern 
rocket propulsion, the assumptions invoked in analyzing the canbustion 
and expansion processes will be discussed. The thermodynamic approach 
must possess a maximum sophistication in order to place the minimm 
emount of limitation upon the theoretical resulta, Therefore, assumptions 
are limited to those justified by theoretical or experimental evidence, 
those introducing a negligible effect on the final results, and those 
representing the iz al limit for chemical and physical phenomena which 
lack experimental verification. These assumptions are the following: 


= we 


(1) The gaseous species behave ideally, i.e. obey the 
perfect gas. law. 


(2) Dalton's Law applies. 


(3) Molar enthalpies of all species are independent 
of pressure... 


(4) Molar entropies of condensed phases are independent 
of pressure. 


(5) Steady-state conditions exist during combustion 
and expansion. 


(6) Condensed phases occupy a negligible volume. 


(7) The entropy of mixing is negligible for condensed 
. phases. 


(8) Vapor pressures of condensed phases are independent 
of total preseure. 


(9) The combustion process is adiabatic, 
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(10) Combustion occurs at constent pressure. 


(11) Thexmal and chemical equilibria exist in the 
combustion chamber, 


(12) Homogeneous mixing is attained. 


(13) The canbustion products in the chamber have zero bulk 
velocity with respect to the walls of the chamber. 


(14) Isentropic, frictionless, one-dimensional flow occurs 
in the nozzle. 


(15) The exit velocity of the combustion products is uniform 
across any section normal to the nozzle axis. 


(16) The gases are insoluble in the condensed plese, and the 
condensed phases themselves are immiscible, 


(i7) Ths particles are thermodynamically macroscopic. 


(4.0) Condensed phases are in thermal and velocity equilibrium 
with the ges stream curing the expansion process, 


2. Discussion of the Asswnptions 


The first five of these assumptions may be considered as justified 
by existing experimental evidence. ‘The only questionable point appears 
to be the assumption of ideal behavior for the vapora of the metal 
oxides. Because of their high boiling points, these compounds could 
conceivably deviate somewhat from ideality at combustion conditions. 
However, the relatively amei}, concentrations of these vepors in the 
combustion products will tend to minimize the significance of any 
possible non-ideal behavior. The ideal gas approximation has been found 
to be quite adequate for the more common canbustion products. 

Assumptions six, seven, and eight are made to simplify the calcu- 


lation procedure. In these cases, correction is not particularly diffi- 


cult, but the effect upon the final results is negligible. 
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Any error introduced hy these assumptions is certainly insignificant 
compared to uncertainties in the basic high-temperature thermodynamic 
data themselves. 


Assumptions nine through fifteen are peculiar to the rocket 


.f 
Z 


propulsion process itself, i.6, an adiabatic combustion followed by an 
isentropic expansion through a converging-diverging nozzle. These 
assumptions represent the ideal, limiting conditions of snergy utili- 
zation and overall efficiency, and are approached experimentally in 


large rocket engines. In practice, the geometry of the rocket engine 





may affect the attainment of squilibrium. A short reaction chamber may 
not afford a sufficiently lohg residences tims for some chemical reactions 
to reach equilibrium, In the nozzle, the attainment of equilibrium could 
conceivably be inhibited by the nature of the chemical reaction and the 
existing physical conditions. For example, associative reactions would 
be less Likely at lower pressures when fragments would be more remote 
and their probability of collision relatively less. The one-dimensional 
flow approximation is an excellent one in short ducts of variable cross- 
sectional area. In this case, the fundamental cause for the change in 
flow conditions is the variation in the area, and other causes, such as 
boundary layer development and skin friction, have only a secondary 
influence. Changes in properties normal to the streamline are con- 
sidered negligible compared to changes in the direction of flow. Hence, 
average properties over the cross-section are evaluated. Errors are 
minimized when the change in cross-sectional area per unit length is 
gmall, when the radius of curvature is large compared to the diameter, 
and wher the shapes of the velocity and temperature profiles remain 


constent with lengtn.® 
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Assumptions sixteen through eighteen represent ideal, Limiting 
approximations for unproved physical and chemical phenomena, The 
questions of gas solubility in condensed species and of the miscibility 
of liquid phases are worthy of consideration, These phenomena could 
bear some significance, particularly for cases in which liquid metal 
oxides constitute a relatively large percentage of the total mass of 
combustion products. However, assumption eighteen is probably a greater 
potential source of error. The degree of thermal equilibrium between 
the gas and condensed phases is dependent upon the rate of heat transfer, 
while the degree of acceleration of the condensables by the gas is 


a rata ~ 
Teme Vv 


£ moms 
greatly affected by the particle sizes of the condensed phases, Ags the 
particle size of the condensable decreases, both thermal and momentum 
equilibrium are approached to a greater degree. Efforts have been 
made to predict the degree of approach to these equilibria by the use 
of empirical rate constants, /*8 The tentative conclusion of these 
studies was that both thermal and momentum equilibrium are very nearly 
approached in two-phase flow provided the particle diameter of the non- 
gaseous phase does not exceed ten microns, 

The expansion process through the converging-diverging nozzle may 


be analyzed by four basic approaches, These are: 


(1) Shifting chemical equilibrium accompenied by complete 
phase equilibrium. 


(2) Shifting chemical equilibrium with frozen phase equilibrium. 


(3) Frozen chemical equilibrium with complete phase equili- 
brium, i.e. equilibrium condensation, 


(4) Frozen chemical equilibrium and frozen phase equilibriun. 
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The first basis is the most ideal one and establishes the maximm 
theoretical performance of a given propellant combination at an 
established set of operating conditions, On the other hand, the fourth 
basis represents the lower limit of performance. Experimental evidence 
indicates that the equilibrium conditions of the chamber are shifted to 
a certain degree during the expansion process. Only ths degree of the 
shift is questionable, and this undoubtedly is dependent upon the 
particular propellant combination. Theréfore, for practical purposes, 
if both the first and fourth approaches are evaluated, the theoretical 
performance may be "bracketed" by establishing its maximum and minimun 


limits. Sew Oke 


3. Uncertainties in Performance Calculations 
Practically all of the fundsmental uncertainties related to 


theoretical performance calculations are concerned with the phenomenon 
of two-phase flow. Until more basic information is obtained on this 
phenumenon, theoretical performance calculations cannot be refined 
further. The particular problem areas related to rocket propulsion 
which are in need of more basic study include: 


(1) The extent of approach to thermal equilibrium between 
the gaseous and condensed phases, 


(2) The extent to which the velocity of the condensed phases 
lags behind that of the gas phase. 


(3) The rate of formation and growth of the condensed particles. 

(4) The effect of miscibility of two or more condensed phases 
and the effect of solubility of the gas phase constituents 
in the condensed phase, 


(5) The degree of non-ideality of the vapors of condensed phases, 
particularly in the region of saturation conditions. 
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C. Thermodynamic Significance of Specific Impulse 


l. Fundamantal Thermodynamic Concepts 
Thermodynamics is the science which deals with the relationships 
between various forms of energy. Applied to chemical reactions, thermo- 
dynamics not only reveals whether or not a reaction is possible but 
also specifies the quantitative change in energy which occurs during 
the reaction. Thermodynamics is independent of the path of a chemical 
reaction; only the initial and final states need be considered. However, 
f thermodynamics has its limitations. Since it is not concerned with 
i time, thermodynamics will yield no information concerning the rate and 


mechanism of a reaction. 


The first law of thermodynamics is the law of conservation of 
energy. This states that the total energy entering a process, plus any 


addition during the process, must equal the total energy of the system 


leaving the process, i.e. 


Tig I ee ee 


dB = dq-dw (10) 
i, where B = internal energy of the system 
q@ = energy absorbed by the system from 


its surroundings 
w = work done by the system on the 
surroundings 
The sscond law of thermodynamics may be considered as a measure of the 
degree to which energy can be transformed into useful work. The point 
property waich measures the amount by which the energy of a system is 


degraded in going from one state to another is called entropy, and is 


defined by the second law as 
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a 
= 2 L 
ola (11) 


where S = entropy of the system 
T = temperature of the system 
4... reversible energy change at 


constant temperature 
Interpreted statistically, the second law states that, on the average, 
the entropy of an isolated aystem increases continuously until that 
system reaches a state of maximum probability. In reality, all spon- 
taneous reactions are to some extent irreversible and are accompanied 
by a degradation of energy. Finally, the third law of thermodynamics 


entropy of any homogeneous substance which is 


a 


states that the absolut 
in complete internal equilibrium is zero at O°K, The first two laws are 
explicitly involvad in the analysis of rocket propulsion, while the 
third law is utilized implicitly. 


Utilizing the first law, enthalpy is defined as 


aH = dE + a(pV) (12) 
where H = enthalpy 

B = internal energy 

P = pressure 

V = volume 


Since enthalpy is defined in terme of functions of state, it ia also a 


function of the state of a system, By re-writing the first law in the 


aq (38), aT + | (3), + p | av (13) 


the specific neat capacities are defined as 


form 


y= (3). (14) 
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cp = | Maen) = (3F) 05) 


for constant volume and constant pressure, respectively. for the 


general case, 


er ee [> $F], (3) (16) 


B 


but for ideal gases this relationship reduces to cp - cy = R. (17) 
In order to express the first and second laws explicitly in terms 


of temperature and pressure, the Gibb's free energy is defined as 
ar = aH - d(TS) (18) 
or QF = dB + pdV + Vdp - TdS - saT (19) 


Under reversible conditions in which all work is accomplished by 


expansion, 

d& = TdS - paV. (20) 
CGanbining these equations, 

aF = Vdp - Sat. (21) 


This is probably the most useful combined statement of the first and 
second laws since pressure and temperature are the independent variables. 


The concept of free energy is extremely useful as a criterion of chemical 





equilibrium, and is related to the equilibrium constant by the equation 


aF = RT d ink. (22) 


At equilibrium, the free energy change of a chemical reaction ia zero. 
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2. Application of Thermodynamics to Rocket Propulsion 





All of these general thermodynamic relationships are employed in 
the analysis of the combustion and expansion processes of rocket propul- 
sion. However, the particular thermodynamic significance of spscific 
impulse is derived from the law of conservation of energy. For the 
macroscopic, steady-state flow process, Bernoulli's general energy 
equation per unit of mass becomes 


| se +22 + 3+ oats =| eta + + PyVyjea- W (23) 


where EB = internal energy of the fluid 
2 = potential energy of position 
vy = fivid veiocity 
pV = work of expansion 
q = energy absorbed by the system 
from the surroundings 
Wy = non-expansion work done by the system 


Position 1 represents the nozzle exhaust conditions while position 2 
represents the combustion chamber conditions. Since the velocity is 
zero at stagnation conditions, vo drops from the equation. Also, the 
difference in potential energy between the chember and the nozzle exit 


is negligible. Thus, the equation reduces to 


2 
Ar+ A(pv) - 2 = a-W, (24) 
2g 
or 2 
Aq - 2 = q- (25) 
2g 


AH = oh (26) 


al 





or v= V2e AH (27) 


For the case of optimum expansion, 


Vv ates i 
I = & = 2r \H 2. 
spe \/2 A (28) 





where Isp = specific impulse 
Ve = exhaust velocity 
@ = acceleration of gravity 
J = mechanical equivalent of heat 
AH = enthalpy change 


Thus, specific impulse may be defined rigorously from the total enthalpy 
change of the system between the combustion chamber and the nozzie exit. 


Numerically evaluating the specific impulse equation in English units, 


Isp = 6.940 V/ H, ~ Hy (29) 
where He enthalpy of chamber gases, Btu/lb 

He = enthalpy of exhaust gases, Btu/1b 

Isp = ld thrust/1lb mass/second 

J = 778 ft-1b/Btu 

B = 32.2 ft/sec 


Similarly, evaluating the equation in metric units, 


Isp = 9.330 Be - He (30) 
where H, = enthalpy of chamber gases, 
calories/gram 
Hy = enthalpy of exhaust gases, 
calories/gram 


Isp = lb thrust/lb mass/second 
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These equations are rigorous and apply to isentropic two-phase flow as 
well as to isentropic compressible flow. In the latter case, another 


form of the specific impulse equation may be utilized. Since (32) =O 








T 
for an ideal gas, the adiabatic change in enthalpy is 
H = feo at (31) 4 
Since en = gr R (32) 
Au= By * fy “fo - Te). (33) 
Under isentropic conditions, : 
qi 
- f pe] Set t 
T= T, [#| k (34) 
RT Pp \ Kel 
~&: . %{[a-(22) x (35) 
ss oe Be Be ft (F) 
Therefore, (36) 
le k=] 
from which Isp = \/2R - KK . Fo [ L- Pe \ k ] (37) 
c Br wt We) | 


" 


where combustion or stagnation temperature 


mean molecular weight of combustion products 
mean ratio of the specific heat capacities 
of the exhaust gas, ¢p/cy 

chamber or stagnation pressure 

exhaust pressure 

specific impulse 

gas constant 

acceleration of gravity 
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when K = 1544 ft-lb/(mole)(°R), g = 52.2 ft/sec*, and Ty = OR, 


Isp = 9-797 (38) 





D. Steady, One-Dimensional, Isentropic, Compressibie Flow 


l. General Principles 
An isentropic process is both reversible and adiabatic and may 

be considered as the limit to be approached by real adiabatic processes. 

as was discussed earlier, this approximation constitutes a limiting 

condition for expansion through a converging-diverging rocket nozzle. 

Compressib)- flow implies variations in density throughout the field of 
flow, resulting primarily ‘rom changes in pressure. The rate of change 
of density with respect to pressure is an important variable in the 
analysis of compressible flow and is inberently related to the velocity 


of propagation of 4 pressure disturbance, !° The one-dimensional flow 





approximation implies that the rate of change of fluid properties normal 


to the streamliue is negligible compared to the rate of change along the 


streamline. If condensed phases are present during the expansion 


process, the concept of compressibility discussed in this section mst 





be modified. Nevertheless, the principles which follow constitute a 





basic, useful approach to the analysis of rocket propulsion. 
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2. Velocity of Pressure Wave Propagation 


For compressible flow, the Navier - Stokes! solution to the general 


equations of change may be written agit 


Bee - We) + p[VWV)] + b [Vivevow | + EAR (39) 


ta the absence of external and viscous forces, the Navier-Stokes’ 


solution reduces to 
222 + 9 (40) 
for one-dimensional flow. Similarly, the general equation of continuity 
for the fluid as a whole is 

Be = - (VV) (41) 


which reduces to 


ay. . 2 3 i 
2 PSE + = H+ St Gu) FO (lz) 
for one-dimensional flow. Under steady-state conditions, these 


equations become” 


‘ 
Mee bw tw) 
nik Pes = - vz of (4h) 


Simultensous solution of the above equations defines the velocity of 
propagation of a pressure disturbance, or the velocity of sound, for 
one-dimensional, steady-state, ideal flow. As a limit, this velocity 


of propagation may be considered as isentropic for an infinitesimal 
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pressure distrubance, end expressed mathematic-lly as 


Ve =O = V (2p /2F), (45) 


For the isentropic flow of a perfect gas, p/p = RT and r/p* is 


a constant. Therefore, 


a . ag 6 
— = k (46) 
2p = 2B = xr, 7) - 
and ( ) k kRT (47) 


The velocity of sound may now be expressed as 


Ve = Ver (48) 


pressure 

density 

mean specific heat capacity 
ratio, Cp/cy 


universal gas constant 
absolute temperature 
velocity of sound 


a 


where 


MD 
phon 


eww 
nun 


3. Basic Flow Equations 


Velocity, temperature, pressure, and density are the four i-de- 
pendent properties selected to describes the state of the system. 
Therefore, four independent equations are required to express the 
relationship between these properties at any two points in the nozzle. 
These equations are: (1) the equation of conservation of mass; (2) 
the equation of conservation of momentum; (3) the equation of conserva- 


tion of energy; and (4) the equation of state. 





Under steady-state conditions, the mass rate of flow through a 


differential length of nozzle is 


Ape = apr (49) 


where cross-sectional area 


fluid density 

fluid velocity 

distance along nozzle axis 
time 


if 
+R OND 


town u 


Therefore, the equation of conservation of mass is 


$ (pr) 20 (50) 


ot San AR a TSS 


te 
i 


Ra 


Which, upon integration, becomes 


= : = PY (51) 
i. 

: or m = A 

: v (52) 
2 

z where m = mass rate of flow 

* V' = specific volume of fluid 

ny 


iae rate of change of momentum at a particular inlet position is 


ape SE = ape 63) 


while the rate of change of momeutum at any downstream point is 


tee ee La 
it 


apv2 + 9. (apy? )ax (5h) 


From Newton's laws of motion, the difference in these rates of change 





of momentum equals the resultant of the forces acting on the inclosed 
volume of fluid. These forces are an upstream force Ap on the inlet 


section, a downstream force Ap + S- (ap )ax on the exit section, and a 


2T 





force Ppodx exerted oy the walls of the nozzle. In this latter case, 
P is the perimeter of the section and &is the half-~angle of divergence 
of the nogzle, Thus, 


-& (a,)ax + PyMdx = 9 (Apv)ax (55) 


However, dA = Pdx tane( or dA = P@dx since his small. Therefore, 


the equation of conservation of momentum may be written 


- © @) + e#@ = & wp) (56) 
or a 2 + o (a pv) = 0 (57) 


The total energy contained by the fluid stream consists of its 
internal energy and its kinetic energy. Using the previously developed 


nomenclature, the rate in influx of energy is 
APY (B+ hv") (58) 
while the rate of efflux of energy at a downstream position is 
apy (a+ hv) + 2 [ aprige vy] ax G9) 


However, external forces are doing work on the gas at a rate Apv on the 
inlet surface and at a rate Apv + ed (Apv)dx on the outlet surface. 


Therefore, the total energy balance through the nozzle becomes 


2 
- = (Apv)dx = a [apr (E thy ) [ex (60) 
or | 4ee wal: + 1 ¥)| = 0. (61) 
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In terms of enthalpy, the equation of conservation of energy is 


co te 
| Aer ee 3%) | = 0. (62) 
The equation-of-state for an ideal gas is B/p = RT, Since 


a& = oyGT and cp - cy = R, 


a (B) = (ep - cy) = | (63) 


or ag = gh a (3) (64) 


Upon integration, 


3 (65) 


4. Significance of the Mach Number 

The equation of conservation of mass, the equation of conservation 
of momentum, and the Mach number are each independent of either the 
adiabatic or the isentropic restriction. This latter term, the Mach 
number, is a dimensionless quantity obtained from inspection analysis 
ef the Navier-Stokes equations, and is defined as the ratio of stream 


velocity to the local velocity of sound, i.e. 


M= (66) 


al 


The Mach number has a profound influence upon the nature of compressible 
flow, as saown by the following rigorous mathematical derivation. 
By combining the continuity equation and the law of conservation 


of energy, ths Bernoulli equation is obtained in differential form as 
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2 
dH = -a (2s -v dv 
From thermodynamics, 
Tas = dH - Sp 
‘i 


which, at constant entropy, reduces to 


aH = ap 
P 
Thus, a, = -pv dv 


(67) 


(68) 


(69) 


(70) 


This is Buler's equaticn of motion for an inviscid fluid, showing that 


pressure decreases with acceleraticen of flow and increases with deceleration 


of flow. Expressing the equation of continuity in logarithmic 


aifferential form, 


dinp+tdlnaAtdlav =0 


or f+ ah + Lg 
Pp A v 7 


Substituting and rearranging, 
— is | fad = gf 


ae Yl. iw 
or ws = [1 - m+ bia. 


Thus, the following conclusions are reached : 19 
(1) For subsonic flow (M<1), 


aa 


a a 
a>? ay <0 
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(2) For sonic flow (M = 1), 


da = 
ane (76) 


oi 


(3) For supersonic flow (M>1), 


i <0 & > o. (77) 


From these relationships, it is seen that the effect of an area change 
is exactly the opposite for subsonic and supersonic flow. This leads 
to the interesting conclusion that an accelerating stream starting 
from rest must first decrease in cross-sectional aréa uati. s velocity 
of Mach = 1 is attained. Then, the same stream must increase in cross- 
Sectional area in order to reach supersonic velocities. In practice, 


this phenomenon is achieved with the converging-diverging nozzle. 


5» Adiabatic Flow Relations 
Since all possidle states of isentropic flow lie on a line of 


constant entropy, one cf these states must correspond to zero velocity. 
The conditions at this point, referred to as the stagnation conditions, 
are independent of any changes which may occur downstream. The stag~ 
nation enthalpy, temperatvre, and pressure are designated Ho, Tp, and 
Ro respectively. Although the expansion process through the nozzle is 
isentropic, several useful flow restrictions are dependent only upon 
the adiabatic restriction. From the first law of thermodynamics, the 
total energy between the stagnation point and dny downstream position 


must be conserved, i.e. 


Ho = H+ o (78) 


iaascas ie iinet wide. 


Also, under ideal conditions, 


AH = Cp AT 


Cp - Cy = R 


ep/cy = k 


= 
o * ger 8 


(79) 


(80) 


(81) 


(82) 


Thus, the velocity at any position may bs expressed in terms of the 


stagnation temperature as 


v= V\2ep (ToT) = iE R (T,-T) 


(83) 


The ratio of the stagnation temperature to any adiabatic downstream 


temperature may be expressed in terms of Mach number and specific heat 


ratio, Again by the first law of thermodynamics, 


t 2 2 
—_— vs ° i 
F =s i * 268 1 + ie ep 


= 2. 
However, ep = pH R and c© = kRT 
Therefore, az = 1 +E s 
ry 
or To = 1 + Kolb w* 
T 2 


Finally, the flow per unit area for any adiabatic process may be 


expressed as 


>|s 
" 
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v= BH 
ir “Vert 
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6. Isentropic Flow Relations 





Under the conditions of isentropic flow, the pressure-volume- 


temperature relationship for ideal gases becomes 
k=] 
Ts Po sete Vy \ kel 
= = ok (89) 
a -(# a) 


Utilizing this, the pressure and density at any downstream position 


may be evaluated from the equations 


fel 
Bo =[.+ | (90) 
oo 
and 5 -[2 + “|e (92) 


Furthermore, the flow per unit area of cross-section now becomes 





m =) k . Do ° x (92) 
A VR ip 
Ro [2 +h | ee 


By differentiating this equation with respect to M and setting the 
derivative equal to zero, it is found that the maximum m/A exists at 
the throat where M=1. 14 Thus, the maximum flow per unit area of cross- 


section is 


(93) 
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where A* denotes the throst area, or critical area, at M= 1. An area 


ratia between the throat and any downstream position may now be expressed as 


kts 
42 mad (Sr) are (94) 


For each equation containing the term we » theze is both a subsonic 
solution and a supersonic solution. I+ has besn shown that the sub- 
sonic solution applies to the converging portion of the nozzle while 
the supersonic solution is valid in the diverging portion. Correspond- 
ing to any given value of Mach number, there is one precise value of 
P/Poo T/Egs (/ Pos and &/A* which can be evaluated fron equations 
already discussed. Shapiro has tabulated a set of working tables in 
which these dimensionless ratios are summarized for an ideal gas 


having k = lke 


Te Choking 

The fact that the masa flow rate per unit area possesses a 
maximum value at the throat leads to a phencmenon called chokinge 
For a given converging-diverging nozzle operation under fixed atagnation 
coaditions, the nature of the flow depends upon the pressure retio Pe/Poe 
if the external downstream pressure is gradually reduced from an initial 
condition of D,=P,; the velocity through the converging section of the 
nozzle increases uitil M =1 at the throat. This is the critical 
pressure ratio. Further reduction in the external pressure will result 
in supersonic flow being maintained beyond the throat. However, true 


isentropic flow beyond the throat can exist only if the external 
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pressure ig equal to or less than the one specific ideel exhaust pres- 
sure defined by the nozzle area ratio, If the external pressure 
exceeds the ideal exhaust pressure but is less than the critical throat 
pressure, irreversible discontinuities occur in the diverging section 
of the nozzle. Although these discontinuities, called normal shock 
waves, may occur adiabatically, the entropy of the fluid is always 
increased. The normai discontinuity always involves a change from 
supersonic flow to subsonic flow accompanied by a pressure rise and an 
energy loss. The magnitude of this energy loss may be evaluated from 
the Fanno and Rayleigh curves for the particular Sxpeast oh prodsiess 
The Fanno line is the locus of points which satisfies the energy 
equation, continuity equation, and the equation of state. On the other 
hand, the Rayleigh line is the locus of points satisfying the momentum 
equation, the continuity equation, and the equation of state. A 
straight line connecting the two points of intersection of these loci 
represents the normal shock wave. Shapiro has also developed a set 
of working tables which relate the pressure ratio, density ratio, 
temperature ratio, and area ratio across the shock wave to the upstream 
and downstream Mach numbers. 1 

If the nozzle exit pressure exceeds the external pressure, the 
fluid will contiaue to expand after being exhausted from the nozzle. 
However, no useful work is derivec from this free expansion. Three- 
dimensional turbulence occurs in vhe immediate region outside the 
nozzle. Again a loss of energy is incurred. Thus, either overexpansion 
or underexpansion of the working fluid results in reduced performance. 
Ideally, the nozzle exhaust pressure should equal the external downstream 


pressure. 
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8, Isentropic Flow Performance Parameters 


By introducing the isentropic flow relationships into the velocity 


eyuations already developed, the exhaust velocity is evaluated as 


Ve = 2k Rte To l Ps\ oa (95) 
e* er 7S iw ~\ Bo oe 


where MW 
R? 


molecular weight 
gas constant expressed in work units 


The maximum possible exhaust velocity is obtained at an infinite 


pressure ratio po/pe, for which the ideal cycle efficiency 


becomes unity. When po/pe becomes infinite, 
 eormanes ih 
Vmax = \per ar: a (96) 
Simultaneously, T/To and (/ fo become zero while A/A* becomes infinite. 
This leads to the interesting conclusion that the maximum attainable 
exhaust velocity is finite even though the working fluid is expanded 
to an infinitely low pressure, temperature, and density. Earlier, the 


exhaust velocity was shown to be related to specific impulse by the 


equation 


Vo = Isp ’gs 


for ideal expansion. 
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Although the characteristic velocity c* is readily evaluated from 
experimental data, its theoretical evaluation requires a knowledge of 
throat conditions. By combining the definition of ¢* and the equation 


of continuity, the theoretical relationship becomes 





ot = —2e (97) 


The conditions at the throat can be expressed in terms of the stagnation 


conditions as 


Be 7 2.\e 

Bo (2x) ¥ E (98) 
2 

T - 2. «(oc (99> 

T> k+] So 

e of2\4 

Po oT k=l (100) 


where k is the average specific heat ratio between stagnation and 

throat conditions. However, an accurate determination of ¢* is diffi- 
cult because of the change of k with chemical reaction. Probably the 
most rigorous approach is to maximize fy by calculating this quantity 


at several assumed throat pressures. 


Since ve = ct = Qn. 
(3) 8 (101) 
and p*- MWe = ponte (102) 
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2 
es = Po , BI* B =) ue Flo ue 
p* kemwe pe \T, kMW* (103) 
ct = 
(104) 





For the case of a nonereacting gas with constant specific heat, this 
equation is readily evaluated, Otherwise, a trial-and-error solution 
is required. In general, c* may be considered as being affected by 
conditions upstream of the throat but as being independent of downstream 
conditions. 


Theoretical evaluation of the thrust coeffielent is also somewhat 


involed although its experimental evaluation is straightforward: 


Cp = 1 + cosh + (pe- (105) 
Pea®* 
= ¥. =pa ° Ag 6 
or Cp st ag re (106) 


‘Unlike c*, cp is a function of exit pressure and of nozzle gecmetry 


downstream of the throat. For ideal expansion, i.e. when pe=pa, 


re s Ye 
FO" 38 (107) 


Expressing this relationship in terms of the flow equations for a 


chemically-reacting gas with variable specific heat, 














average specific heat ratio 

during expansion 

= average specific heat ratio between 
the chamber and the throat 

k* = average specific heat ratio at the 

throat 


where 


w OF! 
" 


The thrust imparted to a missile may be expressed in terms of the 


thrust coefficient as 
F? = A* py cp + (pe-pa)Ae (109) 


This applies either theoretically or experimentally, Finally, the area 


" yatio may be expressed in terms of the pressure ratio and the specific 


heat ratio as 





L 
72 k 
4. = (110) 
er ae ee Se 
(=) ay Ve | re -(2) Kk 
k-1 Po 
or in terms of the thrust coefficient as 
ke) L 
A = bt 2. | Kl | Po (111) 
a* Gg | kel Pp 
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E. Steady, One-Dimensional, Isentropic, Multi-Phase Flow 


1. General Principles 


The previous analysis of isentropic, compressible flow is completely 


rigorous only if all the products of combustion are ideal gases. 
Theoretical analysis of the expansion process becomes considerably more 
complicated if a condensed phase is present during all or part of the 
expansion. Either solids, liquids, or both may be present in the prod- 
ucts of combustion, depending upon the composition of the propellant 
and the operating conditions, 

. In any flow system consisting of more than one phase, the general 
equations of change apply in each phase. Furthermore, the solutions to 
the equations of change in the individual phases must be matched at the 
interfaces. In view of the fact that the mechanisms of the various 
interphase transfers have not been completely defined, a general milti- 
phase transport solution to the general equations of change has not 
been obtainea, +6 The problems of two-phase momentum transfer alone are 
quite complex and not adaptable to a simple, rigorous analytical 
solution. The simultaneous occurrence of heat and mass transfer adds 
to the complexity of the problem. Therefore, discussion of multi- 
Phase flow in a converging-diverging nozzle wili be limited to a 
seiacqumat tative aeeatment: Using this approach, the degree of ettain~ 
ment of thermal and momentum equilibrium between the gas phase and the 
condensed phases can be analyzed. 

By applying the ussumptions which were presented and discussed 
previously, thermodynamic relations can be formulated for: (1) the 
case of a condensed phase existing in constant concentration throughout 


ho 


the expansion; and (2) the case of condensation during expansion. In 
the latter case, the concentration of the condensed phase is a function 
of the temperature, pressure, and rate of condensation. Whereas for 


isentropic flow of a perfect gas, 


pv 


= constant (112) 


the corresponding relationship for isentropic, multiphase flow ist? 


pv? 


= constant (113) 
&)R + @jSy_ + Miace 
3 ¥, * Plan ¥ 
where $ = ae : ; “2 uy) 


es + MWa2cps 


Pp = pressure 
Y= volume 
CV, = constant volume specific heat capacity 


of the gas phase 

constant pressure specific heat 
capacity of the condensed phase 

mean molecular weight of the gas phase 
mass fraction of gas phase 


mass fraction of condensed phase 


Cpg/ev, 
gas constant (cpg-Cy,) 


u 
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This latter equation includes the specific heat capacity of the condensed 
phase waile retaining the form of the equation for the isentropic flow 


of gases. &valuating the relationship between k and 9d, 


hi 
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when a = 0, g=—# = k 
Cy 
& 
when a > 0, k> @>1 


2. Flow with Constant Concentration of Condensables 





Yor the case in which a condensed phase exists in constant 


concentration during the expansion process, the enthalpy balance for 


adiabatic, frictionless flow igi 
aH = c, aT, + £cp, aT, = Vp a, = Rr, 2 (115 ) 
Bg “@ * “Pa 's & “Pp Ep 
where H enthalpy of multiphase system 


Spy specific heat cf gas cer unit mass 

Cp, = specific heat of condensed phase 
per unit mass 

f = mass ratio of solid to gas 

Tp = temperature of gas 

Ts = temperature of condensed phase 

p> = pressure 

Vg = specific volume of gas 

R = R/MW 


Integrating, 


He - Hy = Op (To-Teg) + fepy (Te-Tog) = L (v,? + fv") (116) 


where ec denotes chamber conditions 
e denotes exhaust conditions 


The relatio’.ship between temperature, prescure, and heat capacity is 


at, | at, _ =a 
Dg +f CDs 8 i = R 5? (117) 
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: aT : 
Assuming that —8% varies linearly with pressure, the integrated form 


ef the above equation becomes 


7 Cp + fey dl, 
=e -[2e e * at, (118) 
c ¢c 


For the limiting condition of Tg = tT =T, 


Te 5 Ea ee 





Ty Po (119) 
while for the other extreme of Tg = Tos ats = 0 
and 
R 
e 
Pp, : 
Te -[22| & (126) 
Ty Po 
For the general case, 
° 2 
" . at, } ee es ae dee 
| om * f ope ate | (Fe - Tog) = 2 tet ye (221) 


3. Flow with Simultaneous Condensation 


If condensation occurs during the expansion, the flow conditions 
are altered considerably since the phase change from gas to liquid or 
solid is accompanied by the release of heat and a decrease in the 
quantity of working fluid. assuming thermal and kinetic equilibrium 


for the condensed phase, the enthalpy balance ial? 








aH = cp, aT + f cy aT + Aly ang = V'g dp (122) 
where 
Ay = molar latent heat of phase change 
Dg = moles of gas 


Cp = specific heat capacity of the newly 
formed condensed phase (solid or liquid) 
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By applying the ideal gas law, 
£ Any(® at | at =n, noe = Foe (124) 
Eg 2 Soe AY ar | poe SBT AS 3 


This squation can be solved only if ae the rate of condensation or 
nucleation, is known as a function of temperature and pressure, Although 
a considerable amount of work has been done on the rates of nucleation, 
the equations are too complex to be readily solved. Qualitatively, the 
theory shows that the rate of condensation increases with the ratio of 
the actual vaper pressure to the saturation vapor ppeesunen? During 

a flow process, the rate of condensation will always lag somewhat behind 


that value required to maintain equilibrium. 


4. Approach to Momentum Equilibrium 


In order to evaluate the order of magnitude of momentur and heat 
transfer between the gas and condensed phases, the semi-quantitative 
treatment will be extended to these phenomena, The condensed particles 
will be assumed to be spherical. Using this approximation, the force 


or drag acting on a moving spherical body is 


2 
Fy =mta' = $ Cp at P, (vg-vs) (125) 

where Fy = force acting on the particle 

m' = mass of particle 

at = acceleration of particle 

Cy = drag coefficient 

fe = density of gas 

A' = frontal area of body =Thrr* 

Vg = velocity of gas phase 

Vg = velocity of condensed particle 


ae a Rn a nn — 
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The mass of the particle is 
Bs amp x (126) 
3 fg 


from which the -cceleration of the particle is 


CG 2 

ae ae = 2fe"D (vg = vs) (127) 
8 Ps r 

Although a rigorous solution of this equation requires the experimental 

evaluation of the drag coefficient, an approximate treatment may be 


based upon Stokes! Law, The equation of motion for spherical particles 


under the influence of Stokes! drag is 


2 
1 0x = 6Mayr - 128 
aes Mgr (vg - Vg) (128) 
where m* mass of particle 


viscosity of the gas (assumed constant) 
radius of particle 


wouu 


Bg 
r 
This equation cannot be integrated directly for the general case, 


but an exact sol.tion is obtained for linear acceleration of the gas 


stream, i.e. 


1 
Vg = (tg, + 2a'x)e (129) 
By combining these two equations and using the relationship vg = at 
* 


2, ice ea 
ai + OTT Mgr 3% - GT 4x (vg, *@ax)£=0 (130) 


where Ygo = gas velocity at the instant the particle 
under consideration is injected into the 
gas stream 
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4& more useful expression of these relationships is 


Ss . OMe | te 1 (131) 
ax mV g6 & 
2 
Vv x 2a! 
and (#8 = 1+ Say (132) 
Bo 
XN 
where TK = 1 end SB = 0 at x = 0. 
: "go “Bo 


In rapidly accelerating flow, the Reynold's number for the particle may 
readily exceed the limit for Stokes® drag. The equation of motion then 


takes the non-linear form* 


oe - © op x v3 )° (133) 


However, provided the Stokes’ drag coefficient is cq = 24/Re, the 


acceleration is found to be 


avg 


are (vg-va) (134) 





This equation may be solved with the aid of the equation for the conser- 
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were 


vation of energy and the equation for the rate of tempsrature change 


aT/dt. Maxwell repo:ts that when r = 5x10 cm, V¥g=0.98 Vg, after only 


4% of the particle's residence time in the nozzle has elapsed,“3 


5. Approach to Thermal Equilibrium 
Analysis of thermal equilibrium between the gas and condensed 


phases is somewhat more complex because three major ‘heat transfer 





mechanisms are involved. These are convective heat transfer between 
the particies and the gas , heat transfer within the particle, and 


radiation from the particle. The convective heat transfer from the 
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condensed particie io the gas may be represented by 


dg? = ‘ - } 
aa = bat (1-72) (135) 
where h = average heat transfer coefficient 
A’ = effective surface ares 
qg*' = total heat transferred 


Assuming instantaneous heat transfer within the particle, 
’ art 
ga = > Chg ac = hA'(T,-Tp) (136) 


and for spherical particles 


: aT, a h Top) : 
i Pat 8 Tg) (137) 


This shows that the rate of heat transfer increases as the particle 
size decreases, Altman and Garter report that Tg-Tg will be less than 
160°K provided r does not exceed 16 én diameter.*4 Assuming that 
heat transfer is rapid at the surface of the particle, the limitation 
on the rate of transfer will be thermal conductirity within the 
particle. Altman and Carter have also shown that this is rapid for 
particles of less than 10° em diameter. The borderline for adequate 


3 


conduction appears to lie between diameters of 107 ana 10 4 cm. 


The energy emitted from the condensed particle by radiation is 
- @ = anreo tH (138) 


Stefan~-Boltzmann constant 
emissivity constant 
enexgy transferred by radiation 


where a 
€ 
Et 


nittw 
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For spherical particles, 


~a 2. a¢tlo 
at Parte (139) 
Altman and Carter report that very little energy is lost by radiation 
=—2 
from particles which exceed 10 cm diameter. On the other hand, the 


5 


rate of the temperature drop of particles less than 10~ cm diameter 
ean exceed that of the gas during expansion.” Maxwell reports that 
particles not exceeding 5x10 7 cm diameter can decrease in temperature 
as rapidiy as the gas with a temperature difference of 20°K being maine 
tained between phases.-° 
From the previous diseussion, it may be concluded that both 

momentum and thermal equilibrium between the gas phase and the con- 
densed phase are rapidly approached provided the particle diameters 


4 


are less than 10 cm. As the particle size increases, the degree of 


attainment of equilibrium decreases. The following table, prepared by 
Altman and Carter, shows a tabulation of specific impulse values for 


both complete thermal equilibrium and complete lack of thermal equilib- 


‘rium at various assumed particle velocities corresponding to a mass 


fraction of 20% solid:27 






Ys 
Vg 


Isp at Tg=Ty Isp at Tg=T, 





0 203 
0.25 214 
0250 221 
1.00 227 
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From a practical viewpoint, this table shows that the attainment of 
momentum equilibrium has more significance than the attainment of 
thermal equilibrium, Under Altman's assumed set of conditions, a maxi- 
wum Isp change of 1 1/2% can result from non-thermal equilibrium while 


a maximum Isp change of almost 10% can result from non-momentum equilibrium. 


6. Fundamental Equations of Flow 


Assuming the validity of thermal and momentum equilibrium between 
the gas phase and condensed phase, thers is no difference between the 
proverties of the multiphase mixture and the ideal gas considered in 
the previous section on compressible flow. The general equations of 
conservation of mass, momentum, and energy remain valid. The equation 
of state is now expressed as 

z= J +2 (140) 
g-l f/f 
where @ has replaced k, Similarly, the isentropic flow relations and 
the isentropic performance parameters for compressible flow may be 
adapted to multiphese flow by repsacing k by the modified specific heat 
ratio }- 

Since pis less than k, the volume of the system becomes more 
sensitive to adiabatic pressure variations when condensed phases are 
present. This result applies only to adiabatic pressure changes when 
both compressibility and heat transfer are involved. As the concen- 
tration of the condensed phase is increased, the nozzle discharge 


coefficient increaseS, This is to be expected since the condensed 


4g 











phase possesses a much greater density. Hence, the condensed phase is more 
readily transported through the nozzle throat where the rate of flow is 
limited by choking. The presence of a condensed phase has a relatively 
minor effect on the nozzle thrust coefficient. In fact, this i3 also as 
expected since the thrust coefficient is comparatively insensitive to 
Variations in the adiabatic exponent k for pure gases, The presence of 6 
condensed phase necessitates a samewhat larger nozzle area ratio in order 
to achieve optimum reaults. From a practical atandpoint, phase separation 
ids more likely to occur as the concentration of condensablea is increased. 


Therefore, a smaller angle of divergence is desirable in the nozzle. 


7. Uncorteintios in Milti-Phase Flow 
Earlier it was pointed out that the only fundemental uncertainties 


related to theoretical performance calculations are those concerned with 

the phenomenon of multi-phase flow. Additional fundamental studies ere 
needed on the attainment of momentum equilibrium between phases, particularly 
when the condensed-phase content of the system is 80-90% by weight. Exieting 
data are confined to the range of 20-35% by weight soniensables. Conceivably, 
the energy lowses dus tc drag eould be more appreciable at the higher 
condensed particle losding. Also, upon expanding such syatems to sub- 
atmospheric pressures, phase separation in the nozzle may be n more severe 
problem. Another fundemental uncertainty of potential significance involves 
the rate of formation and growth of condensed particles. This factor . 
affects not only the size of the condensed particles but also the attainment 
cf chemical and phase equilibrium during combustion and expansion. Ths 
questions of miscibility between condensed species, gas solubility in 


condensed phases, and non-ideality of condensable vapors are, in all 


probability of lesser consequence, 
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F. Ballistic Importance of Specific Impulse 


From the previously developed equations, it has been shown that 
specific impulse rigorously defines the amount of useful energy which 
may be obtained from the combustion and expansion of a fuel and oxidizer. 
AS an approximation, specific impulse is directly proportional to the 
squere root of the combustion temperature divided by the mean molecular 
weight of the exhaust products. The effect of the specific heat ratio 
upon specific impulse is not as significant as the effect of the com- 


bustion temperature and the mean molecular weight of the exhaust products. 





In addition to its thermodynamic significance, specific impulse 


also possesses an inherent ballistic importance, The burnout velocity, 





the vertical range, the horizontal range, the payload weight, and the 
physical size of a missile ere all influenced by the specific impulse 
of the propellent. Tormey lists the general relationships between 


specific impulse and these various vehicle performance characteristics 


as 28 
; W 
c = Ispege 
(1) Burnout velocity spegeln nd : (141) 
(2) Vertical range = bp + K Isp in EB (142) 
f 

= K’ Wy 
(3) Horizontal range = 2Spo + K Tap ol | (143) 
(4) Payload = f (Isp)(Wy)(W,) (144) 
(5) Physical size = f (Isp)( / bulic) (145) 

1 





where Isp = specific impulse, seconds 
Wy = total initial mass of the vehicle and 
its propellant 

Wp = total final mass of vehicle after burnout 
g = acceleration of gravity 
Dpo = height of missile at burncut 
Spo = horizontal @istance of the missile from 

; its launching point at burnout 
K = physical constant 
ln = natural logerithm to the base e 


These five equations are necessarily quite general and are presented to 
convey the concept rather than the quantitative relationship between 
specific impulse and the performance parameters. In practice, exact 
values for range, payload, size, and vertical height must be determined 
by a detailed analysis of the powerplant and the rocket vehicle. These 
values are dependent upon a series of arbitrary design choices, e.g. 
thrust-to-weight ratio, staging, trajectory, stc. Nevertheless, the 
preceding equations point out the overall importance of specific 


impulse in rocketry and that increases in specific impulse will Lead 


to iacreases in missile performance. The first equation of the series, 


which relates specific impulse to the burnout velocity of the missile, 


is examined in more detail in the next section. 


G. Relationship between Specific Impulse and Propellant Density 


1. Introduction 
Specific impulse is the criterium generally used for comparing 


the relative capabilities of propellant combinations. However, in 


applying rocket propellants te a specific unit of hardware, the propellant 


density must also be considered in order to optimize performance. 
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Both of these parameters influence the incremental boost velocity 


imparted to a missile by the combustion of a propellant. The actual 
relationship between specific impulse and propellant density is a 
complicated one, particularly for miltistage missiles, and depends 
upon the specific missile application. Such a detailed analysis is 
beyond both the theoretical and calculational scope of this report. 
However, a simplified theoretical discussion is presented to point out 
the roles of both parameters in the application of propellant systems 


to actual missiles. 


2. The Masseto-Volums Retic Approach 


The combined effect of specific impulse and deasity on the 
incremental boost velocity imparted to a missile is obtained by applying 


Newton's laws of motion to the missile 330 


é ° a U / 2 
ma =P = ~Isp’g = + geae(pe-De) - Wmrg' sinB- uy = (146) 


where Wa total mass of missile 
Y%y = velocity of the missile 
t = time 
Isp = propellant specific impulse 
& = acceleration of gravity (32,17h ft/sec“) 
&' = local acceleration of gravity 
Ae = nozzle exit area 
Re = nozzle exhaust pressure 
Ps, external surrounding pressure 
B angle of trajectory with the 
horizontal 
y = empirical coefficient dependent upon 


physical conditions of the system 
and surroundings, e.g, missile surface 
area, density of surrounding medium,etc. 


23 





This generalized equation includes the effects of pressure thrust, 
gravity, and aorodynamic drag, in addition to momentum thrust. However, 
for the ideal case, only the momentum thrust must be considered. Then, 


the previous equation reduces to 


av aw, 
a aE = - Tepte (147) 
W, 
Integrating, = Isp*g*ln & = tep’g’ln =D (148) 
Si We Wn-Wp 


s where vp = incremental increagze in missile velocity 
: Wy = total mass of missile (hardware and propellant) 
‘ %p = macs of propellant 
; We = mass of missile hardware (payload and inert arts) 
ie or ¥ 
Vp = Isptgrin| 1 + 2 
i (149) 
! Since Wp = YP (158) 
y, 
vy, = Isp*g*in | 1 + ie (151) 


or ¥p = Ieptg?la [ + i | (452) 


P= density of the propellant 
Vp = volume of propellant in missile 


We = mass of hardware 
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This equation enables the ideal boost velocity to be expressed in 





terms of the propellant density and specific impulse for a giver mass- 
to-volume ratio, i.e. for a given mass of hardware per volume of propellant, 


This, of course, represents a simplified situation, In practice, a 





change of propellant usually is accompanied by a change in hardware 
design. Nevertheless, the previous equation shows that the incremental 
boost velocity becomes independent of propellant density when the mass- 
to-volume ratio becomes zero. On the other hand, density becomes 
equally significant with specific impuise only at an infinite mass-to- 
volume ratio. Therefore, a change in specific impulse is always more 
significant than a corresponding change in propellant density on over- 


all missile performance. 


r 3» The Relative Effectiveness Approach 


a For thé seme situation in which various propellants are compared 





in a given unit of hardware, the relationship between specific impulse 


and density may also be evaluated by Geckler's approach, 2 Neglecting 


the effect of aerodynamic drag and incremental changes in gravity, 





ao ie ee Cr, i 


2 a! = 2 4, 
Vp + Isp*g In nd Iapeg in me (153) 
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burnout velocity of missile 


specific impulse of propellant 
acceleration of grevity 
total mass of missile 


mass of propellant 
Wp = final mss et burnout (W#,_-W)) 


where ‘pb 


oo 
fo ow 


Differentiating, 
avy = Ieptg a| ln 2 + ln = ° a (Isp°g) (154) 


ov avy, = reprera| ta] a why; + 1n 2 ea (Isp*g) (155) 
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Dividing by vp, QB = Aspe , al #| + hae + In@a/"t). a (Isp * 8) (156) 


av ~ a(I e L * J w 
wi Ma ee eB) tee 


¥or a constant volume of propellant, changes in mass arise only froma 


change in propellant density, so 


dilp = Vad? (158) 
and al fl. wate * arog. f (159) 
ok & - a(iepsa) + tebe = ee (160) 
The term aes . z 


must be unity in order for the density and specific impulse to be 
equally influential on burnout velocity. This term, which may be 


considered a “relative effectiveness coefficient", can be expressed as 


“2 

“p — 1 
In |g 

on 





where Wp/ iy is the propellant mass fraction. Upon numerical evaluation, 


the following relationship is found ;3 
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Propellant Relative Effectiveness 





Mass Fraction of Density to Specific Impulse 
0.95 0,317 
0.90 0.391 
0.85 0.448 
0.8u 0.497 
0.75 0.541 
0.70 0.581 
0.66 0.656 
0.50 0.721 
0.40 - 06783 
0.30 0.839 
0.20 0.896 
0.10 02958 


For example, for a missile with a pronsllant mass fraction of 0,90, a 
chabge in density is 39.1% as significant as a corresponding, percentage 


change in specific impulse. 


4. The Arbitrary Parameter Approach 
& third semierigorous appreach has been presented by L.J. Gordon 


in which an arbitrary parameter is defined in terms of specific impulge 
and density as? 





b = Ispe p® (162) 
where Isp = specific impulse 
Pf = density 
b = arbitrary parameter proportional 
to boost velocity 
8 = empirical constant dependent upon 


missile design conditions 
The parameter b is approximately proportional to the boost velocity 
for a propellant under the migsile design conditions used for computing 


the exponent s, which is always contained within the limits 0 @,H1. 
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Also, the exponent 3 is assumed to be independent of density. In log- 


arithmic form, the previous equation becomes 
Log b = log Isp + 8 log P (162) 


This is the slope-intercept form of the equation for a straight line and 
enables propellants to be compared graphically as a series of straight 
lines for a fixed design condition. The ordinate corresponding to s = 0. 
is the log Isp. 
Previously, the ideal boost velocity was shown to be 
vp = Ispeg in = = tepre in 2 * ahs | (163) 


From this, the exponent s is defined as 


Ba, £ (2) 
isp 


of\ | %\ dP (164) 

P Psp 
For the simplified case of a constant ratio of hardware mass to 
propellant volume, the exponent s may be evaluated as 

\ v. f 

= © 
(165) 
8s phe SE ee an as =. B=) 


V; V; 
[> Bp] BP | 


e 


Bln B 


58 


ast 





where are es ae P (166) 
We 


When Wp/Vp = Q, the exponent s = 0. When We/Vp =00, the exponent s = 1, 
By evaluating s for varicus propellant densities, it can be shown that 
8 possesses a relatively constant value at “ixed design conditions, 
Another fundamental design situation is that of constant hardware 
and propellant weight, with the volume of propellant being a variable. 
In this situation, rigorous evalintion of 


: - £(23) 
. g Yb yg Isp (167) 


yields a cumbersome mathematical expression. However, by assuming 
differential increments of ¥, and ? over a small range, it can be 


shown that the exponent s is mich less than for the previous case of 


era 


ine constant missile volume, Both of these cases have considered single 
stege missiles only. 

Calculations of single stage beost velocities based on constant 
volume of propellant ignore the fact that the stage weight varies with 
density. Neglecting structural features, thrust levela, and other 
variables which may have to be changed if a stage weight is changed, 
Gordon has also examined the effect of changing the propeliant density 
in the second stage of a two-stage missile on the mass ratio of the 
first stage.” It i8 assumed that the two-stage missile has been 
designed to give equal boost velocity increments when both stages are 
loaded with the seme propellant. However, this condition will not 


be maintained if different propellants ere considered in each stage. 
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The total boost velocity may be expressed aa 


W, 
Vy = Yb) + Ybo = (Iep)z e in(B) : + (Isp) g (Be), (168) 


8 
Applying the concept of b = Isp- P e the exponent s must now be defined as 


This implies that s is the rate at which a change in propellant density 
affects the boost velocity relative to the rate at which a change 

in specific impulse affects the boost velocity. A rigorous mathematical 
solution of the previous equation would be quite cumbersome. However, 
by assuming incrementai numberical chan.ss in specific impulse and 
density, Gordon has shown that the resultant value of s ia about one~ 
hel? thet for a single stage vehicle at the same ratio of We/V)-36 
Hence, the effect of propellant density on overall missile boost 
velocity is only about one-half as great in the second stage of a two- 
stage missile as in a single state vehicle having the same fixed design 
criteria, This furtner emphasizes the relative importance of propellant 


specific impulse as the key parameter in overall missile performance. 
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CHAPTER IV | 


THERMODYNAMIC ANALYSIS OF PROPELLANT PERFORMANCE 


In the previous discussion of the fundamental theory of rocket 
propulsion, specific impulse was shown to be measured by the change in 


enthalpy of the combustion products during an isentropic expansion 





through a converging-diverging nozzle. Also, the basic assumptions to 
be employed in performance calculations were presented and discussed as 
were the fundamentals of both compressible and miltiphase nozzle flow. 


in this section, two rigorous methods of conducting performance calcu- 


lations are presented. The basic thermodynamic properties of temperature, 

enthalpy, and entropy of the combustion products are calculated at the 

chamber and at the nozzle exit in each method. By using these properties 

and the chemical composition of the combustion products, other properties, 

such as the moles of gas, the specific heat capacity, the specific heat 

ratio, and the molecular weight of the products may be evaluated. | 
In the first method, thermodynamic equilibrium constants are ; 

utilized to express the concentration of certain chosen species in : : 

terms of the concentrations of other more abundant species. The solu- | 

tion is a trial-and-error one in which the chamber enthalpy balance and i| 


the nozzle entropy balance are conducted as assumed chamber and sxit i 


temperatures, respectively. The second method, which was utilized in 





the calculations presented later in this report, is a more general one 


and is applicable to all equilibrium reactions. Instead of utilizing | 





arbitrary equilibrium constants, solution is obtained by minimizing the 





total chemical potential, or free energy, of the complex mixture of com- 
buation products. The combustion and exhaust temperatures are obtained 
as part of the eclution so that these quantities do not enter into the 
trial-and-error iteration. The accuracy cf both methods is dependent 
only upon the accuracy of the basic thermodynamic data for the react- 
ants and the products of combustion. Sample performance caiculations 


for both methods are included as Appendices. 


A. Method of Thermodynemic Equilibrium Constants 


1. General Approach 
This method may be subdivided into the following steps, each of 


1 





which will be presented and discussed in detail: 


(1) The propellant composition, the combustion pressure, 
and exhaust pressure are selected, 


(2) The products of combustion are qualitatively selected, 
and independent mass balance end chemical equilibria 
equations are written to allow simultaneous solution 
of the syste. 


(3) At an assumed combustion temperature, the casrosition 
of the chamber products is calculated to a predetermined 
degree of accuracy. 


(4) dn enthalpy balance is made at the sesumed canbustion 
temperature by comparing the enthalpy of the chamber 
products to the enthalpy of the propellant. Steps (3) 
and (4) are repeated until the combustion temperature ia 
defined, i.¢. the enthalpy balence is satisfied. 


(5) The entropy of the combustion products is calculated 
at the correct chamber temperature. 


(6) At an assumed exhaust temperature, the composition of 
the exhauat products ia calculated to a predetermined 
degree of accuracy. 
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(7) An entropy balance is mede at the assumed exhaust 
temperature by comparing the entropy of the 
eanbustion products in the chamber and at the nozzle 
exit. Steps (6) and (7) are repeated until the 
exhaust temperature is defined, i.e. the entropy 
balance is satisfied. 


(8) The enthalpy of the exhaust products is calculated 
av the correct exhauss temperature. 


(9) The various other thermodynamic and rocket performance 
parameters are evaluated, 


2. Selectjon of Operating Conditions 
To determine the optimm performance of a given propellant 


system, detailed calculations must be made at several fuel-oxidizer 
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ratios. In practice, the operating and exhaust preésures are dictated 


by the rocket hardware and the mission to be performed. However, since 


uo 
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specific impulse is dependent upon the combustion and exhaust pressures, 


@ stendard basis of comparison must be chosen for theoretical calcu- 


ietions. By tradition most solid propellant performance calculations 
are reported at chamber pressures of lOOOFSIA expanded to 14.7; PSIA. 


Qu the other hand, liquid propellant perfarmance calculations are 
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frequently reported at chamber pressures of 300 PSIA, 500 PSI&, or 600 FSIé 


expanded to 14.7 PSIA. However, recently the desire has been expressed 
to standardize all performance calculations at chamber pressures of 
2G00 PSIA expanded to 14.7 PSIA. This practice has been followed in 
the calculations presented later in this report, particularly in view 


of their basic prototype nature, 








3. Mass Balances and Chemical Equilibria 
The composition of the combustion products is obtained by the 


simultaneous solution of independent equations expressing mass balances 
and chemical equilibria. The total number of equations to be solved 
must equal the number of individual species present in the combustion 
products, One masa balance may be written for each element present in 
the propellant. The remainder of the equations must come from equilib- 
ria relationships. Although the total moles of gaseous products is also 
a variable, this quantity may be obtained from the summation of the 
individual contributions of the various gaseous species. 


For the general chemical reaction 


ak + bBog==cC + aD + OF + fF (170) 


the equilibrium constant may be expressed as 


_2 Soy qb ° up f 
Pa +PB 
e¢,d,e,f 
Dp. Dg oy 
or K, = ORR aE = nme (172) 
n 
ees | 
where &, B,C, etc. = chemical reactants or products 
a, Dd, c, etc. = number of moles or the respective 
components A,B,C, etc, in the 
balanced equation 
Pas PB» Pce otc. = partial pressures of the 
components at the assumed reaction 
temperature 
Dae Dg, Og, etc. = actual molar concentrations at 
the assumed temperature 
K = equilibrium pressure constant 
and Kn = equilibrium concentration constant 








While mass balances are preferably written on a mole basis, thermodynamic 
equilibrium constants are usually expressed in terms of the partial 


pressures or activities of the chemical species. However, use of the 
equilibrium pressure constant enables all unknowns to be expressed in 


mole units. According to Dalton'’s Law, 


af 
?j ° dp PE (173) 


where Ry partial pressure of the j th product 
= moles of the j th product 


nj 
ap = 
2 


total gaseous moles 


= total pressure 


Then, K, © (174) 


(175) 





= «(th) (176) 


where An = the number of woles of products minus the 
number of moles of reactants 


Knowledge of the chemical species existing at the combustion 
temperature is a prerequisite to making accurate performance calculations. 
In general, every compound that can theoretically exist will be present 
in some concentraticn and must be considered until it is proven to exist 


in negligible quantities. For example, in a propellant system containing 
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carbon, hydrogen, oxygen, and nitrogen, at least the following species 


mast be considered as products of combustion:3/ 
HO, Hoe Go, Cb. Oo, No, OH, No, H, anf QO; 


The appropriate mass balances may be expressed as: 


(1) op ByO + BHy + DG0 + O¢9, + Op * BN * Boy 
Dyo + By + Dg 

(2) ay = 2ny, + QnHoS + Ooy + By 

(3) By = 2x9 * Bgo + 2ngo9 + 2Moy + BNO + Nog * No 
(4) ag = 360 * B¢o2 

(5) on = 20yp + ONO 


(177) 


(178) 
(179) 


— (280) 


(181) 


The following chemical equilibria may be chosen end their equilib- 


rium constants expressed as follows: 


(6) C02 + Ho Co + B20 


= 7G0-PH,0 = nG0 et, = a4 
* 2602 ePig zi ®C02-"Hp »| 4 





(7) HpQ + 1/2 Noug Hp + NO 


_ PiaePAO ae 1/2 
* ve ="? [3 
PO". 


(8) 2H,0Gg>— 2H + Op 
1 


2 


(282) 


(183) 


(184) 
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(9) 1/2 Huge (185) ! 
_ 7x - xf 2p’? 
oon eal 
Pip 
(10) H20mge "1/2 fe + OH (186) 
1/2 
re Bip” (? Po 
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Since ten independent equations have been expressed in tea unknowns, 

the product composition can be determined by simultaneous solution. 
The use of snecific chemical equilibria, as shown above, is quite 

desirable for hand calculations. However, this approach is less 

desirable for programming a high speed digital computer since the 

incorporation of a new product of combustion into the study requires a 

modification to the program. To avoid this, a more general equilibria 

approach is to utilize the equilibrium constants for the formation of 


each product of combustion from its monatomic gaseous apeciss. For 


example, 





PaIcl 
” Pal-PoL | 


(187) 





he t Iteration 


Since the previous composition calculation may involve fifty or H 
more simultaneous equations, direct algebraic soiution is not generally 


feasible. Therefore, some type of iterative solution is most conveniently 
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used, depending on the chemistry of the system, Typical methods of 
iteration are discussed in the literature.3%» 39» 40 Iteration is 
continued until the quantity of each product converges within a pre- 


6. 
gram-moles per gram of 


determined degree of accuracy, 6.g2 1,0 x 10” 
propellant.42 In view of these detailed, repetitive procedures, high- 
speed digital computers are profitably utilized for obtaining rapid 
rigorous solutions. In fact, the complexity of the problem almost 
dictates the use, both economically and technically, of a large digital 


computer, such as the IBM 704 or the Remington-Rand 11034. 


5. Bathalpy Balance 


after determining the chamber composition at the assumed combus~ 
tion temperature, an enthalpy balance is made, The enthalpy of the 


products of combustion may be obtained by integrating the energy equation 


te feerg, *(4F) (188) 


along an isotherm between a very low pressure p° and any pressure p at 


which the gas mixture is ideal. This gives4= 
‘ rk 
wee of) [v-2(-34) |o (169) 


and subsequently 


H = Inj H, (19a) 


where H enthalpy of the complex mixture 
nj = moles of the j th product 


Hj 


i 


enthalpy of the j th product 
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The enthalpy of mixing is zero for a mixture of id-1l gases. When the 


enthalpy of the producta of cambustion above a atandard reference temperature 
equals the heat of reaction of the propellant, the adiabatic combustion 
temperature is defined. Mathematically, 
— Enj(Hy), = ABR = Eny (Aves), -Eny (ABs), (193) 

where. H, = enthalpy of the j th component at 298°K 

w =heat of reaction at 298°K 

ny = moles of the j the component 

AE, = heat-of-formation of the j th component at 298K 


and Dp = producta 
r = reactants 


Recanbining the above equati oi, 
t, 


This states that the summation of the heat~-of-formation at the standard 
temperature plus the change in enthalpy above that temperature of all the 
combustion products must equal the heat~-of-formation of the propellant 
at the standard temperature. For each asswaed combustion temperatune, the 
chember composition and the chamber enthalpy balance must be calculated 


as described above, These steps are repeated until the enthalpy balance is 
satisfied. 


6. i Chamber Entropy 
After the combustion temperature, composition, and enthalpy have been 


defined, the entropy of the complex mixture is calculated. Rach individual 





product of combustion contributes to the total entropy of the nixture, as 
in the calculation of total enthalpy. However, the direct mixing of gases 
at constant pressure ani temoeraturs is ar irreversible process, so there is 


an additional entropy-of-mixing contribution. In order to evaluate thia 
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entropy change, a reversible process must be devised so that the second 
law of thermodynamics applies, i.e. 

ds = faz (193) 
The total work done by each gas specie as it is expanded isothermally fron 


the initial total preasure to its equilibrium partial pressure is 


we -HIZn; in nj (194) 
per mole of mixture. By the first law, 
AE = aw (195) 
but AE = 0 for all isothermal changes of ideal gases. Therefore, 
q@ = -FE Xn; ln nj (196) 
or AS = -Rin; In a; — (197) 


Although the effect of pressure on ideal enthalnoy is zero, thie contribution 
must also be sdded to the total entropy. For cne mole of mixture, 

AS = Rin pp (198) 
Canbining these three contributions, the total entropy of the combustion 


products in the canbustion chamber ia 
Ss, = in, Ss, - ngs Ella) Bp R ln pp (199) 


where c chember entropy 


Ss. = 
ny = moles of the j th product 
8; = ideal entropy of the j th product at T, 
np = total moles of gaseous products 

py = total chamber pressure 

R = universal gas constent 
When solids or liquids are present as products of combustion, they are 
assumed to occupy no volume. Consequently, they are included in the ns; 
summation but do not contribute to the entropy of mixing or the effect of 
pressure on entropy. As was presented earlier, the entropy of mixing of 


liquid products is considered negligible. 
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Te Exhaust Composition Iteration 
The products of combustion are expanded isentropically through the 


converging=diverging nozzle to an established pressure and an unknown 





temperature, Again, using mass balances and equilibria relationships, the 


exhaust composition is calculated at the exheust pressure and an assumed 
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exhaust tempsrature, As before, iteration is continued until the quantity 
of each product converges within a pre-determined degree of accuracy. 
8. Entropy Balance 


After determining the exhaust composition at the asaumed exhaust. 








temperature. an entropy balance is made. The belance is satisfied when the 


— 


change in entropy between the chamber and HozZl]e sxit is ZETO, 166s 
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In all probability, several aseumptions of temperature will heave to be made 
before a solution is obtained. Both the exhaust composition iteration and 
the eutropy balance are repeated until the solution is reached.. 

9- Exhaust Enthejoy 


After determining the exhaust temperature and composition at constant 


entropy, the enthalpy of the exhaust products is evaluated an 


Tt. 
H, =|{n x, + In, A (201) 
8 | 1"tdogg 3 F298 


Previously, the enthalpy of the chamber products was evaluated as 


r 
H, = [20 *slase + Zny Pt (202) 
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It is important to note that the By are the actual equilibrium concen+rations 
of the various product species at the exhaust and chember temperature, 
respectively. Thus; the totel enthalpy change through the nozzle is 


AH = H, ~ 5, (293) 


10. Rocket Performence Parameters 
From the previously calculated properties, specific impulses is 


rigorously evaluated aa 


(a) Ten-| [2x Gig =H) (204) 
& 


The volumetric specific impulse, the apecific heat capacity, and the mean 


molecular weicht of the products may be rigorously evaluated as 


(b) Ivap =f: Isp (205) 

(c) co, = Enjoy, (206) 

(a) MW = 100 prems proveljunt (207) 
Rp 


it should bs acted that 100 sg of propellant is a atendard basis upon which 
to calculate propellant performance, 

The preferred method of evaluating the other parformance parameters 
may depend upon the nature of the propellant system, When all products 


of combustion are gases, the following approach is rigorous; 


Me? 2 =zxkes i 
Sy T- in t/t, (208) 
ln Pg/Ps 
(fe). alt fi-y (209) 
°F k2.7 k+1 Te 
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eer aca NE ner nce = ERE 
Tyan = A Pa RE emma - 


| mat 
a ETE) (2 (210) 
A™ Cy kt] t Pe 
(bh) c* = Igpee (211) 
Cy 


However, these equations require some adjustment when multiphase flow existe. 
i.e. when solids or liquids are present as products of combustion. In this 


case, the following approach is preferred: 





(z) Cy = Cp (212) 
Cy Cpeige 
ce Fin 


= 2(100 gu propeiie is (2th) 
0 Aa assy aaa) = 2+567 
1/2 ar) 


(4) g =al/2r y ; ES (215) 
c* | a [ kt] 

(3) +) = Tepe 

3) cy (approx depeg (216) 

(ik) dy (approxe) = 6 Dy Pate (217) 
a* c* Py a 


Since different values of specific impulse, exhaust temperature, and specific 





heat ratio will exist for shifting and frozen flow, the parameters DoAg/m, cy: 
and A,/'A* will also vary. The characteristic exhaust velocity c* can be 
rigorously evaluated only if throat conditions are accurately known. However, 
this requires sane detailed point-to-point calculations. A good approximation 
to c* is obtained by using the chamber temperature and specific heat ratio in 
the above squation, 
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B. Method of Minimization of Free Energy 


1. General Concepts 
For the determination of the equilibrium composition of a complex 


mixture, the only thermodynemic property required is the Gibbs' standard 
free energy. Earlier, this wes shown to be 

a¥ = Vdp - SaT (218) 
under adiabatic, reversible conditions in which al] work is accomplished by 


expension, This is probably the most useful combined atatement of the 


the independent variables. However, a more direct approach to the 


calculation of propellant performance may be achieved through the basic 


first and second laws of thermodynamics since pressure and tempereture are 


definition of free energy 
F=H-TS (219) 


Substituting the expressions of enthalpy and entropy developed in the previous 


section, 
Fe inj, - Tins, * Rlnop 2 in ny + RTngln pp £220) 
np 
where = total free energy of the combustion products: 


F 
nj = moles of the i th product 
H j = enthalpy of the j th product 
S5 = entropy of the j th product 
T°= temperature 

R = universal ges constant 

Dp = total gaseous moles 

Pp = total pressure 


The equilibrium composition may be determined by finding the eet of positive 
values nj which minimize the previous equation while simultaneously 


satisfying the mess balances, 





A procedure for applying the method of the minimization of free 
energy to propellant performance calculations was first presented by White, 


Johnson, and Dentzig.43 In their approach, the free energy function F/RT 


ode OF 


is minimized at constant temperature instead of utilizing equilibrium 


ie 


equations as wes done in the previously discussed calculation method. 
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Another feature of their approach was the derivation of the method of 


"“gteepest descent" for solving a symmetric set of linear equations by the 


matria vounaique in the iteration procedure. Only m #1] linear equations 


in m * 1 unknowns need to be solved where m represents the number of 


elements préessnt in the mixture. Thus, the number of simultaneous equations 
40 be solved is usually mch amelier than the number of products of 


combustion. 
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This advantage is perticularly sdaptable tc linear progremming 
for a lerge, high-opeed digital computer. However, the principal limitation 
of the White, Johnson, and Dantzig method is that, as developed, it is 
applicable only to complex mixtures in which all the constituents are gases. 
This is unsatisfactory for the evaluation of many prototype propellant 
systems in which solid or liquid products of combustion are formed. 
More recently, an improved technique for the minimizetion of free energy 

vas developed by A. B. Roberta. Instead of minimizing the free energy of 

@ complex mixture at constant temperature, Robexte’ method allows both 
temperature and canposition to vary while maintaining either constant 

enthalpy or constant entropy. This mettiod has been progremmed for tne 
Remington-Rand 11034 canputer and ie ‘applicable tc evaluating the equilibrium 
composition of complex mixtures conteining solid, liquid, and gaseous 
constituents. Since this method was used for conducting the performance 


calculations presented later in thie repart, the important progremming 
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concepts as well as the canputational technique will be discussed in detail, 


2. Establishing the Mess Balances 
In programming notation, Roberts’ method desigi.ates the elemental 


mass balance ant 


z ve M, = Ey (221) 


where i = ingredient element 
j = exhaust specie 
E, = total gram-atoms of tiie i th element present 
in the propellant formation 
. = moles of the j th exhaust specie present in the 
equilibrium mixturs 
ay <o atoms of the i th element present in the j th 
exhaust specie. A; ; must be a positive, whole 
integer. 7 


M 


The maximum possible concentration of any given specie is defined by the 
minimum quentity of any ingredient element that is available from any of 
the applicable mass balances, 1.€. 
0 <M; S$ min By = Bj (222) 
aij 
Thus, the most restrictive mass balance defines B;, the maximum number of 
moles of a given specie j which can possibly be formed. For example, the 


value of Bj for CO, would be smaller of the following: 





0 Srco, 2 = By (223) 
For programming convenience, let x5 = Mie Then, the mass balance may be 
Tre-stated as Bj 

= Aq ay xj = By (225) 

j 
or 

24,5, B, = 

po 2 4 (226) 





E, 16 


Now, let 


Ss es aa ee (227) 
Ayj 


Both xj and 8, j very in value between O and l,. 
Thus, the mass balance may be expressed in standard normalized form as 





for which the coefficients &ij are stored in matrix form. The quantities x4 
f are proportional to the mole fraction of the respective exhaust species, 


The total moles of gas are expressed as 


| mn 
Mo= Z My (229) 
t jel 
where m= gaseous elements 
n = gaseous molecules 
Pp = condensed phases 


The maximum possible number of gaseous moles would cccur if ths entire 


propellant formulation existed as gaseous elements. Thus, the limits for 
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| 
| Mo are 
| | O<My BE, = By (230) 
| Letting My = By Xo. 
| = (231) 
Bo Xp = By xj 231 
=) 
| & 
or 
mtn By 
X= & = x (232) 


between 0 and 1. It should also be noted that Bj xj Tepresents the moles of 
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The quentity x, is ano.ker progremming convenience which varies in value 
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any species j present in the equilibrium mixture while By Xp represents the 
total gaseous moles at any one time, From Dalton's Law, the partial pressure 
contributed by each gaseous specie j is 

xj By 
To Bo 
where By, By» snd pp remain constent while x, and x, ere variables. 





Py = Py (233) 


The various x;'a indicate the degree to which a specie contributes to 


completing each masa balance. To insure operability, the x;"s are ordered 


in magnitude for each trial solution, 1.6. all species j which contain element 


i are listed in order of decreasing concentration. Since the x5*s are 
proportional +o the mole fraction of each respective specie j, this erdering 
procedure establishes priority for choor’... a basis upon which to define 

the denendent variables. The species possessing the greatest x5 are chosen 
as the basic variables, with these species being selected ss the besie in 
the mass balance for which they possess the greatest A, 4° For exemple, any 
concentration of COp contributes more to the oxygen balance than to the 
carbon balance. ‘Thus, C85 would be comsidered as a besis only in the oxygen 
balance. Then, it could be eliminated fran the cerbon balance. In this 
manner, a set of majors and minors is inherently established, 1.6. a limited 
number of dependent basis variables are chosen while the remaining epecies 
are independent variables. The number of basis variables ig limited to the 
number of elements in the propellant formulation. The concentration of any 
of the independent variables may be changed by an arbitrary amount. Then, 


the concentrations of the basis variables are determined from the mass 


balances and the concentrations of the independent variables, When the order 


of priority changes, the basis is revised. Both the mess balance equations 
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and the total gcescous moles equation ore maintained throughout the trans- 


formations of the matrices, Previously, the free energy equation was written 


tk ¥ = Znjily ~ T2n8, + Hg yo moe Ring ln pp (234) 


Incorporating the program notation, this becames 





intntp mentp mtn ) B i j . 
@ F= Sf H4(Byx,)-T LY S,(Byxs)*RT Y = (Byx,) In B +R B,yx; 1n pp 
: fog te fat TS Foy 5 % Fai J 
(235) 
. or 
; _ mint aa 27 
. Fez B jx (Hy - 1S,;) + RP s Byx,(1n Byx,) + RT Box, (1n pp) (236) 
ie j=l jeirvs Box 
e o~o 
i where F = molar free energy 
7 H = molar enthalpy 
S = molar entropy 
er T = temperature 
a; = Olea of exhaust specie j 
ft B; = maximun possible moles of specie J 
1 £ x", = ratio of actual moles of specie j to the 
La J maximum number of moles of apecie § 
. L R = universal ges constant 
7 Np = total gaseous moles 
By = maximum possible total moles of ges 
Xx, = ratio of actual total moles of gas to the 
| 9 = ratio of actual totel moles of gas to th 
: maximum number of moles of gas 
|; pr = total pressure 
| i m= gaseous elements 
E n = gaseous molecules 
7 Pp = condensed pheuses 
if 
if 36 C) a ational A 
Fram the basic relationship 
a¥ = Vdp - SdT (237) 
it follow: that, for an ideal gas, 
2E| =-s (238) 
oT > 
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For a mixture of ideal gases, the additional restriction of constent 
concentration must be applied to the previous equation, This may also be 


shown from the definition of free energy 


FH -TS (239) 
Then, ak = dH - TaS ~- Sat (240) 
ar aT af af 
ak = 6, - 18 -S (241) 
dt av 
However, as = Cp (z42) 
ar? 
Therefore, ag=-S <0 (243) 
at 
a°F = -aS = -C_ <0 (244) 
2 
atari 


These equeti ons express the effect of a change in temperature on the ideal 

free energy of a complex mixture at conatent pressure and concentration. 

The first derivative dF measures the slope of the free energy curve while the 

aT 

Second derivative é*y messurea the curvature. Both the slope and the change 
a 

in slope ave always negative with the second derivative being numerically 

smaller than the first. Therefore, the curvaturs is always concave downward. 

The slope of the free energy curve is a good approximation to the actual 

curves over a fairly wide renge of temperature. Furthermore, this approximation 

is a conservative one in that the tangent to the curve always lies above the 


actual curve, aa illustrated below: 














F(T + AT) R(T) -S(T) AT (245) 
Are-s(t) At (246) 

Thus, free energy may be readily estimated over small temperature intervals 
by the first derivative. It is important to note that these ecuations 
apply rigorously only to the ideal free energy at one atmosphere pressure. 
Neither the effect of mixing nor the effect of pressure is included. 
Whereas chenges in the partial molal free energy arise from changes in 
temperature, changes in the entropy of mixing and pressure arise fron 


changes in concentration. 


Both of these latter quantities may be incorporated into the total free 


energy summation through «he use of a supplemental equation of the form 
y=xing (247) 


where y = free energy of mixing and pressure 
x = concentration of the various gaseous species 


Yor the normalized form of the equation where OS@xie1, the absclute minimum 
y occurs at x = 1/e and has the value y = -l/e, where e is the base of 


natural logarithms. 
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The slope at the y-axis is infinite while the curve interaects the x-axis 
at O and 1. Since the calculation of the x ln x function is a relatively 
lengthy procedure on a digital computer, Roberts has developed a quick 
approximation which permits a reasonable decrease in free energy for a 
modest calculati onal effort, 4° This approximation is made by Soiateueting 
a@ parabola which is tangent to the y = x ln x curve at the initial point x 


and which passes through the origin, as shown below: 


y = x(ax-b) 





This is also a conservative approximation since the change in free energy is 
never es great on the parabola as cn the x In x curve except at the point of 
tangency. One such approximaticn is made for esch specie at each trial 


composition. Furthermore, 100% convergence is mathematically guerenteed. 





For the true free energy curve 


y=xine (248) 
‘ah ay fy kee (249) 
while for the parabolic approximation 
y = x(ex - b) (250) 
ae dy/ax = 2ax - b (251) 
caabining these into two simultaneous equations. 
= in x-= ax” = bx (252) 
x(1n x +1) = 2ax® - bx (253) 
the constante are found to be 
a=} (254) 
x 
bel-iInx (255) 


Veing this approach, the approximate change in y = x ln x corresponding to 


[ =e el + ln | (256) 


(y +Ay) = (x +Ax) In (x +Ax) E(x +Ax)| x 


a change in x ia 


(y+ Ay) ¢ (x +hx) (444 1n x) (257) 
Ay = A(x in s)Ax[1+1nx+ 3] (258) 


It is amportant to note that x, the concentration vuriable, is actually the 





quantity By, used in the program notation. Also, during the progress of a 
calculation, both T and x are interim numbere. The A values are the emounts 
of change employed in the next trial solution. Changes in coticentration are 


always made for one spesis et a time, choosing that specie which will refleot 
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the greatest change in the total free energy of the system. For each change 
in concentration, there is a corresponding change in temperature, and vice 
versa, By conducting a number of smll, incremental changes, the free-energy 


of the system is rapidly minimized with « minimum of calculetional effort. 


4. Details of the Calculation Method 


The caleulational routine is initiated by assuming that the entire 
propellant composition exists as monatomic gases. This, of course, is a 
hypothetical situation which defines the maximum free energy and the maximum 
total meles of gaseous products for the system. Also, at constant enthalpy, 
the corresponding temperature of the monatomic gaa mixture is a very low 
hypothetical one. In turn, the various products of combustion are introduced 
incrementally into the complex mixture. Each addition serves to reduce the 
free energy and increase the temperature of the mixture until equilibrium 


is attained. As shown previously, the free energy of a mixture of canbustion 


products is 
F Fae (H, -TS,) + Ping aay) Mie 
= n - + ln + a mo 
jr Jodi J = nT Dy (259) 


For a single gaseous specie of the mixture, the free energy is 


mtn 


ree ny (Hy - 18, + Rfny In nj) * RP ng In 7 (260) 
where F = free energy 
H = enthalpy 
S = entropy 


T = temperature 
n, = moles of specie j 
R°= universal gas constant 
= total pressure 
= total moles of gaseous products 
m = gaseous elements 
n = gaseous nolecules 
p = condensed phases 
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For any incremental change in the concentration of the combustion products, 
there will be three contributions .o the chenge in free energy of the 
mixture: (1) the free energy contributed by the formetion of the new 
product; (2) the free energy lost by the decrease in concentration of the 
reactant(s); and (3) the free energy chenge introduced by the variation in 
the total number of gaseous moles. The method of handling these changes in 
the calculational approach may be best shown by an exemple. Consider the 


formation of nitrogen gas from ita atans: 


2N = N,, 
Any = ~2Any, (261) 
An, = An, * Any, x. ~ Any, (262) 


4 change in the number of moles of N, creates an inherent change in the 
number of moles of Ne. At a given temperature and concentration, the free 


energy of the mixture may be expressed as 


~ ¢ 
Fe ty | ~ TSy cn] * RT In ny + mnt, t-255, 009] Fin, 18 Dy 
*eoe + RT ed in br (263) 
oe 


Since constant enthalpy must be maintained, any change in concentration will 
necessitate a change in temperature of the mixture. Over a short temperature 


range, Gc, Av? AH. 


Therefore, 
A= 0= Amy, | Hy, (7) - 2447) + Zajc, Ae (264) 
Hy, ~ 2 Hy 
Ate a A (265) 
LM | TM : 
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Thus, a coefficient is obteired which relates the change in temperature of 
tha mixture to the change in concentration of an individual specie. The 
dneremental free energy change created by e decrease in the concentration of 


N is 
tT +AtT r+ At 


AT = ny ki = my) - 2 Any, [Bs = 15, (266) 


+ R(T+ AT) [on - 2Any,,) In (ay - 2 Any.) ~ Ty ln my| +R - AY (ny In ay) 


This includes the ideal free energy and the free energy of mizing. Similarly, 


the incremental change in free energy caused by the formation of N, is 


12 +AT tT +AT 


2 -T : 
AF = my, [i Sy, | tat Eo 7,,| (267) 
+ R(T+AT) (on + Any, ) In (Dy, + Any, - By, ln t,|* Re AT(ny, in ny.) 


Finally, the incrementel change in free energy due to the variation in total 


moles of products is 





P, P 
Ar = x(t + At) (np - Any.) 1m — = - (268) 
2 


However, this may be more conveniently re-stated as 





Av = -R(t +2) (Amy) In pq + BOAT np an PE + R(T *A2)| op In a, - 


ap (269) 
(ng = Any) in (np - dm) | 


Now, by applying the AF approximetions previously developed, a number of 

coefficients may be Obtained which are functions of AT andAn. However, 
? 

since AT =k An, the approximate change in free enezgy may be reduced to 


& cubic function of An, 








First, consider the ideal contribution to free energy. From the 
curve of free energy versus temperature at constant concentration, 
Ar = #1 + AT) - F(T) S S(T) AT (270) 
Therefore, 


o o 


Constitutes a contribution to the total change in free snergy. This 


aot. contribution is the product of a constant and AT. Also, 
: F(t + Av) SFr (r) -s(t) At (2723 
Therefore, 
oAn. (BOY + A (rr \ 
~< “ON, SAN‘ p Ny, \ No/p 
and 


a a 
~24n,, (3), Art Any, (Sy) AT 
are also contributions to AF. 


Second, the entropy of mixing will be considered. From the parabolic 


approximation to the function y = x ln x, 
A(xinx) < ax[1 +10 x+ Os | (272) 
Therefore, 
2An 
AF SRT +AT) (-2Any,)| 1 + In ny - = N21 + R(my in ny) AT (273) 


+R (t+ AT) (Any) + 1n By, + Ady, + R(ny, in ny.) AT 
2 


Fran this expression, the following additional contributions to the change 


in free energy are obtained: 


(1) R(ny + In ny) AT + R(my, ln Py, ) At 
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(2) RT(-2 Any, ) (1 # In ny) + RT (Any) (1 + In ny) 


No 


(3) RE(-2Any) (-2An, ) + RE (Any) ( Any,) 
Py "Np 


(4) Re AT\-2 Any) (1 + ln ny) + ReAT(Any,) (1 + In ty) 


we 
™y ny, 


(5) RAM -2Any,) (-2Any,) + RAT (Am) (Amy) 


Finally, consider the effect of a change in the total moles of gas on 
the free energy. This contribution contains a term of the general form 
-A(x In x) where x = np and Ax = - Any,» Beceuse of the minus sign, the 
seme approximation used in the entropy of mixing contribution will not 
apply here. The desired quadratic must satisfy the requirement 

-A(x in x) = A(-x In x) quadratic in Ax (27h) 


eas shown by the illustration 
y 


dae 


; x 


This inverted quadratic approximation becomes 


arenr+ An)|-Cn by 1) ( Brg) Ena 2) * eg 1 haa (Am9?| 
(Mypax > Bp)? (275) 


By applying L' Hoapitels’ principle, the coefficient of ( Ang)* basomes 
1/2 when np = ape Thus, the following contributions to free energy 
erise from the effect cf tots] moles: 


(1) -RT ( Ang) In pp 
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(2) -Re AT ( Ang) 1n py 
(3) *Re AT ng In pefny 
(4) -RE (in op +1) Ang 


(5) -Re AT (12 ny + 1) Any 


ie 
(6) srt agg, ~ Ba? * tage tna ( Ang)? 
(ae x - nf) 
; 3 bor § 
Ota eae "Tina ” "Tex 


(7) +n°-At — 2 
(=, > np) 


Now, the total incremental change in free energy, es a function of temperature 


(Ang)? 


and concentration may be represented by 


" Cm * y C] 2 r 2 P 
AF (t,x) § sls AT + Bi Ax, *Gy Ax, AT * DC Ax.) + EB Az,) as] (276) 
However, utilizing the approximatioa 
T=k'Azx (277) 
the above equation reducesz to one with omly three coefficients and one unknown: 


' ae eee : 
Ar<é ae Ax, + Bi Ax,) « Cy Az] (278) 


The value of Ax which minimizes AF is found by differentiating this equation 
and setting the derivative equal to zero: 


Ax = By + a4)? - 3 ayes (279) 
3 o§ 
Thus, a conservative approximation has been developed which can be solved 


directly to give Ax. Furthermore, the true change ik free energy exceeds 
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that found by thie approximate solution. In order to maintain the mess 
balances, the concentration of any specie cannot become negative. At 
the seme time, the concentration cennot exceed the maximum defined by 
the wess balances. These two restrictions serve er limits for the 
walues of x found from the solution of the cubic. 4s equilibrium in 
approached, the eapproxamate solution comes closer to and finally 
coineides with the rigorous evaluation of sanposition, temperature, and 
fres energy. in initiating the paretolic approximation, it is assumed 
that z = 2/3 xyey instead of x = 0 so that the constants a end b 
will be finite in value. The position x = 2/3 x,,, ie the point of 
tangenoy for which, om the average, the parabolic curves best spproximetce 
the y = x ln x curve. The important concept is to utilize a parabola 
whiech always lies above the 7 = 2 lu x or the y = -x ln = curves. 

‘The summation of the various contributions to an incremental 
change in free energy yiclds, fa the 2N—" Np exemple, the following 
values for the constants of the cubic equation: 


AY Teta( 5p * ERR ny Im ay) = 207i * (FR Dp = ORE (2 + Lo my) 


# RE (A © In ty) - RE (ln np #2) ~ RE Am py # KIB oy in (280) 
Bt =< — - 2k"( Bye * k!( 3 )e - 2k'K1 # In my) + k’R (282) 
(2 © dn my) RPA Lo mp +2) © RAR Im pp « Ret (MMe © Be) * Pyar 18 get 
(Saaece > He)? 
or ARR. EER. «oy POE” MT) * Maes 12 
a i) 


een: | (262) 
(nex = ty )* 











5e Introduction of Thermodynomic Data 
Thermodynamic data for the various products of canbustion ere 

conveniently introduced im equation form as a function of temperature over 
@ specified temperature interval. Of course, reference values of enthalpy 
and entropy are also introduced at a reference temperature within each 
interval. In order to insure an accurate fit of the thermodynamic deta 
over ths complete range from 298°K to 6000°K, several equations may be 
utilized, if necessary. However, in general two equations are sufficient 
to provide a highly accurate fit. Over the low temperature range fram 298°K 
to about 1200°K to 1500°K, specific heat capacity data are expressed as an 
ascending polynaniel of the form 

eS a+ bP + cr + at « ert (283) 
Fron thia, the enthalpy, entropy, and free energy of the individual species 


are evaluated as 


Hs fojat= Ht + or + yi? + on) + att + op5 (284) 
2 3 4 5 
s'= [Cg =s' # ain T+ bf + CT aT? + oft (285) 
‘i is sae ae 


F's H-TS'= Ht = TS’ + eM - In T) - br? of = att ~ of) (286) 
2 12 20 


where H = enthalpy at any general temperature T 
8 = entropy at any general temperature T 
F = free energy at any general temperature T 
C_ = specific heat capacity at any general temperature T 


H' = reference enthalpy 

8' = reference entr 

T = temperature x 107. 

@, Db, Gy dy © = empirical constants for the curve fit 
Over the high temperature range extending fran about 1200°K to 1500°K, the 
specific heat capscity date are expressed es e descending polynomial of 


the form 
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c seer + or”? « at73 + et (287) 
P 


Then, the enthalpy, entropy, and free energy equations over thia interval 


become 
H CaT=Ht+aT+binT-G-ad-¢e (288) 
| P T ore 383 
g=fcat=s':salnT-b-o -a-e (289) 
[oe tT or 3P 3r 
F=H-1S =H! -1S' + (aT-b)(i-ln Tr -g-g.- (290) 
2r 6r 1203 


Although any reference temperature within the corresponding temperature 
interval may be utilized, 298°K and 3000°K, respectively, were chosen a” 
convenient references for most species, The former ig particularly desirable 
since it is a standard reference for many experimental data. By using the 
two reference points, dusl integration is not necessary over both ranges 

of the heat capacity curve. 

The functional method of introducing data into the calculation hes 
both a thermodynamic and a progremming advantage. Since the equation 
constants and the reference data are known, temperature and canposition 
are the only variables involved in the calculation of enthalpy, entropy, 
or free energy. For a trial ccmposition, temperature may be computed 
directly at either constant enthalpy or entropy. The progremming advantage 
arises fram the great saving in data storage and from the rapid solution of 
a polynomial equation on a high-speed digital computer. The actual storage 
of tables of thermcdynemic data at 100°K temperature intervals would be 
quite cumbersome, Furthermore, the process of table look-up and interpolation 
is slower thar the calculation of point properties from the polynamial 


equation. 
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6. Evaluation of Thermodynamic Properties and Performance Parameters 





After evaluating the chamber temperature and camposition by ths 
minimization of free energy, the entropy of the system is calculated as 
was done in the method of equilibrium constants... Then, the exhaust 
temperature and composition ere svaluated by minimizing the free energy 
of the system at constant entropy, and the exhaust enthalpy is calculated, 
From these data, the various performance parameters are calculated in 


the same manner as for the method of equilibrium constants, 


@} 
Wwe 














CHAPTER V 
SUE THEORETICAL ANALYSIS PROGRAM 


A Selection of Prototype Propellant Systems for Evaluation 


The object of the theoretical calculation work presented in this report 
is to evaluate the performance of several prototype rocket fuels and oxidizers. 
In general, good fuels and oxidizers should possess relatively endothermic 
heats-of-formation and yield combustion products with large exothermic 
heats-of-formation. Furthermore, only the most exothermic chemical reactions 
are worthy of consideration for application to rocket propulsion. Since the 
most exothermic reactions occur betweer, the elements in the upper-left and 
the upper-right corners of the periodic table, thie study was confined to 
fuela end oxidizers consisting of the elementa carbon, hydrogen, oxygen, ni- 
trogen, chlorine, fluorine, lithium, beryllium, boron, aluminum, and magnesiun. 


The following prototype materials were chosen for evaluation: 


Fueis 
Name Forma Name Formula 

Liquid Hydrogen Ho Liquid Oxygen Qo 
Hydrazine NH), Liquid Fluorine F 
Unsymmetriceal Dimethyl Oxygen Difluoride OF, 

Hydrazine CoN5Hg 100% Hydrogen Peroxide H,0, 
Polyethylene CH, Nitrogen Tetroxide N, 0), 
Dicyanoacetylene C),No Nitrogen Trifluoride nf 
Pentaborane pote Chlorine Trifluoride cif, 
Aluninum Hydride Perchloryl Fluoride C10 
Beryllium Hydride Bel? Nitronium Perchlorate re 10), 
Lithium Hydride Li Ammonium Perchlorate 4,010), 
Magnesius Hydride MeH, 


4 





en ni, 





These materials were chosen in order to obtain the maximum amount 
of basic information for a modest calculational effort. By evaluating 
Several fuela with a given oxidizer, a relative order-of-merit may be 
established for the fuels, Similariy, the evaluation of several oxidizers 
with each fuel reveals the relative merits of the oxidizers, Each of 
the yreviously mentioned fuels was evaluated with each of the oxidizers. 
While only oxygen and fluorine were considered worthy of general study 
as oxidizers, five light metals, in addition to carbon and hydrogen, 
were considered as prime fel ingredionts. The clemente nitrogen and 
chlorine were essentially considered as carrier elements. 

The study was established systematically to reveal the effects of 
various fundamental chemical moieties on overall performance, For 
exemple, the interactions of the fuel groups N-H, C-H, C-N-H, Al-H, 
BeH, Li-H, Ba-H, and Mg-H are evaluated with the oxidizer groupe N-F, 
O-F, Cl-F, Cl-6-F, N-0, Cl-0 and H-0. 411 caloulations were conducted 
in accordanse with the assumptions and procedures discussed eerlier, 
Yor each fuel-oxidizer combination, several mixture ratios were 
eveluated in order to define the optimum performance ani provide sane 
indication of the influence of propellant composition on performance. 
Although many of the fuel-oxidizer casbinations in this report are 
not workable propellant combinations, this approsch satisfies the 
ovjective of obteining basic theoretical combustion information. 
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B. Thermodynamic Properties of Reactants and Products of Combustion 


1. General Thermodynamic Requirements 


When performance calculaticns are conducted rigorously as described 
previously in the method of equilibrium constants and the method of 
minimization of free energy, the accuracy of the results is governed by the 
accuracy of the basic thermodynamic data utilized. For reactants, only the 
heat-of-formation at a reference temperature must be known. However, this 
quantity must be known accurately because of its importance in the enthalpy 
balance. If either volumetric specific impulse or ideal boost velocity is to 
be calculated, the density of the reactants must also be known. For products 
of combustion, additional thermodynamic date must be available, Both the 
entropy and enthalpy must be known as a function of temperature in addition 
to the standard heat-of-formation., With these data, a consistent set of 
thermodynamic functions can be camputed. Usually, the conventional 
thermodynamic functions of C,, S, ~(F-Ho9g)/T, H- Hogg» and H-H59/T are 
tabulated at temperature intervals of 100°K, Sometimes the heat-of-formation, 
free-energy-of-formation, and equilibrium-constant-of-formation are also 
included in the seme tabulation. 

The availability of accurate thermodynamic date for application +o 
rocket combustion studies deserves some special camment. As shown by the 
performance calculations, combustion temperatures of 3000°K-/,\000°K are common. 
These obviously exceed the temperature range of experimental studies and 
require the use of thermodynamic data which have been computed from statistical 


mechenics, When the molecular structure and the basic moleculer translational, 
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rotationa), vibrational, and electronic contributions are known, the statis- 
ticel mechanical approach is quite accurate. However, a fundamental protlem 
in high temperature thermodynemics ia molecular stability. New apecies 
which are insignificant or unknown at room temperatures may be quite stable 
in high temperature complex mixtures. An exataple of this phenomonen is 
metaboric acid, HBQ., which replaces B20, as the dominant boren specie in 
the B-0-H system above 2000°K, Frequently, insufficient data are available 
on the high temperature decomposition and chemical squilibris of species 
which are quite stable at lower temperatures. This shortcoming can only 


be rescived by sxperimentel research. 


2, Heat-of-Foymation of Reactents 
The standard heats-of-formation for the reactants used in thiz 

theoretical study are, for the most part, quite accurately known, These 

are listed along with the densities of the reactants in Table I, Only the 

: heats-of~formation of 41H, and BeH, are reasonably in doubt. Since these two 

hydrides have not been prepared in pure quartitios, sscurate experimental 
heat-of-formetion data are lacking. Estimates of the AH» cf AlH, and Bel, 

l vary from between #20 Keal/nole and ~-15 Koeal/mole. Whereas the correspondiag 

hydrides of lithium and magnesium are readily synthesized from the elements, 

attempts to prepare AlH, and BeH, in this manner have been unsuccessful. 

These data, along with the 4l-H and Be-H bond strengths, indicate that the 


qrecee ea Se a 


AH, of both hydrides should be endothermic, Therefore, a Ale of 
#10 Koel /aole was chosen for Alls, and a AHp of +10 Koal/mole wes also 
selected for BeH,. 


ee hep ee ee 
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TABLE I 


HEATS-OF-FRMATION AND DENSITIES OF REACTANTS 


Reactant Heat-of -formation Density Ref. Temp, State 
Keal/mole g/ec OK 

Ho ~] .887 0.071 20.4 Liquid 
Nou), 412,05 1.004 298.16 Liquid 
CoN5He 412.,7h 0.786 298.16 Liquid 

CH 6.50 0.920 298.16 Solid 
c ie +118.0 0.970 298.16 Liquid 
BoHy +7 7h 0.620 298.16 Liquid 

Al +10,0* 1.72* 298.16 Solid 

Be +10,0* Unknown 298.16 Solid 

; Lin’ 21.6 6,820 2928.16 Solid 

MgH, -17.9 LAS 298.16 Solid 
: Fo “3.17 1.52 85.2 Liguid 
; OF +326 1.53 128.0 Liquid 
e: N20), 5h 1.43 298.16 Liquid 
& NF, 32 02 1.532 144.0 Liguid 
iy C1l¥3 ~hy ol 1.810 298.16 Liquid 
4 C10,F -10,1 1.43 298.16 Liguia 

NO3C10 +8.0 2.22 298.16 Solid 

NH),C10,, 69.42 1.96 298,15 Solid 


i 
; 
i 
J 
& 
¢ 


* Estimated Values 
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3. Tbernodynemic Properties of Combustion Products 
The current status of the thermodynamic properties of combustion 


products varies from excellent, as in the case of hydrogen gas, to purely 
theoretical, as in the case vf such metallic oxyhalides as AlLOCl. In 
most cases, the properties of the dominant species are reasonably well 
defined for engineering applications. This is particulerly true for 
products of combuation formed from the elements carbon, hydrogen, oxygen, 
nitrogen, chlorine, and fluorine. These data have been reported by Huff, 
et. al.o47 Gordon4®, stuli and Sinke49, Rossini, et. al.,°° and others. 
In general, somewhat less is known about the thermodynamic properties of 
combustion producta containing the métais lithium, boron, beryllium, 
Magnesium, and aluminum. However, two relatively recent reports by Evans, 
Prosen, and Wagman provids a thorough summary of the properties of boron 
compounds te 52 Still another recent report by Beckett, et. ales 
provides a summary of the known thermodynemic properties of lithium, 
beryilium, eluminum, and magnesium and their compounds with oxygen, 
hydrogen, fluorine, and chlorine.” Additional compilations of therm- 
dynamic functions neve been made by Gordon 54+ 55, Stull end sinke5®, 
Henderson, cet. ale!’ 58 ond others which are applicable to the vere 


formance calculations presented in this report as well as to calculations 





involving other elements. All thermodynamic data used in this report 
are swmmarized in Table II, III, and IV at the end of this section. In 


Table II, stenderd values of heat-of-formtion and entropy are as mmarized 





at reference temperatures. Constants for the heat capacity oquation 


Cy zat wT + cf? + ATF + etd (291) 


are presented in Table IIY end constants for the heat capacity equation 


Sn Lene ae Set Wey Oe 


Cpe & + oe? + oP? + are « et-4 (292) 


are presented in Table IV. These curve fits were provided by Messra. C.B. 


Pane oan 
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Henderson and R, S. Scheffea? Unless otherwise specified by subscript, 
all species listed in Tables II, -II, end IV are gasewme molecules, 

In somewhat more detail, the current status of the thermodynamic data 
on the metal-containing species may be summarized for various temperature 
rengese Fran 0°K to I000°K, experimental data are rather complete and fairly 
accurate for the species stable over this temperature range. Heats-of- 
formation and heat capacities have been measured for practically all of the 
aignificant solids and liquids. However, the chemical instability of the 
hydrides end the hydroxides, with the exception of Li0H, has precluded 
accurate measurement of their vapor pressures end heats of vaporization.©° 
Low volatility of same other compounds has hindered vapor-phase studies in 
this temperature range. 

Experimental difficulties increase in the range of 1000%-2000°K, but 
these are not sufficiently formidable te prevent most thermodynamic 
measurementa. Specific heat capacity data are camplete for same refractory 
materials, such as MgO, but usually dic not extend to the melting point far 
substances wiich melt between 1000°K and 2000°K, such as 1150.5 The avail- 
ability of vapor-condensed phase equilibrium data for the stable species helps 
to compensate for lack of other data. The sewer techniquesz of optical and 
mass apectrametry and molecular beams have greatly assisted in identifying 
the composition of high temperature gas mixtures 52 However, mere accurate 
afalyses of such mixtures is needed. Experimental heat-of-fusion deta ere 
needed on most campounds, 

Above 2000°%K there are very little data on condensed phases, Specific 


heat capacities have been extrapolated from lewer temperatures. Thermodynamic 
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properties of the gas species are quite accurately known for the mona- 
tomic and many of the diatomic species. However, in the polyatomic species 
and somes of the diatomic species, moleculer constant deta are needed. Also, 
dissociation energies need refinement in most cases .o2 

A number of polyatomic gas species have been observed at low pres~ 
sures by vaporization of solids and liquids in mass spectrometers. In a 
few cases, the electron-diffraction and infra-red spectra are aiso avail- 
able. However, the configuration and molecular con¢g.vants have been estab~ 
lished with reliability for only a few of these species. 

Basic low temperature data are available on the more significant alum- 
inum compounds, such as 41,096 ALCL ys and ALF 46 At high temperatures, the 
existance of ALP, and AlCl, has been proposed on the basis of indirect evi- 
dence. Also, the uxyhalides AlOF and A10C1 and the interhalides AlF,Cl, 
AIFC1,, and ALFCL may be stable molecules, although no direct evidence for. 
their existence is available. These molecules were first proposed by 
Gordon. 4 Their fundamental vibrations and bond distances have been es~- 
timated from values for zelated molecules 65 Equilibrium studies of the 
Al-Q system indicate that 41,0 is the predominant vapor species and that 
42,03 does not exist as a gas. The heat-of=formation and stability of 
crystalline aluminum nitride are probably the most questioneble of the 
light metal nitrides. Experimental neats-of-formation vary between -57 
and ~77 Keal/mole, with ~-70 +10 Kcal/moie being the best evailable es- 
timate for AlN(, )° However, #ince AlN is a minor constituent in a fluor- 
ine or an oxygen atmosphere, the uncertainty in its properties is relative 


ly insignificant in thie study. 
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In gencral, data for the beryllium compounds are not as well known. 
Both low and high temperature heat capacity data are needed on Be(0H).., BeFo, 
BeCl,, BeF, and BeCl, but molecular constants are available. The mixed 
halide BeFCl hes bien preposed but not ver fied experimentally. Data on 
beryllium oxide are somewhat more complete, including the high temperature 
vapor phase. Gaseous polymers containing two to six BeO molecules have 
been identified by Chupka in the mass wpectrocrapic’? Although this is of 
basic interest, the effect on rocket performance is not significant since 
the condensed phase is the predominant oxide specie. Solid BeO is more 
volatile in the presence ef water vapor, suggesting the formation of Be(0H)o. 
The properties of BeOH ars estimated values. 

The basic low and high temperature properties of the more significant 
magnesium compounds are fairly well defined. Some refinements would be 
useful in the molecular constants of the halides, and direct evidence hes 
not been found for the existence of the mixed halides. Since magnesium 
oxide is an extremely high-melting compound, the gaseous oxide species are 
relatively unimportant in this study. 

Experimental low temperature deta are available on all important stable 
lithium camnpounds except LiCl. Some data are also available on these species 
in the high temperature range although the number of investigators has not 
been plentiful. Probably the most questionable data are thoge near the 
transition temperatures of LiF, LiG1, and Lin0. The polymeric gaseous 
halides (LiF)5, (LiF), (LiF), s(LiGl)o, and (LiCl), have also been reported 
by Chupka.°7 The thermedynemic properties of these molecules need to be 


further refined, particularly since the lithium helides are geses at normal 
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rocket combustion tenperatures. 

The products of combustion of boron have been atudied more extensive- 
ly than those of either lithium, beryllium, alwninum or magnesium. How- 
ever, the problem is a more complex one, particulariy in the case of the 
boron oxides, B203¢ the principal stable oxide, melts and vaporizes at 
much lower temperatures than 41203, 3s0, or MgO. This means that the 
vapor species piay a more daminsnt role in the combustion of 2 borone 
oxygen system. Also, the possible namber of high temperature boron 
species excseda that of lithium because of the greater valence and 
ccordination number of boron. Frobably the most controversial questions 
concerning the thermodynamics properties of B20, arise fran the 
incvasistensy of the molecular constant data with the thermosiemicel 
data for that molesule; Also, the structure of the crystalline state 
has not been absolutely resolved. The most probabis sources of error 
fin the thermochemical data are the heat capacity of the liquid and the 
heat~of fusion. 8 The gaseous species 1,30, and HBO: are significant 
at low and high temperatures, respectively. In fact, HBO, is the 
predominant boron apecios in the 5-0-H system at temperatures above 
2000°K, Although the heat-cf-fcrmation of HBOp has been questionable for 
sums time, recent data have fairly well established a maximus AHe of 
140.8 Koel /mole.°9970 Bven more recent studies by White indicate the 
Ally of HBO2 to be -234.5 Keal/mole at 298°K, The data for the halides 
B61, BCl3s BF, end BF are considervd reliable. The existence of the 
gubhelide BCl2 hes been suggested by the existence of the Cimer B2Ci,. 
Data for BF, heve slao been estimated similarly to those far BCla, 

AlFo, ond AlC1,- Direct evidence is available for the existence 
of BPo0l, B¥Clg, and BFCl, although the data on these 
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molecules need refinement, (“9 The oxyhalides BOF and BOC] are hypo- 
thetical molecules whose properties have bean estimated assuming a 
linear molecule. 

In general, accurate heat-of~formation data are needed most for con- 
ducting theoretical performance calculations. However, another area of 
almost equal importance is that of high temperature properties. Both 
qualitative and quantitative deta are needed. Experimental measurs= 
ments cf high temperature heat capacities and heats~of=transition are 
needed for solids and liquids. Although these studies may require new 
experimental techniques and equipment, they will undcubtedly refine the 
lata now available. More detailed discussions of the thermcdynamic prone 
arties of the reactants and products considered in this study may be 
found in the original literature citations listed in the footnotes and 
the bibliography. 

In Table II, the enthalpy data at 298° are the standard heats-of- 
formation of the various species, i.e. the heats-of~reaction for the 
formation of the species from their elements -xisting in their standard 
‘states, At higher temperaturea, the enthalpy consists of the standerd 
heat~of-formation at 298°K plus the heat content above that temperature. 
Of course, the entropy data are absolute values at e pressure of 1 atmos- 
phere. In the two heat capacity equations, T is expressed as degrees 


Kelvin times 1072, 
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STANDARD VALUES OF ENTHALPY AND ENTROPY FOR PRODUCTS OF COMBUSTION 


Reaction 
Product 


Bs 


H,0 
HEL 
bet 


Cs) 


o(e) 


CO. 
0 


of 


No 
NO 
F 
Fo 
cl 
Clo 
CF 
OF 
or 
CF), 
CH 
CH, 
CH 
2 
CHF 
CHOFo 
NH 


NHQ 


aw? 


Gok 
HC. 
COFs 
Al A 
ni) 
Al 
A1C1 
AlCl, 
alcL 
ar? 
AlFo 


AHp 
Keal /mole 


+52 .089 
0 


-57-798 
“22,063 
6h 2 

Q 
#17089 
26,416 
94,052 
-? #55 


#933 
#113,00 
Q 
+21.60 
+#18.903 
G 


#2892 
0 


#830 
=18 0 
#120,0 
-218,0 
+143,0 
+660 
“17.89 
+5422 
162 a 6 
-105,.5 
+77.0 
+h, .0 
-11.04 
+9420 
+738 
+31.20 
-150.5 


#775 
+6104 
-11.5 
=76.0 


s 
cal/mole/°K 


27.393 
31.211 
15.106 
4h 0645 
412526 
1.361 
37.762 
W7 214 
51.070 
38.469 
49.011 
13.888 
36.61 
452763 
50.339 
37.917 
48.477 
39.457 
536291 
50,889 
58.25h 
62.352 
62.635 
k2 34 
Kha 
4.59 
48,00 
62.04 
58.92 
43-30 
"7.96 
46.034 
48.41 
57-91 
48 eel 
61 82 


39-31 
He Bh 
5443 
63.025 
75073 
51-450 
57 +330 


TABLE IT 


Temp, 
°K 


298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 
298 


298 
298 
298 
298 
298 
298 
298 


105 


iat 
Keal/mole 


465.511 
+21 208 
27 598 
=O e 198 
~h3 * 121 
+14.78 
+184 .44 
= 4.059 
57.520 
+73,071 
423 455 
+30,786 
4#126.437 
422.160 
+. 32h 
+33.641 
424.192 
+l) .066 
42) 416 
4106.45 
417.405 
-70 142 
#153243 
+16)..9 
*97.0 


#101.0 
-101 5 
-h7 056 
+98.77 
$71.19 


4117.2) 
4124.95 
+66.79 
-101 el 
421.330 
+90.88 
84.79 
+12 .90 
=36.671 
-87. 61 
~35 419 
oD 3 874, 


Ss 
cal /mole/°K 


38,862 
48.470 
68.419 
62 364, 
58.709 
12.30 
49.29 
65.370 
79 84 
50.095 
67.980 
61.376 
48.089 
63-756 
68.854, 
496791 
68.436 
51.623 
73-761 
70.122 
86.974 
102 .163 
113.842 
60,03 
68,88 


84.08 
109,16 
102 26 
60.93 
71.4348 


66.98 
98.77 
75 86 
106.56 
25.217 
50.842 
63.842 
75 «120 
96.914 
120,086 
71.552 
90.768 


Temp. 
°K 


3000 
3000 
3000 
3006 
3000 


3000 
3000 
3000 
3000 
3000 
3000 
3000 


3000 
3000 
3000 
3000 
3000 


3000 
3000 
3009 
3000 
3000 
3000 
3000 
3000 
3000 
3000 
3000 
30c0 











Reaction 
Product 


ALF 
AlFCl> 
ALFoo1 
ALFCL 
A10OC1 
BLOF 
ALG 
A120 
41205(s) 
Ab205(2) 
LG 
rr ae 
<g) 
‘anes 
(s) 


e 


TABIEZ II CONTINUED 


Ale s Temp» H 
Keal/mole cai/mele/°K oK Keal/mole 
-283.0 66-11 298 #231012 
185.0 72.430 298 -132 2523 
-228,0 70.885 298 -176.619 
-121,0 61.962 298 ~82 .035 
56.0 60.580 298 “14,818 
-100.0 57317 298 61.677 
#16.8 52.16 298 +40.61 
039 oh 61 2200 298 42160 
-,00.2 lz e 172 298 =339 3 
~313.63 
105.0 63.765 296 50.56 
+104..0 51-707 298 8=— #12826 
35567 15.89 298 -346.19 
~70.0 5.20 298 63 .896 
+9880 
: 422,630 
4135 02 36.649 298 4148.62 
4204.0 48.232 298 4227 «607 
4117 41 O47 298 4139.63 
~60.3 3.57 298 2 6239 
4155-0 50.712 298 #178 169 
446.25 50.944 298 #70 2233 
=-2).8 592134, 298 “14.807 
9761 69.320 298 45 6342 
2 .87 47 6892 298 -19-583 
-154,0 53-235 298 116.99 
-270,0 60.71 298 -220.31 
89.0 572376 298 -51.130 
-210,.0 65 688 298 -159-597 
-154.0 68.784 298 -102 .859 
-62.0 56.706 298 24.164 
-118.0 5h 3h2 298 =80 505 
+16 .0 48.609 298 #38015 
-105 14 57233 298 53 058 
-1h0.8 572660 298 ~94..269 
281.65 49.38 1006 
=300.98 18.55 298 =29h4..36 
cy) 2-28 298 
410.7135 
*77 +9 32 055 298 #912338 
+76.0 42.19 298 498.744 
413.0 51-993 298 +37 2036 
84,0 57.87 298 hh e762 
“16.5 49.151 298 +6 .963 
-192.10 52.256 298 -154.65 


S 
cal/mole/°K 


109.17 
116.643 
114.976 
95 e341 
9.216 
90.855 
71.897 
89.710 
69.736 
86.977 
108,54 
71-776 
34.6907 
152738 
11.15 
16.47 
48.121 
67.850 
59-314 
15.10 
69.616 
70,881 
91.10) 
111.940 
66.935 
82 2547 
100.32 
87.936 
106.312 
110.488 
87.50 
84.674 
67.026 
98 Ah 
93769 
131.11 
111.70 


33-15 
208,82 


13.684 
402 

60.782 
711.920 
89.828 
68.399 
82 294, 


Temp. 
OK 


3000 
3000 
3000 
3000 
3000 
3000 
3000 
3006 
2313 
2313 
3000 
3000 


906 


AN: 
oun 


3000 
3000 
3000 
1200 
3000 
3000 
3000 
3000 
3000 
3090 
3000 
3000 
3000 
3000 
3000 
3006 
3000 
3000 
3000 
3000 
3000 


600 
3000 


1556 
3000 
3000 
3000 
3000 


3000 
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TABLE II CONTINUED 








Reaction AHe s Temp. H s Temp. 
Product Keal/mole cal/mole/°K xX Keal/mole cal/mole/PK °K 
BeF (4) 3222 93 34.468 1076 
BeFa(s) -241.7 9.000 298 
BeFC “134.0 55.960 298 =95 «yO 87 630 3000 
Be0(5 “14301 34.370 298 =111.243 28.754 2823 
Boots} 90.213 340776 2823 
Be0¢ a) +30.1 47.208 298 453.282 66.140 3000 
(Bev75 ~67 .026 85.570 22h2 
(Be0)3 -197 2539 115.771 2242 
(Be0)), 2301.052 146.668 22h2 
(BeO) 404.565 173-868 22h2 
(BeO): 514.078 196.907 2242 
Be(OH)o = =2.559.72 54.964 298 ~134.463 87.365 1500 
BeOH -28,0 54.202 298 +3..787 T.G27 3000 
Be( OH a “216.1 13. z 258 
NV #121. 67 800 
rig. sige aSc140 298 #524965 hk 5 3000 
Li($) +968), 19.78 1600 
Lis #50 .47 47.021 298 #75 421 672998 3000 
Lil +3220 40.78 298 455 253 59.97 3000 
Licl 49.0 50.834 298 ~24,.623 71.191 3000 
(LiCl)o 88,50); 111,891 3000 
Lil(,y © -80.3h 38.074 1200 
LiF -77460 47.754 298 53233 67.67h 3000 
(LiF) -176.21 96.463 3200 
(LIF), -260.66 155.421 3000 
Lio 13.8 46.65 298 +*36..99 65.60 3000 
Ling) 242-4 9.06 298 
Ls50¢5 ~68.68 63.61 3000 
L420(,) — - 35.80 55.820 298 “1.19 83.430 3000 
Lidn © 3.0 54.370 298 25.95 83.871 3000 
LioH =53.8 53.366 298 ~22 238 77.923 300¢ 
(Li0H)> -87.76 118.619 3000 
LiF(s) 1.46 23 8.57 298 
LiF(y) «115.7 3954 2000 
LigN(s) 47.5 8.90 298 
Mg( 1) +7.01 18.29 1000 
5) +35 ..6 35.51 298 #49 .038 46.97 3006 
+41 46.11 298 +64.573 65 443 3000 
MgCl +2,0 55234 298 #28.646 76 hah 3000 
MgCla, g) -103.7 61,65 298 64.83 94.324 3000 
MeC1.5( 1) 129.76 5459 1000 
MeF 41.0 52.814 298 -16,362 72 948 3000 
MeFo ( gj -169.7 5592 298 -130.66 88.344, 3000 
MeF (1) ~215.91 58.32 2000 
MgFC ~137.0 60.054 298 -97 822 93.607 3000 
MgQ(p) 13-9 50.9k4 298 437.991 70.992 3000 
MgO) -89.402 40 234 3073 
£23 -107.90 34.214 3073 
MeCl2(s) ~142.06 42,00 900 
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Reaction 
Product 


MEeFo(s 
2 ) 


Mg( 0H 2 3 


Me3N2(s) 
c1lF 
¢10 


AH, 
Keal/mole 


“24410 
~56.0 
~221,.0 
“110.3 
~13.423 
+4,3..00 


TABLE II CONTINUED 


s 
cal/mole/°K 


ue ,82 
56.981 
15.09 
21,0 
52,062 
540379 


Temp. 
°K 


1500 
298 
298 
298 
298 
298 


108 


3 
Keal; mole 


"95-99 


-0 e 8510 
453.501 


s 
cal/mole/°K 


83.126 


48,098 
65.808 
68.217 


Temp. 
°K 


3000 
823 


1500 
1500 








TABLE III 


SUMMARY OF CONSTANTS OF Ba HEAT CAPACITY EQUATION 
6. = a9 + bT + CTS + aT3 + ef+ FOR PRODUCTS OF COMBUSTION 








Pp 

Reaction Temp. Range 
Product °K a b a a e 
H 298-1200 46491 2.1104 =449036 4.7721 #1.6529 
Hy 298-1200 4.9782 12.6667 -28.6774 = 27.5445 = =9.2560 

Q 298-1200 7 eT 4G “0.4000 4.6918 “2.7618 0.5388 
Hi 298-1200 57282 = 78535 -18.3714 18.8757 +-6.6140 
HF 298-1200 5.3188 21.0072 -25.9848 25.7227 ~-8,8189 
Crs) 298-1200 2008755 1121421 9 =3.1995 =hehi8, 244909 
o(6) 298-1200 5e53K5 305997 Selah6 7.7652 2 644o 
Cc 298-1200 6.5135 2.367) =5 4030 8.1190 326364 
CO, 298-1206 4.3090 20.6041 =20.9840 11.8419 2.7832 
0 298-1200 55202 -1.1625 0.2062 1.0710 0.6299 
Op 298-1200 8.1634 ~10.067) 28.6774 -27.5958 9.1338 
on 298-1200 7.2173 069874 -3.7090 52660 -2.20h7 
N 298-1200 4.9556 0.0600 -0,.0889 0.0410 i) 
N 298-1200 722990 =2.2393 0750 ~0.6623  =-0,6229 
N6 298-1200 8.0151 6.7019 15.9410 12.5923 3.4646 
F 298-1200 6.8092 ~8.0466 16,0915 -14.8332 5.0394, 
r 298-1200 520585 12.5570 18.0551 12.6196 = =3..4646 
Cc 298-1206 1.5072 2362522 9 4902232 he 3739 1 SG 
Clo 2982200 5.1523 17.0695 31.3058 26.3340 9-8. 2644 
cr 298-1200 6.7024 ~0.3563 8.2324 -9.1886 3.0680 
cr, 298-1200 721238 16.4083  +2.9395 -3-6248 = 1.9745 
CF 2981200 1 07148 50.8439 ~61.7657 36 4601 8.5419 
cry 298~1200 0.6148 66.6450 =77.7637 =e 3853 9.9271 
CH 298~1200 72683 “1.8710 3.1512 ~1.0671 0,02430 
Cc 298-1200 8.0892 2.6220 9.9827 6.3896 1,318) 
Co 298-1200 2.2834 = 8264 5927199307790) =12.6.2479 


CHF 298-1200 0637K7 5904746 6926683 4322455 = ~10.9599 
CHD¥ 5 298-1200 327552 21.6690 5.1725 “17-4742 6.8712 
NH 298-1200 70346 0.2360 9 -0..6454, 203496 1.0510 
298-1200 8.1021 2.5106 93232 ~6.4932 1.6039 
298-1200 701234 2.3491 9.7882 -7.5632 1.7072 
298-1200 706062 -541455 12.8420 “10,0281 2.7218 
ON, 298-1200 75571 2966794 -=39420639 28.2055 ~8.1166 
298-1200 4e691h = 18.4680) 22.2652 15.0556 9 +. 3682 
COF2 298~1200 206503 = - 3967223. 23295 = 22 57K6 =~ 8420 
al, 2981200 309677 79856 “20.0784 20.3835 7.2441 
} 900-6000 7G 0 0 ) ) 
rN 298-1200 7.7836 = =7.6791 0. 21.6722 919.2239 57959 
ACL 298-1200 6.4388 9.9343 715-7794 11.5560 ~3.1835 
AICly  —- 298-1200 8.2309 25.9667 41.3807 30.1295 8.2564 
acl, 298-1200 10,5301 35-6621 56.1667 40.5592 11.0510 
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Reaction Temp. Range 
Product °K 


AIF 298-1200 
4lFo 298-1200 
AlF, 298-1200 
AlF¢l, 298-1200 
AlF,C1 298-1200 
AlFC1 298-1200 
A1061 298-1200 


ALOF 298-1200 
AlQ 298-1200 
A190 298-1200 


Al50 298-1200 
4rz030$) 2100-6000 
A102 298-1200 
AIN(G) 298-1200 
AIN( 5) 298-900 

AIPS(g) 298-727 

ALFS(s) 727-2000 


B13 2100-6000 
Brg) 298-1200 

298-1200 
wfic8) 298-1200 


BN(g) 298-1200 
BN;,) 298-1200 

{ 298-1200 
Bl, 298-1200 
BCL, 298-1200 


BF 298-1200 
BF, 298-1200 
BF 298-1200 


i 298-2200 
BFZCE 298-1200 


298-1200 
BOC1 298-1260 
BOF 298-1200 
EO 298-1200 


Oy 298-1200 


298-1200 
ap, pe 
0 298-1200 
52098) 298-606 
( 35 298-1200 
Bers) 298.1700 
Bec) 14300-6000 
Be(,) 2982500 
fe 298-1250 


TABLE III CONTINUED 


501754 
4.9599 
57138 
6.7952 
6.8416 
7.1218 
8.3706 
5 +9524 
4.9828 
524563 
6.5714 
31.15 
7.7923 
58327 
547 
2.580 


‘20.93 


75 

364917 
528345 
TehO53 
4.5298 
701326 
5 02489 
4.9226 
6.5126 
5e314h 
6.1558 
326440 


10.4212 
4.7597 
0.1952 
5 04370 
52 030 
~] 2880 
413 
75 
497 
10.6340 


b 


12.1638 
3.0210 
4.9180 
46.2482 
4304470 
26.3110 
1964358 
25.5204 
11.9091 
22.9840 
13423794 
2.0 


41.8087 
11.4582 
7.80 
86.775 
1.50 


0 
9.8630 
527108 
=3.2771 
27 8125 
-35555 
11.3558 
32 5208 
4321581 
8.3476 
73366 
39 00299 
52 e624 
~8.0392 
40.4030 
22.5827 
10,6706 
706457 
19.6261 
24.2014 
791232 
391950 
488.83 
161.8829 
24 

0 

0 
22.5709 


110 


¢ 


“16.5747 
4966949 
-61.7929 
-69.0741 
62.1709 
-38.1612 
=25, 4114 
332.2289 
“16 6 4703 
2767457 
=203.2832 
0 


57.9975 
#175701 
3] 

157.717 


0 
-23.20913 
#2 ohh51 
721383 
26.4321 
14.0143 
“lh. 82 33 
46.0491 
~60.2 744 
11.7021 
12.9662 
«2.9112 
40.5830 
47.1216 
5221268 
=29.6029 
-3-9988 
19.9985 
-19.6396 
19 e 4200 
“91.5570 
3404995 
-1190.0 
189.244 
0 

0 


¢] 
47.4823 


a 


10,8250 
3400750 
4064308 
48,0722 
41.4930 
26.0522 
16,0854 
19.7106 
11.3880 
16.2804 


14.0387 
0 


38.2374 
12 Ahh 
0 


14525 
0 


0 

22 6147 
1.7599 
=4.0601 
14.4284 
~13.0869 
9.2787 
30.9468 
39-7437 
929308 
-20.8438 
2307216 


“2.2705 
11.6047 
8.2347 
55925 
15.0008 
1016.67 
111.294 
0 


Q 
Oo 
-38 07665 








-2 7339 
~8.9490 
-10.2369 
-J2 07522 
-10.8748 
6 8226 
33-9854 
4.7327 
~351496 
~3-7667 
38.7424 
i) 


-9. 7448 
- 3414 


“4805 
0 


0 

-7 o8740 
1.2598 
0.6299 
306332 
41677 
2.25140 
8.0194 
-10.1276 
-3 0466 
7.8372 
3a 
17.4848 
-729951 
5.3220 
1.5553 
503543 
=-3.0012 
“1.4398 
13 e 7120 
Bg 05699 


25 69599 
Q 
0 


Q 
11.2880 


t 
re 
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TABLE III CONTINUED | 


Reaction Temp. Range 





Product °K a b c d e 
Beol 298-1200 52314 1144345 14.9863 904369 ~2.3086 
BeCl> 298-1200 51303 31,6832 =i 6868 29.9397 7-700 
BeF 298-1200 6.1152 323155 0.9714 7.0498 1.1968 
BeFo(g) 298-1200 2 4903 31.1361 -3469557 = lve 2050 9 =. 1.859 
BeFo(7) 1000-6000 21,0 Q 0 fo) 0 
BeFapq\ 298-1200 8.4 12,0 0 Q 0 
BeFO2X"’ 298-1200 506566 26,5295 735-2020 22.7686 5-716 
Bes 298-1200 =526027 61.5775 +93.2205 58,0057 18.8518 
Beats} 2700-6090 1h.5 ) Q 0 i) 
Be0, 298-1200 70077 “2.8128 12.7682 12.5039 349368 
Be( fi}. 298-1200 -0.6491 79-4228 -119.3368 86.4588 #23 7485 
Be OH €98-1200 74132 0.5027 141321 1406854 467205 
Bes) 1200-3000 11.93 0.91 0 ) @ 
Be3N2(s) 298-800 732 30,8 0 0 o 
4( 2) 298-1200 he97 0 6 0 0 
Li(y) 500-6000 6.87 Q © co) 8) 
Lig 298-1200 703496 6.8209 -10.6540 9-79.35 21993 
Li 298-1200 6.6048 -0,7687 10.0314 “10.8285 3.5662 
Licl 298-1200 5.6195 11.9363 --18,.0048 12,6189 -323657 
LiF 298-1200 5.2328 10,6854 12.9571 70580), “1.7254 
Lio 298-1200 6.9059 ~1.8393 10.5388 -10.6287 323779 
Li2Q(g) 298-1900 2.4011 51.2102 62.7011 38.2906 8.5721 
1490(}) 1500-6000 24.0 @ 0 ty) Ny) 
1440.) 298-1200 6.1995 961561 1.7129 726438 31166 
11,8" 298-1200 2.6280 «2942237 3305974 =. 045 =l06172 
140 298-1200 8.3729 6.3251 2322359 -2025663 6.1118 
LiF(y) 1906-6000 15.5 0 0 G a 
1ACl(} 880-1200 1720691 9-152 ) i) o 
Lt} 298-800 11.73 23.0 0 0 0 
Mg) g00-6000 8.1 0 0 Q ) 
Me ) 298-2500 497 0 9 Q a 
nei 298-1200 6.9 -3.2669 14.7073 “1404920 4.6052 
Mgol 298-1200 6.5499 905397 915.1524 1241001 -3-0559 
MgClo(g) 296-1200 8.1795 26.3921 202976 = 30.9229 <3.4994, 
42(1) 900-6000 22.1 0 o 0 0 
298-1200 54028 12,1908 9-17-5182 = 12.8730 =3 0802 
Mel 2( ¢) 298-1260 41208 38.3454 5ST 5026 = HO 12 -10.6744 
MeFo(7) 1400-6000 2226 to) co) 9 ) 
maret 298-1200 5.8100 3327002 52.0245 37.0813 -10,0122 
MgO 2) 298-1200 5.2592 12.0747  -16.6978 11.0012 2.7946 
MgO 2900-6000 14.5 0 © (3) 0 
mgoee 298-1200 5 8757 13.4434 13-5188 7o22ly 1.067 
Mes) 300-3200 9.965 1.805 Q 0 ) 
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Reaction 
Product 


MeFo(s) 
"2(S 

Me3No(s) 

clF 

G10 


Temp. Renge 
°K 


400-2000 
1061-2000 


298-1266 
298-1200 


TABLE III CONTINUED 


16.47 
20.77 
20.07 
28.5 

5 05350 
=2,0832 


b 


2276 
11.2¢ 
10.66 
Q 


10.5214 
620494 


112 


aaca 


13 e 6836 
367874 


m& 


QaagQgo0 


8.2632 
=2 3555 


® 


oeaeos 


-] 28898 
0.5893. 











TABLE IV 


SUMMARY OF CONSTANTS OF THE HEAT CAPACITY SQUATION 
C= a + ptt + cT™* + aT"3 + ef4 FOR PRODUCTS OF COMBUSTION 


P 

Reaction Temp. Range 

Product °K a b c ad e 

H 1200-6000 4.9680 0 4) G 0 

Hy 1200-6000 13216.5 <27.4278 67.2718 ~64.8607 40,7688 

sft 1200-6000 13.7762 1.8182 -19,5586 21.5310 7.7672 
1 1200-60600 10.2217 9 =7.2079 «=:15..688 22,0974 11.6068 


1200-6000 12.3357 =25.0605 67,8036 ~91.7611 45.9920 
og 8) 1200-6000 8.9014 “14.7665 31.8293 =3329478 13.2068 


o¢e) 1200-6000 507187 0.3176 “13.1090 27.1718 -16.0 
co: 3200-6000 8.9548 3.1902 ~17.8004 28.4934  -16.0 
CO. 1200-6000 17.0619 = 10.1588 9 ~-17.6949 -18.8447 72457 
0 1200-6000 6.7126 #1325630 39-3302 49.8513 23.2136 
G, 1200=6000 11.0832 ==3.3193 -S..4192 22.2363 9 ~12,.6224 
off 2200-6000 11.1146 12,8166 27.9832 36.0147 17.8456 
N 1200-6000 11.2577 52.8970 162.330 =214.4616 102.7204 
No 1200-6000 924233 1.7978 0.1459 0.8427 1,0156 
NO 1200-6000 10.1998 -8.4006 23.0720 31.863 15.8143 
F 1200-6000 49810 0.3996 «1.3628 2 03392 1.3058 
F. 1206-6000 8.8708 0.8054 =30422h 50525 <2 .9017 
c 1200-6000 5ei251 -1.7062 9-7-2114 “10.6191 5.6584 
Cl 1200-6000 11.0049 = 12.3026 = 33.9264 “43.7925 21.0372 
cr 1200-6000 9.7813 ~3.6223  7+1708 “8.0166 3.2911 
CF. 1200-6000 13.9133 -0,03111 -1,0264 -0.1186 00,2110 
cr 1200-6000 23.5138 -20.5833 52.6576 66.2369 30.0422 
oF, 1290-6000 30.2094 = =21.8359 70 4 96 “52.9712 21.6034 
CH 1200-6060 10.5378 8.2322 16.2034 -20,1049 9.2827 
C 1200-6000 13.8598 0.5278 -9.6784 7.0023 ~1.3080 
G 1200-6000 21.0765 <2.7760 2.0961 “17.1159 13.3333 


Cc 1200-6000 25.8030 0.3565 -8.018), 4.2406 ~0 6329 
1200-6000 10.5927 8.44398 15.7380 -19.057i 8.8608 
CN 1200-6000 9.4034 25-5485 129.9354 177.2701 86.6245 
COND 1200-6000 21.3002 -4.0233 8.6897 -15.9300 8.4383 
HCN 1200~6000 18.0749  -15.9472 31.9731 39.9297 18.4810 
Aig) 1200-6000 6.7710 =16.8320 56.1522 -79.2737 40.0435 
1200-6000 11.1719 9 --11.7417 27.5952 -32.2290 13.8397 
A1C2 1200-6000 10.4069 =7.7297 = 18.7755 ~21.3683 9.0717 
AlCl 1200-6000 14.9053 0.01525 ~0.1658 -0,07478 0.04219 
A113 1200-6900 19,8693 0.01952 -0.3761 0.04007 0 
Ar 1200-6000 5 elk 12.1638 16.5717 10.8250 +-2.7339 
Al¥, 1200-6000 14.9017 00,0238 -0.6166 0.1412 =0,042 
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Reaction Temp. Range 


Product 


AlF 
AFCl> 
al Fct 
ALFC1 
aiocl 
ALOF 
A10 
Alod 
4103(8) 


0 
'3( 8) 
moe) 1) 
BoOy g1) 
pel 
BeC1 
BeClo 
BeF 


erei® 


Belg) 


°K 


1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-2500 
1200-6000 
1260-6000 
900-3000 

1000-2400 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
14200-6000 
1200~6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-£000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
600-6000 

6N0-6000 

1200-6000 
2500-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 
1200-6000 


TABLE IV CONPINUED 


a 


19.8763 
19.8743 
19.8738 
14.9100 
14.9088 
14.9030 
8.6384 
13.9092 
55 20332 
20,8688 


8.9423 


14.699 
10.3983 
563886 


18.0317 
9.9021 
10.2413 
14.9077 
19.8782 
8.7610 
14.9054 
19.7823 
14.9062 
1928719 
19.8794 
14.9077 
14.9050 
10.2272 
20,867) 
407597 
3726781 
25.8298 
30.5 
30.5 
67.4342 
21.9304 
10.9483 
10.2809 
14.9066 
10.5587 
14-9061 
14.9098 
9.9108 





b 


-0,0514 
-0.0300 
-0.0247 
-0,06719 
605508 
0,0087 
303156 
6.011336 
84.0657 
~0,03878 
9,001028 
0 
~7.0962 
4.5609 
~1.4373 
2.4833 
«37 -909k, 
5.0205 
6.8557 
-0,0398 
=0.0711 
1.9896 
=-0,0149 
0.9185 
-0,022h 
~0 0004 
-0 0826 
-0.0334 
0.0029 
9.1175 
0.0048 
24.2014 
1.21163 
0.0601 

(0) 


6 

1.7731 
“148. W244 
-10.5468 
“Te 0561 
=0,0294 
-10 25803 
-0.0226 
-0.0631 
5 ) 0763 
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c 


0.7156 
~0 04986 
-0 6365 
0.1884 
063275 
07949 
-12.9668 
105.1027 
-0,6914 
0.1281 
-1.789 
4.206) 
16.9404 
3.0706 
95 +3357 
1163352 
16.4790 
-0.5047 
#05586 
8.2963 
-) .6857 
5 03915 
1.1104 
-1,6333 
=0.8952 
-1.4548 
=] ,8188 
25-9717 
~3-1005 
19.4200 
=20 0253 
1.6868 
0 

0 
-28.7114 
435 «4418 
24.1330 
16.9462 
~0.5268 
28.9426 
“1.3434 
-0.6087 


11.5275 


a e 
—0,2812 0.2116 
“0.1797 0.1266 
=0,1499 021266 
-0.3579 0.2110 
-0,3052 0.2110 
0.0626 9.0422 
19.9806 11.0265 
0.07978 0.04219 
“KZ. 432 0 
«0,206, 0.1668 
0.008015 0 
6 Q 
15.8978 923332 
15.2569 ~7-8346 
=£11.5386 47.6000 
#1325906 5.8650 
-18.9494 8.0591 
0.1671 0.1266 
-0,.3812 0.2532 
12.1788 -6.6739 
0.0081 0.2532 
501252 0 -2 4G 
-0,0553 0.1688 
0.0604 0.1688 
0.4174 0.3376 
0.2612 0.2532 
-36.2764 18.1357 
0.5032 0.4641 
8.2347 21-4398 
14.5691 -3-2912 
0 0 
0 ) 
22.4K54 62447 
~428..4375 5200 
28,6829 12.4894 
“19.4619 8.2700 
0.1476 0.1266 
~36.5533 16.6245 
-0,0556 0.2110 
-0.3033 0.2532 
“1307952 569494 
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TABLE IV CONTINUED 


Reaction Temp. Range 


Preduct °K a b c a e 
(Be0)o 1000-600 19.981 0 -2.207 0 a 
(Be0), 1000-6000 31.957 0 -3.278 0 0 
(BeO) 1000-6000 43.870 ) -3.083 ) 0 
(BeO)« 1000-6000 555«846 © hel 54 0 a) 
(Be0)z 1600-6000 67.760 0 -3959 0 7) 
Be( 0H)o 1200-6000 26.7677 0.8694 -14,.3231 12.1060 -3-2489 
Bed ~ i2d0-6000 23.8695 0.6035 -9.5250 8.8324 267004 
Be(GH)2(s) 298-6000 35.211 0 -1.903 ) i) 

0 -3.067 0 0 


Be No(s) 800-2600 36.753 
Liz 1290-6000 13.64by 24 7845 60.3401 -68.3563 28-9873 
iif £000 1248739 “B5eS9T 364999? 6958 18.1856 
Licl 1200-6000 8.9450 0.0290 -0.0397 “0.1542 0.0844 
(LiG1)> 1000-6000 19.969 0 -1.969 0 | 
LiF(,) 1200-6000 10,8294  -9+9162 23.8040 -27.1839 11.5190 
(Li 2 1000-6000 19.969 0 -1.969 0 0 

(LAF )3 1000-6000 31.788 0 -1.788 5 a 

Lio ~ 1200-6006 8.9422 0,0040 -0 «7480 @.0737 0-084). 
Li20( a) 1200-6000 13.9102 0.00539 =1.7148 0.1367 0.2110 
Li0 


1200-0000 14.9121 «= -0 0823-15936 0.2421 0.4642 
Lio# 1200-6006 1328676 0.6212 -9 9390 8.8958 ~2,6160 
(LA0H)p 1000-6000 31.735 0 <2 £735 0 0 
LiCl’, ) 1200-6000 15-2 C0) rs) 


0 6 
Mg ,,) 2590-6000 25 8927 -188,3282 566.1888 -566,0000 0 
3200-6000 13-8517 #1503516 36.6664 ~j2..3675 16.1013. 


MgC1 1200-6009 10.5230 8.3132 20.1329 ~22 48263 926624 
MgClo( 2) 3200-6000 14.9119 -0,0885 0.1132 <0.4B42 0.2532 
1200-6000 10.2757 70273 16.9432 ~19.3802 38.2276 


Meo 2) 1700-6000 14.9101 0.0675 -0+2527 0.3625 0,211 
mares 1200-6000 14.9049 000114 -0-3259 00553 0.0422 
Mgo(,) 1200-6000 10,3480 724032 17-8102 “20,3834 8.6498 
1200-6000 13-8733 0.56hl © -B 8562 8.6828 -2. 7426 
MaClo¢g) 298-1200 19.917 © 0.258 0 0 
(cis 3)298=6008 39.392 23-946 24293 0 e 
(8)3200-6000 20.7734 97-8628 295.5712 385.2596 162.2271 
10 120026000 -135-59H2 -750.9014 1852.353 -2111.363, 8991563 
14g) 1200-6000 1566816 -82.3779 2321-0344 287.9056 132.0276 
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CHAPTER VI 
RESULTS 


A. Sunmary of the Results 


The performance of 100 prototype propellant systems has been evaluated, 
with each of the ten selected fuels being considured in combination with each 
of the ten oxidizers. All calculatione were made at a combustion pressure 
of 1000 pei expanding to exhaust pressures of 14.696 psi (1 atmosphere), 2.0 
pai, and 0.2 psi. The latter exhaust conditions were chosen to simlate 
high altitude operation where large nozzle area area ratios are desirable, 
The temperature, enthalpy, entropy. heat capacity, and molecular composition 
of the combustion products were determined for the chamber and each exhaust 
condition. Then, assuming shifting chemical equilibrium accompanied by 
complete phase squilibrium, the following performance parameters were 
evaluated: 

) Specific impulse. 
) Specific heat capacity ratio. 
) Volumetric specific impulss. 
) Thrust coefficient. 
) Nozzle area retio. 
Also, the following parameters were evaluated for the condition of complete 
frozen equilibrium: 
(1) Specific impulse. 
(2) Characteristic exhaust velocity. 
(3) Thrust coefficient. 
(4) Nozzle area ratio. 

Since the prime purpose of this study was to evaiuste the chemical 
performance of prototype propellant systems, the shifting specific impulse 
has been considered the most important single parameter. As shown eerlier, 


this represents the maximum energy attainable from a chemical propulsion 
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system. This, and some other significant parameters, have been plotted 
in a series of three grapha for each propellant system. The first pre~ 
sents both the shifting and frozen specific impulse as a function of the 
propellant composition for each of the thee exhaust pressures. The 
second presents the corresponding exhaust semperatures at each exhaust 
pressure. Finally, the third graph presents the chamber temperature, 
the chamber molecular weight, and the characteristic exhaust velocity 
es a function of the propellant cumposition. These graphs constitute 
Appendix C. Unfortunately, the volume of output information from the 
computer was much too large to be incorporated in its entirety into this 
report. Two pages cf output were obtained for each chemical formulation 
considered. However, since the complete output contains additional data 
which was not plotted, inc]uding the molecular composition, the completes 
output is presented for one fuel-to-oxidizer ratio of each propellant 
syatem. This constitutes Appendix D. The molecular composition is of 
particular importance since it indicates the relative concentrations of 
the various combustion products as well as those which were considered 
in making the calculations. it should be pointed out that 298°K conati- 
tutas a minimum temperature limit in the calsulations. During the 
frozen flow expansion to 0.2 psi, the temper ture dropped to this point 
in a few inetances. In these cases, the true frozen exhaust temper- 
ature lies somewhere below 298°K. Theze is no maximum temperature limit 
to the calculations, A linear heat capacity is autometically aasumed 
above the maximum temperature for which thermodynemic data are provided, 


if this temperature should be exceeded. 
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In addition to the compilation of data in Appendices C and D, 
the maximum shifting specific impulses of the various systems are sum- 
marized in six tables within this chapter. Each of the various oxi- 
dizers is listed, in order of its performance, with each of the fuels 
in Tables V, VI, and VII. These are for exhaust pressures of 14.7 psi, 
2.0 psi, and 0.2 psi respectively. Conversely, in Tables VIII, IX, snd 
X, the various fuels are tabulated in order of their performance with 
each oxidizer. The combustion temperature at the point of maximum 
apecific impulse is also listed. These tables provide a concise pre« 


sentation of the relative energy available from the various systems. 





Following theses tables is a brief discuasion of the more significant 


findings of this study. The discussion is subdivided by fuels, 
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Table V 


Summary of Maximum Shifting Specific Impulse 
pe = 1000 psi 
Pe = 14.696 psi 





{ 
Wt. % Shifting Isp, 

Fuel Oxidizer Oxidizer Seconds Tc, °K 
Hp OF 2 85.5 a 3591 
Ho Fo &9 410 3964 
Ho Op 8 391 2769 
Ho NF3 93 351 38°76 
Ho NO5C10;, 83 349 2713 
Ho C103F 8h 344 a7uk 
Ho WoO), 84 342 2660 
Ho Ho00 88 322 240) 
Ho ClF3 92 318 34.03 
Ho NH,C10), 91 287 euus 
NoHy Fo 69 363 4688 
Non, OF 2 60 345 4037 
NoHy NF3 3 322 haha 
Nowy 02 48 333 3409 
Now, NO2C10}, 56 295 3363 
Nol, C103F 59.5 295 3466 
Nom, C1F3 73 293 3882 
Nou, N20), 5T 291 3257 
NoH, Hp02 67 286 2923 
NoH, NHi,C10), 3 265 2969 
CHo OFo 19 350 4650 
CHo Fo 72.5 325 4.382 
Co 02 3 300 3683 
CHa NF 78 287 3781 
CHo cid F 81 280 3722 
CHo Hope} 79 278 3586 
CH H50a 86.5 2T7 2999 
GHe N20), 80 276 uhh 
CHo clr val 257 3541 
CHo NH,C10), 90.5 252 3016 
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Fuel Oxidizer 
UDMH OF 5 
UDMH Fo 

UDMH O02 

UDMH NF 

UDMH Clb sy 
UDMH NOpC10,, 
UDMH No0h 
UDMH Ho0e 
UDMH C1F3 
UDMH NH,.C1LOy, 
Ci,No Op 

CyNe No0k 
CiNo H202 
Ci,No OF 2 
Ch.No NO2C10), 
CyN5 N,C1O}, 
ChLNo C103F 
Che Fo 

ChNo NF3 
Chia clF3 

BS OFs 

Bs Fo 

Bs 02 

BSH NF3 
B5H9 

B5Hg C103F 
55 N20), 

Bs NOsC1O, 
BS one 
BS NH,010, 
MgHo F2 

S & 
MgHo Hobe 
MgHo C1F3 
MgHo C104, 
MgEo No 

NgHp NO2C10}, 
MgHo C103F 


Table V Continued 


Wt. % 


Oxidizer 


Shifting Isp, 
Seconds 


351 
3hh 
310 
309 
290 
289 
286 
283 
280 
259 


302 
282 
280 
278 
2°76 
258 
255 
253 
232 
203 


367 
360 
327 
326 
316 


275 


260 








Oxidizer 


Fo 
OFo 
He0e 
M3 


02 

Noo, 
No2C1LO}), 
NH,C1O, 
CLO3F 
C1F3 


Fo 
OF2 
Ho02 
Op 


NF 

N20), 
NO2C10, 
NH,C10, 
C103F 


Table V Continued 


Wt. & 


Oxidizer 


136 


29 
48 
78 


42 


Shifting Isp, 
Seconds 


393 
333 
326 
323 
318 
308 
305 
302 
301 
291 


395 
3€3 
315 
371 
359 
351 
346 
340 
339 
329 


363 
337 
319 
293 
272 
263 
262 
258 
250 
250 


Tc, CK 


5317 
LhO6 
3923 
4.956 
4299 
226 
4175 
3837 
ho29 
4593 


9331 
4347 
3574 
4338 
4731 
3403 


3283 
4173 
4500 


4886 
4093 
1,364 
4109 
3351 
3309 
2711 
3292 
2754 
3172 











Table VI 


Summary of Maximum Shifting Specific Impulse 
pe * 1000 psi 
Pe = 2.0 psi 





Wt. % Shifting Isp, 
Fuel Oxidizer Oxidizer Seconds Tc, °K 
Ho OF 5 87 458 3756 
Ho Fo 92 457 462 
Hp 0 82 437 3182 
He wf; 94.5 390 4212 
Ho NOSC10, 86 389 3080 
Hp C103F 87 383 3273 
Hp NoO}, 86 381 2892 
Hp He0p 91 359 27k0 
Ho clr 93.5 352 3690 
Ho NH,010), 92.5 319 2660 
Nom, Fo 69 408 4.688 
NOR OF 2 63 392 4100 
NOH, NF3 73 357 yoke 
Non, hg 355 3425 
Now, NO2C10, 58 333 3375 
Noh, C103F 61 333 3478 
Nem, N20), 59 329 3262 
Nok), ClF3 a 325 3862 
NoRi, 68 323 2916 
NOR), NB,C10, 73 297 2969 
CHe CF, 19 393 4650 
CH Op 73 344 3683 
Cp NF Wis) 326 3781 
CHp cl 82 319 3734 
CHo NO2C10), 81 318 3603 
CEs NoO), 82 316 3453 
CHp Ho 88 325 2979 
CHp clF T7 291 3541 
CHo 1,210, 91 287 3002 
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Fuel 


UDMH 


C.No 


Oxidizer 


OFo 

Fo 

Oo 

NF 
C1OaF 
NO2710), 
NoOL, 
Ha02 
C1F3 
NH,C1Oy, 


Op 

No0y, 
OF 5 
B20. 
NOnC10), 
NHi,C10}, 
F 


Table VI Continued 


Wt. % 


Oxidizer 


123 


Shifting Isp, 
Seconds 


400 
393 
354 
348 
329 
328 
325 
321 
314 
293 


342 
320 
320 
319 
343 
293 
292 
289 
269 
235 


hoo 
413 
379 
374 
361 
353 
353 
350 
333 
330 


379 


3h6 
327 
316 
315 
311 
310 


305 





fa, CX 


4519 
43h2 
3633 


3686 
3567 
3433 


3794 
3026 


4609 
4127 
4061 
3320 
4168 
3264 
3891 
3708 


3404 


5169 
5165 
4518 
4841 
3400 
447 
4388 
4502 
4487 
3000 
4991 
4LO} 
4620 
3073 
3216 
465 
3073 
3265 
3124 
3390 





Oxidizer 


C103F 


Op 
NOoC1Oy 
NHi,C1LOy 


Table VI Continued 


Wt. % 
Oxidizer 


1 
62 
48 


124 


Shifting Isp, 
Second 


405 
386 
378 
372 
370 
358 
355 
350 
350 
333 


453 
WAS 
Lhe 
436 
413 
412 
L06 
398 
396 
378 


410 
386 
360 
330 


To, °K 


5346 
LLL 
3923 
4299 
5006 
4226 
he52 
3837 
eok 
4693 


5331 
i137 
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Table VII 


Summary of Maximum Shifting Specific Impulse 
pe = 1000 psi 
Pe = 0.2 psi 


Wt. % Shifting Isp, 
Fuel Oxidizer Oxidizer Seconds Te, °K 
Hp OF 2 91 Kol Ti Fxe) 
Hp Fp oh 489 4809 
Ho Op 8h 170 3374 
Ho NO2C10), 87.5 416 3259 
Ho NF Oh.5 413 412 
Ho C1O3F 90 hog 3502 
Hp NoO), 89 4o8 325) 
Hp Ho02 92.5 385 agi2 
Ho C1F3 95 37 3985 
Hp NH,C10), oh 340 2883 
Non, Fo 69 4.36 4688 
NoHy OF 2 63 yok 43.00 
Noh, Oo 4g 364 3425 
NOH, NF 13 379 hohe 
Nol, NOSCLOY 58 359 3375 
NoHy, C103F 61 358 3478 
NoHy N20}, 59 354 3262 
NoHy Ho05 68 348 2916 
Non, ae 73 345 3882 
Now, NH,C10y, 3 318 2969 
CHo OF 2 19 420 4650 
CH Fo 7 4O5 4hOS 
CH Qo 79 379 3699 
CHo NF3 ¥ (3) 352 3761 
CH NO2C10), 82.5 349 3587 
CHo C103F 85 347 3700 
CHo NoO} 82 345 3453 
CHo Ho0p 88 344 2979 
CHo ing TT 315 3541 
CHo NH, C104 gi 309 3001 
125 





Fuel Oxidizer 
UDMH OF a 
UDMH Fo 
UDMH Oo 
UDMH NF 
UDME C102F 
UDMH Nost10), 
UDMH NoO;, 
UDMH 20. 
UDME c1F 
UDMH wH,020,, 
CiNe 

Cho OF 
CiNp No0} 
ChNo 7308 
CiNo NOoC10;, 
CiNo F 

Ci,No ¥f,C10), 
Ch.No C103F 
CN NF 

C1, No ci; 

Bs OFa 
=e Fo 
me 
Bi ro8o 
me ClO 
Be NOot10), 
3 CLF 

BS NHC1Ol, 
MgHo Fo 
MgHo OFo 
MgHe NF3 
Nigro Fave 
MgHo 02 
MgHo Nok 
MgHo NEYC1O}, 
MeHe cir 
MgHo H02C10) 
Melo C103F 





Wt. % 
Oxidiz 


81.5 


Table VII Continued 


er 
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Shifting Isp, 


Seconds 


435 
426 


4543 


4OO1 
4337 
4620 
3073 
3216 
3265 
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Table VII Continued 





Wt. & Shifting Isp, Ps 
Fuel Oxidizer Oxidizer Seconds Tc, K 

Fo TT 7 5368 
OF 62 431 Wh 7 
Ho05 48 yok 3923 
Cp yd 416 4374, 
NF 81 OT 5006 
N 3, 5h 4OL age 

C10, a1 39T 4252 
NH, C10}, 58 393 3837 
C103F 55 390 kek 
C1F3 80 367 4693 
Fo 19 501 5194 
HoQe 61 498 308k 
OF'o 70 497 4527 
Op 53 492 4004 
NoO), 52 462 3446 
NF3 81.5 456 4731 
NOZC1O, 64 456 4132 
NH,C1Oy, 62 46 3283 
C105F 68 bYG 4173 
C1F3 81.5 418 4500 
Fo 82 Wks 4886 
OF > 77 yok 4093 
vk 85.5 392 436) 
ClF3 8h 360 4109 
C103F 76 345 3354 
Op 66 342 3309 
H20. 81 341 271). 
NO2C10), 73 333 3292 
NAC1.0, 63.5 326 2292 
N20), 7TH 325 3172 
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Oxidizer 


Fuel 


Lif 


Table VIII 


| 
Summary of Maximum Shifting Specific Impulse 
Pe * 1000 psi 
De * 14.696 psi 
Wt. % Shifting Isp, 
Oxidizer Seconds To, OK 
89 410 3964 
80.5 395 5331 
82 363 
69 363 4688 
82 360 5101 
1365 353 5317 
Tl 34h 4342 
Te 329 4963 
72.5 325 4382 
85 253 3691 | 
78 391 2769 || 
54.5 371 4138 | 
70 327 4518 | 
ho 318 4og9 
48 313 34.09 
63 310 3614 
4g 302 4816 
73 300 3683 
31 270 3076 . 
66 263 3309 
85.5 hi 3591 i 
67 383 4347 iT 
80 367 5169 ; 
Te 351 4460 
79 350 4656 
60 345 4037 
TT 337 4093 
59 333 4207 | 
79 305 4516 
68 278 4061 at 


Ci,No 


FS ae Rk ee nee 
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#202 
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Table VIII Continued 


Wt. % 


Oxidizer 


56.5 
48 
88 
73 
67 
81 
39 
16) 
86.5 
81 
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Shifting Isp, 
Seconds 


375 
326 
322 
316 
286 
283 
280 
280 
277 
262 


352 
342 
308 
306 
291 
286 
282 
276 
265 
250 


359 
351 
326 
323 
322 
319 
309 
302 
287 
232 


329 
318 
"293 
293 
291 
250 
280 
275 
257 
203 


Te, 


3574 
3923 
2hoh 
3400 
2923 
3008 
3073 
3320 
2999 
2711 


3403 
2660 
4226 
426 

3257 
3429 
hook, 
lal 
3073 
3171 


4731 
38°76 
4809 
4956 
hohe 
430h 


964 
4620 
3781. 
3680 


4500 
3403 
4109 
3882 
4593 
L487 
3794 
465 
3541 
3hok 





tava DO 


Oxidizer 


C1O3F 
C1O3F 
C103F 
C103F 
C103F 
CLO3F 
C103F 
C1O3F 
C1O3F 


NOoC10), 
NO2C10), 
NOoC10), 
NO2C10L. 
NOpC1Oy, 
NO2C10), 
NO2C10L, 
NOoC10,, 
NOpC10, 
NOoC1LOy, 


NH.C10), 
NB}CLO, 
NH,C1O, 
NH,C1O, 
NH,C1LO} 
NH,C10), 
N,C1O}, 
NHC1O, 
NEC1O, 
NH,C10, 


Fuel 


Ho 
B 
ici 
Nowy, 
Liz 


MgHp 
Ch.No 


BsHg 


Table VIII Continued 


Wt. % Shifting Isp, 
Oxidizer Seconds 
8h 3hh 
68 339 
719 306 
50 301 
5905 295 
81 280 
6 272 
46 260 
62 255 
83 3h9 
55 346 
ho 305 
78 302 
56 295 
70 289 
79 278 
61 276 
38 263 
73 258 
62 340 
58 302 
91 287 
82 285 
50 266 
B 265 
85 259 
80 258 
90.5 252 
82 2hg 
130 


Te, °K 


27h 
4173 
ANT 
o29 
3466 
3722 
3351 
3237 
3708 


2713 
3394 
4175 
430 
3363 
3540 
3586 
4o72 
3124 
3292 


3283 
3837 
2hh8 
3484, 
3073 


3038 
3333 
3016 
2687 
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Table IX 


Summery of Maximum Shifting Specific Impulse 
Pe = 1000 psi 








Pe = 2.0 psi 
Wt. % Shifting Isp, é 
Oxidizer Fuel Oxidizer Seconds Tc, K 
Fo 92 4ST 462 
: F2 Bete ag W533 33° 
. F B 13 5165 
‘ Fe rae 62 ha ee 
my F WV 
x Fo rai 5 405 5346 
Fo UDME TL 393 43h2 
Fo MgHo TH 379 4991 
Fe CHp 7205 372 4.382 
Fo C),No 85 292 3891 
82 437 3182 
0g rete 53 436 404 
02 BSH 70 379 4518 
02 Aix yb 371 L209 
Op Ne kg 355 3425 
02 UDMH 65 354 3633 
02 CHe 73 344 3683 
Co Chto 52 342 ee 
MgHo 33 316 2 
02 Li 66 306 3309 
oF. bs! 87 458 3756 
ore Beko 67 Wis W347 
OF2 Bok 80 422 5169 
OF voMe 15 4.00 4519 
OF2 CHp 19 393 4650 
OF 2 NoHy 63 392 4100 
OF 5 Lik TT 386 093 
OFo AlH3 62 386 47 
OF MeHo 69 352 ok 
OF2 ChNp 68 320 4061 
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Oxidizer Fuel 


Ha05 Belo 
He0e ALH 
Hp0p B 
#200 m° 
Has Mate 
Ha05 Nou, 
Hoo UDMH 
F302 CyNo 
He02 CHo 
Hp0o Li 
N20) Bes 
NoO), Ho 
Nod), ALE 
No0, BsHg 
NoO), Noh), 
N20}, UDMH 
Na0h CyNo 
NO}. CH 
No0) MgHo 
NoOh Lil 
NF3 BeHo 
mS Bo 
WS pee 
NF3 Liz 
NF3 No, 
NF3 UDMH 
NF3 MgHo 
NF3 CHp 
NF3 Ch.No 
C1FS BeHo 
C1F3 Ho 
C1F3 BH 
C1F3 Lin 
C1F3 A1LH3 
C1F3 No, 
clr NgHo 
C1F3 UDMH 
C1F3 CHo 
C1F3 C,No 


Table IX 


Continued { 








Wt. % Shifting Isp, . i 
Oxidizer Seconds Te, K | 
58 4he 3708 | 
48 378 3923 | 
73 361 34.00 
91 359 2740 
42 327 3073 | 
68 323 2916 
81 321 3008 
19 319 3320 
88 315 2979 
81 30h 2711 
50 411 3403 
86 381 2891 
52 358 4226 
19 353 4388 
59 329 3262 
TH 325 3433 
64 320 4127 
82 316 353 
yh 310 3265 | 
Th 290 3171 
81.5 413 4731 
9h..5 390 4212 
88 374 4BH1 
81 370 5006 | 
85.5 360 4384 
73 357 ake | 
76 348 3964 
78 346 4620 
78 326 3781 
85 269 3680 
81.5 378 4500 
93.5 352 3690 
&8 333 HUST 
84 330 4109 
80 333 4693 
3 325 3882 
778 315 LESS 
16) 324 379% 
TT 291 3541 
87 235 34-04 
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Oxidizer 


C103F 
C103F 
CLO3F 
C1LO3F 
CLO3F 
C1LO3F 
C103F 
CLOSF 
C103F 
C1O3F. 


NOpC10), 
NOacioO), 
NOoC10), 
NOoC1O), 
NOoC1O), 
NOpC10b 
NO2C104 
NOoC10), 
NOeC10), 
NOoC1Oy, 


NH,C1LO), 
NACLO), 
NH,C1O), 
NH),C10), 
NHYC1LOY, 

C10, 
NHC1O, 

c1Oh 
NH),C1O, 
WH, C10), 


Fuel 


Beko 


B 
ict 
Nowy 
UDMH 


CHa 
Lik 


CyNo 
BeHo 
Hp 
ALH3 
Bs 
UDME 
CHp 
C,No 
Liz 
BeHo 


BsHo 


Nol, 
UDME 
CNo 
Li 


Table IX Continued. i 
| 


Wt. % Shifting Isp, 

Oxidizer Seconds Tc, OK 
68 396 4173 
87 383 3273 
19 353 Yhh7 

5 350 4204 
61 333 3478 
5 329 3686 
82 319 3734 
76 312 3351 
hg 305 3390 
62 289 3708 

2 4O6 3866 
86 389 3080 
a1 355 ho5e 
19 350 4502 
58 333 3375 
3 328 3567 
81 318 3603 
64 313 4168 
ho 308 3124 
73 299 3292 
62 398 3283 
58 350 3837 
%3 330 3000 
92.5 319 2660 
50 311. 3073 
73 297 2989 
86 293 3026 
83 293 3264 
82 288 2687 
91 237 3001 
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Table X | 
Summary of Maximum Shifting Specific Impulse 
Pe = 1000 psi ‘ 
Pe = 062 psi | 
Wt. % Shifting Isp, 
Oxidizer Fuel Oxidizer Seconds Tc, °K 
Fo BeHo 73 501 5194 
Fo Fp oh 489 4809 
Fo B5Ho 8h 460 5165 
Fo aan TT Wher 5368 
Fo Noh, 69 436 4688 | 
Fo UDME 72 426 3h2 
Fo MgHp Ts 420 4991 
Fo CHo Tf 405 4LOS |! 
Fo CiNo 8 325 3891 | 
Og BeHp 53 492 4004 fe 
02 Hp 8h 470 3374 ya 
C2 te 72 heL 4620 
Op 3 Ih 416 4374 
Op UDMH 68 388 3626 
0g Now, 4g 384 3425 
02 CH 73 379 3699 
02 CyNo 55 37h 75 
7) MgHo 33 357 3216 
Oo Lil 66 342 3309 
OF2 BeHo 70 LOT 4527 
OF a Ho 91 491 hik9 
OF p noe 81.5 466 5215 
OF 73 435 45210 
OF's ALE. 62 437 Ahh 
OF wh 63 kok 4100 
OF Lik TT 4oh 4093 
OFo CHp 79 420 4650 
OF 5 Melo 65 394 4337 
OF Cho 65 355 4112 
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Table X Continued { 





Wt. % Shifting Isp, 3 
Oxidizer Fuel Oxidizer Seconds Tc, K { 
H205 Be 61 498 3884 
Hp05 AE, 48 hol, 3923 | 
B20 BsHg 73 399 34.00 
Ha0p HE 92.5 385 2912 
HpGp =e ko 368 3073 
505 C,N5 78 350 3257 
Be0a UDMH 82 349 2992 
H502 No, 68 348 2916 
B05 CHp 88 Bh 2979 
Hp00 Lik 81 341 2711 
NA0 Be 52 462 3446 
W20) Bp 89 408 325) 
NoO}, Ali 5h 402 4292 
No0) 35 79 391 4388 
No0), p 59 354 3262 
NoOy, UDMH TH 354 3433 
NO}, CyNo 67 350 03} 
N20, Malo dlp 3h9 3265 
NoO, CHo 88 344 2979 
No0, Lik TH 325 3171 
IF 3 BeHs 81.5 456 4732 
NF3 Hp O45 413 4212 
IF3 B5Hg 89 413 4854 
NF3 AnH 81 OT 5006 
NF3 Li 85.5 392 438) 
IF3 Melis 78 380 4620 
NF3 NoHi, 73 379 hale 
NF3 UDMH 76 373 3964 
NF3 CH 78 352 3781 
IF3 ChNo 85 301 3680 
C1F3 Belo 81.5 418 4500 
ciF Ep 95 373 3985 
C1F3 B 89 368 4543 
ClF3 3 367 4693 
cir Lif 84 360 4109 
clr 7 347 465 
C1F3 NOH, 73 345 2969 
C1F3 UDMH ie) 337 3794 
CiF3 Chie 7 315 35h2 
C1F3 CN. 87 261 34.04 
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Oxidizer 


C1O 

closF 
Cl02F 
C103F 
C10 

C103F 
C1LO3F 
C103F 
c10 

CLO3F 


NOoC10);, 
NOoC101, 
NOoC10), 
NOoC1OL, 
NOoC10), 
NOoC10), 
NOoC10}, 
NO2C10), 
NO2C10), 
NOoC10, 


NA,C10L, 
NH,C 

NH,C1O), 
NH, C10), 
NH,C1O;, 
NH,C10}, 
NH,C1O4 
NR,CLO), 
NH,C10, 
NH,C101, 


C1,.No 


Table X Continued 


TH, % 


Oxidizer 


136 


Shifting Isp, 


Seconds 


4h6 
409 
390 
390 


—_——— 


at 


B. Discussion of the Results 


1. Liquid Hydrogen Systems 
Both liquid flv-rine and liquid oxygen difluoride were found to be 


excellent oxidizers for liquid hydrogen, with the latter yielding slightly 
higher maximum performance. On the fuel-rich side of maximum performance, 
oxygen difluoride is the superior oxidizer. However, the performance 

drops ~ off more gradually on the oxidizer~rich side using liquid fluorine, 
Combustion temperatures are comparabie in both systems up to the point of 


maximam specific impulse. Beyond that, on the oxidizer-rich side, the 


formation of HF, with the formetion of Ke0 also contributing in the Hp~ 
OF> system. It is interesting to note that the maximum performance | 
occurs on the fuel-rich side of stoichiometric combustion, indicating the 
value of hydrogen gas as a working finid. In fact, disscciation is re- 
latively slight in both systems on the fuel-rich side of peak performance, 
causing only a slight difference in the frozen and shifting specific 


Ho-Fo system burns with the hotter flame. The prime energy source is the 
impulse. However, as the syatems become oxidizer-rich, the degree of | 


ree ec 


chamber dissociation increases rapidly, end the spread between frozen and 
shifting performance increases accordingly. 

Liquid oxygen rates closely behind oxygen difluoride and liquid fluorine 
as an energetic oxidizer for liquid hydrogen. Again, the maximum performance 


occurs on the fuel-rich side of stoichiametric combustion. Of course, the 





principal source of energy is the formation of water. Dissociation, and 


the corresponding spread between shifting and frozen impulse, increases 








as the system becomes oxidizer-rich. The most significant difference 
between the Ho-0> system and previous two is the lower combustion tempera- 


ture of Ho-02. Although somewhat higher performance may be achieved by 


i el 


oxidizing hydrogen to HF rather than to Ho0, a considerably hotter flame 
mist be tolerated. This characteristic is also displayed by the other less- 
energetic oxidizers. in which nitrogen and chlorine serve as carrier 


= etoma for oxygen and fluorine. 





Next in order of performance are NF 3 NO,C10),, C103F and N04 These 
four oxidizera are closely grouped, yielding maximum apecific impulses of 


351, 349. 34h, and 341 seconds respectively at an exhaust pressure of one 


to five seconds at an exhaust pressure of 0.2 psi. The halogens, chlorine 


and fluorine, appear as HCL and HF in the products of combustion. Both 


constitute energy sources, although HCl is considerably less vaiuable in 


| a. adit bol sacslec es aka eee 





atmosphere. Furthermore, the shifting specific impulse spread is narrowed 
this respect. Aiso, its higher moleculer veight mkes it relatively less | 


desirable as a working fluid. In each system, the peak performance still | 
oceurs on the fuel-rich side of stoichiometric combustion. This further 


illustrates the value of mly partially oxidizing the hydroger in order 


system, moet of the nitrogen appeers as N» in the combustion products. 
However, upon reacting fuel-rich, considerable amunts of NH3 are formed | 


to obtain the maximum useful energy. In the H5-NO2C10) and the Ho-No0), | 
at the lower exhaust pressures. The lower pressures, coupled with the | 





corresponding lower temperatures, promote the formation of this larger | 
molecule. In each case, the spread between frozen and shifting specific | 


impulse is considerably greater on the oxidizer-rich side of the performance 
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curve. As before, this is attributed to the greater degree of dissocia~ 
tion at the higher combustion temperatures. 

Hydrogen peroxide and chlorine trifluoride constitute the next oxidizer 
performance plateau. Chemicaliy, the H)-Ho0> system differs from the Ho~ 
Oo system only in that the oxygen is partially reacted with hydrogen in 
Fo02. Whereas the heat-of-formation of liquid oxygen is -3.08 Keal/mole 
at ite boiling point, that of Ho 4, 48h Keal/mole at 298°K. This 
means that a substantial amount of useful energy has been lost in Khe. 
The loss in shifting specific impulse is 69 seconds at an exhaust pressure 
of 1 atmosphere, 78 seconds at 2.0 psi, and 85 seconds at 0.2 psi, 
Similarly, ClF; differs from NF; in that chlorine has replaced nitrogen 
asa carrier for fluorine. The more exothermic heat-of-formation for C1F3 
indicates that the Cl-F bond is stronger, and therefore less desirable 
in a propellant ingredient. The loss in maximum specific impulse is 33 
seconds when exhausting to one atmosphere, 38 seconds at 2.0 pai, and 40 
seconds at 0.2 psi. The lower molecular weight of nitrogen also contri- 
butes to these performaace spreads. 

Ammonium perchlorate was found to te the least energetic oxidizer 
for liquid hydrogen. This was to be expected in view of its relatively 
exothermic heat-of-formation of -69.\2 Keal/mole and the presence of both 
nitrogen and chlorine as carrier elements. The formation of Ho0 and HCl 
constitute the principal energy sources, On the other hand, Ho aad No are 
the principal unreacted working fluids. However, again iarge concentra- 
tions of NH, are obtained at the sub-atmospheric exhaust pressures when the 


stoichiometry is fuel-rich. The C10), moiety is common to both ammonium 
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perchlorate and nitronium perchlorate. However, the marked superiority 
of NOpC1L0), as an oxidizer is quite apparent from the performance curves. 
This stems from the substantially more endothermic heat-of-formtion of 
NO5C10), as well as the greater quantity of available oxygen. Ammonium 
perchlorate, like hydrogen peroxide, may be thought of as low-energy 
monopropellant since both hydrogen and oxygen are available in the same 
molecule for subsequent chemical reaction. 
2. Hydrazine Systems 

Using hydrazine as the rocket fuel, liquid fluorine is the most energetic 
exidizer. Although Nou), has an endothermic heat-of-formation, there is a 
substantial drop in performance caused by the introduction of nitrogen as 
a carrier element. Whereas the Ho-Fo system yields a maximm Isp of 410 
seconds at pe/pe = 68, the NoH,-Fo system yields only 363 seconds. Chemi- 
cally, the NoHy,-Fo system is comparable to the Ho-NF3 system. Of these, 
the Noi) -Fo system is superior because of the endothermic heat-of-formation 
of hydrazine. This points mt the desirability of utilizing the most en- 
dothermic reactants for propellant ingredients, providing other considera- 
tions are equal. Furthermore, the endothermic heat-of-formation of hydra- 
zine causes a higher flame temperature than using liquid hydrogen as the 
fuel. Subsequently, dissociation is greater and the spread between frozen 
and shifting Isp is greater, on a percentage basis. These generalizations 
apply to the use of each oxidizer with hydrazine. 

Oxygen difluoride was found to be the second most energetic oxidizer 
for use with hydrazine. The principal exhaust products for this system are 


BF, 10, Hp, and No. Next in order of performance are nitrogen trifluoride 
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and liquid oxygen. Like the He-NF3 and the NoH})-Fo systems. the NoH)-N¥3 


system is composed of only nitrogen, hydrogen. and fluorine. However, in 


this latter system nitrogen serves as a carrier for both hydrogen ani fluorine. 


Conseauently, the performance is lower. The still lower performance of NoH)~ 


O25 relative to NoH)-NF3, indicates that the presence of nitrogen in the 
fuel is more detrimental when oxidizing with oxygen than when utilizing 
fluorine. This trend is mintained with the other oxidizers considered. 
Both Cl03F and C13, the other fluorine<containing oxidizers, ranked re~ 
latively more potent with hydrazine then with liquid hydrogen although 
the absolute performances were lower. Very little amuonia was present in 
the combustion products of the hydrazine-oxygen system, even on the fuel~ 
rich side, despite the presence of large quantities of nitrogen and hydrogen 
gases. This was due to the higher exhaust temperatures in this system. 
In general, NHj is a significant exhaust specie only below 600°K at a 
pressure of 2-0 psi and only below 500°K at a pressure of 0.2 psi. 
Perchloryl fluoride and nitronium perchlorate yielded equivalent per- 
formances with hydrazine, The presence of fluorine in C103F offset the 
more endothermic heat-of-formation of NOoC10,. With liquid hydrogen, the 
latter oxidizer was superior by 5-7 seconds, Furthermore, the presence 
of fluorine in C103¥ accounted for a slightly higher combustion temperature, 
agein offsetting the higher heat-of-formation of NOsC1C,. Neither system 
was very bighly dissociated, even on the oxygen-rich sides. The maximum 


combustion temperatures occurred near the stoichiometry of maximum specific 


impulse. 
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Closely following perchloryl fluoride and nitronium perchlorate in 
order of performance are chiorine trifluoride, nitrogen tetroxide, and 
hydrogen peroxide respectively. The most noticeable contrast in the 
performance curves of the three systems is the sharp drop-off on the 
oxidizer-rich side of the Clg system. In general, this is typical of 
the behavior of C1F 3 systems and, to a lesser extent, NF3 systems. How- 
ever, this is more pronounced with some cf the other fuels, as will be 
shown later. Ammonium perchlorate is the least energetic oxidizer, for 
the same reasons cited previously. 

3 Polyethylene Systems 

Polyethylene, with the empirical forma (CHo)x, was chosen as the 
prototype fuel to represent the hydrocarbon family. Having the highest 
H/G ratio of any of ‘the aliphatic hydrocarbons, CHo represents the maximm 
energy source from this family. Furthermore, its stoichiometry corresponds 
to that of diethylcyclohexane, CycHod. elthough the performance of the 
latter is somewhat less because of its slightly more exothermic heat~of-= 
formation. 

Oxygen difluoride was found to be the best oxidizer for CHp. One of 
the prime reasons for this is the ideal stoichiometry for the formation of 
HF and CO as products of combustion. At the point of maximum specific 
impulse, the HF is slightly dissociated into H, F, and Ho because of the 
high combustion temperature. In the CHo~0Fp system, greater performance is 
obtained by forming the maximum amount of HF and CO than by under-oxidi-~ 
zing in order to obtain hydrogen. By under-oxidizing, large quantities c* 


solid carton are formed in the exhaust. Over-oxidizing tends to promote 
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the formation of carbon dioxide as well es allowing F and Fo to emerge 
uncombined. ‘he spread between frozen and shifting impulse reaches a 
maximum at the point of maximm performance, 

Next in order of performance is liguid fluorine. In this case, HF 
alone constitutes energy source with all the carbon emerging as the un- 
reacted solid element. This further emphasizes the value of HF as a 
product of combustion. On the oxidizer-rich side of the system, solid 
carbon is still the predominant carbon-containing specie. Only small 
quantities of fluorocarbons were present in the exhaust. Liquid oxygen, 
the next best oxidizer, is distinctly inferior to either OF> or Fo. In 
the CHo=02 system, the predominant chemical equilibrium is the water-gas 
shift reaction. During expansion, some of the carbon monoxide and water 
react to form carbon dioxide and hydrogen. Under-oxidizing reduces per~ 
formance because CO and COo do not supply sufficient energy to the low 
molecular weight hydrogen gas. On the other hand, over-oxidizing lowers 
performance because too mech C0 and Hed are present. Thus, more energy 
may be obtained from a hydrocarbon by its combustion with a combination of 
oxygen and fluorine than with either oxygen or fluorine alone. 

Nitrogen trifluoride, perchloryl fluoride, nitronium perchlorate, 
hydrogen peroxide, and nitrogen tetroxide ~re closely grouped in the next 
energy level. At pc/pe = 68 the performance spread for this group is 11 
seconds while at pe/pe = 5000 it is only 7 seconds. In view of this 
relatively small spread in specific impulse, the considerably higher density 
of NOpC10), makes it the best of the group for overall missile performance, 


Chemically, the C103F-CHo system is the most complex with HF, M1, Hed, COs, 
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and CO all contributing useful energy. Completing the ranking of the 
oxidizers are C1F3 and NH), Cl0,. Maximim performance occurs in the C1F3- 
CHp system near the point of mximm formation of HF and HCl. However, the 
large quantity of unreacted carbon makes the system unattractive from an 
energy standpoint. 

With the exception of OFg, every oxidizer behaves poorer with CH than 
with NoH)- Of course, the performance ‘s still lower than when reacting 
each oxidizer with liquid hydrogen. Thus, carbon appears less desirable 
than either hydrogen or nitrogen as a propellant ingredient, fram the 
energy viewpoint. Also, oxygen appears to be essential for the energetic 
oxidétion of carbon. 

h. DME Systems 

The use of unsymmetrical dimethyl hydrazine as a fuel introduces carbon, 
nitrogen, and hydrogen simitaneously into a propellant system. Although 
two of hydrazine's four hydrogen atoms are replaced by methyl groups in 
UDMH, there is less than 1 Keal difference in the heat-of-formation of the 
two molecules. Thus, the relative performances of hydrazine and UDMH pro-- 
vide a good measure of the effect of carbon on performance, 

Again, OFo is the hest oxidizer with Fo following closely. Liquid oxygen 
ranks third, but considerably below the first two. In fact, the UDMH-0O- 
system exceeds the UDMH-NF3 system by only one second at pc/pe = 68. Howevery 
at pe/pe = 5000 this spread is 15 seconds. This further shows that the 
presence of nitrogen is more detrimental to the C~-H-O system than to the C. 
H-F system. In general, the same comments made previously on stoichiometry 


and reaction products apply here. In a C-H-O-N system, maximum performance 
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usually lies somewhere between balancing the reactants to CO and He0 and to 
COp and Hod. In » C-H-F system, the greatest energy release occurs when 
all the hydrogen and fluorine combine to form HF, 

The other oxidizers, in order of their performunce with UDMH, are 
perchloryl fluoride, nitronium perchlorate, nitrogen tetroxide, hydrogen 
peroxide, chlorine trifluoride, and afmonium perchlorate. It is interesting 
to note that the performance of the various UDMH systems is only 3-6 
seconds lower than that of corresponding hydrazine systems at pc/pe = 68 
for any oxidizer containing oxygen, At pc/pe=5000, the difference is even 
less. In fact, the performance of UDMH-0> exceeds that of NeH)-02 at pe/pe 
= 5000. The same holds when NFz, Fo» and CiF3 are the oxidizers. Although 
the Igp spread between UDMH and hydrazine varies from 13-19 seconds at 
pe/pe = 68, the variation is only 6«10 seconds at pce/ne = 5000. Chamber 
temperatures ars slightly higher when UDMH is reacted with oxygen-contain- 
ing oxidizers than when hydrazine is the fuel. The opposite is true when 
fluorine is th@ sole oxidizer. Although the presence of carbon in UDMH does 
not appreciably reduce the performance below that attainable from hydra- 
zine, two compensating factors mist be considered, One is the smail change 
in heat-of-formation of the fuel. Usually, the addition of alkyl groups 
to a fuel mkes its heat-of-formation more exothermic. The secend factor 
is the greater percentage of hydrogen in UDMH. This tends to offset the 
deleterious effects of carbon. 

In the UDMH-Fo, UDMH-NF3 9 and UDMH-C1F3 systems, practically all of the 
carbon emerged as the unreacted solid element. Even when operating oxidizer- 
rich, only small concentrations of flucrocarbons appeared. Of these, Cry 


was the predominant one. This increased in concentration as both temperature 


145 














sssstoomsegel 








and pressure decreased, but was always appreciably less than the concentra- 
tion of solid carbon. The fluorocarbons appear to be very poor products 


of combustion, both as an energy source and as a working fluid. 


5. Dicyanoacetylene Systems 


ee 


Dicyanoacetylene constitutes a unique prototype fuel in that it has an 
unusually endothermic heat-of-formation and is composed of only carbon and 
nitrogen. Thus, eight of the ten CyNo systems are void of hydrogen. With 
029 Ho02, Ng0)» NO2C10;,. and NH)C10), the performance of CyNp was found 
comparable to that of CHs. However, the combustion temperatures were 
considerably hotter because of the relatively slight dissociation of ths 
combustion products. When hydrogen wee absent, CO and No were the principal 
combustion products, and the maximum specific impulse occurred when re- 
acting slightly oxygen-rich of this stoichiometry. The carbon nitrides CN 
and CoNo were formed in only minor concentrations, being practically non- 
existent except on the fuel-rich side of the CyNo-No0), system. 

The presence of fluorine in these hydrogen-free systems caused a marked 
drop in performance. This was true even for OF> and C103F in which oxygen 
also is present. Using Fo, NF3_ and ClF3 48 oxidizers, the fluorocarbons 
were the principal energy sources. In the chamber CF3 was most predominant 
while CF), became increasingly more important es the temperature and pressure 
were reduced. Although combustion temperatures in excess of 3000°K were 
obtained, the performance was very poor. Thus, the CyNo series emphasizes 
that oxygen is essential for the energetic oxidation of carbon ard that a 
large endothermic heat-of-formtion cannot surpass the value of hydrogen 


in a propellant formation. 
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6s Pentaborane Systems 
Pentaborane, one of five light metal hydrides selected for this 


study, represents the boron hydride family. Both oxygen difluoride and 
liquid fluorine were found to be excellent oxidizers for pentaborane. 
The former system yields a specific impulse of 367 seconds at an exhaust 


pressure of one atmosphere while the latter yields 360 seconds under the 





same conditions. However, the performance of the OF - B5Hy system is 
contingent upon the formation of BOF as a product of combustion. Although 
this oxyhalide has never been isolated, the trimer has been synthesized 
and its heat-of -formation evaluated. Because of its relatively simple 
structure, the probability cf BOF occurring in a high temperature equi- 
libruim mixture is good. A heat-of-formation of =-150 Keal/mole at 298°K 
was used in this study. If BOF is not considered as a product of com- 
bustion, the maximum specific impulse of the OFo-B5Ho drops to 350 seconds. 
This peak occurs at 77% by weight oxidizer and still constitutes a very 
energetic system. Both oxygen difluoride and liquid fluorine produce 
combustion temperatures in excess of 5000°K at their points of maximum 
specific impulse. Because of these high temperatures, both systems are 
highly dissociated. Yor the Fa-B5Ho system, HF and AF are the principal 
species in the chamber. During expansion through the nozzle, BF3 and Hp 
also become major products of combustion. This is particularly true at the 
sub-atmospheric exhaust pressures. 

For the OFo-BsHg system, BOF and HF are the principal combustion 


products, with hydrogen gas increasing in concentration during expansion. 
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Both liquid oxygen and nitrogen trifluoride are alse good oxidizers 
for pentaborane. Closely following is hydrogen peroxide. Although some- 
whet less energetic, the Ho02-B5Hg system has a combustion temperature 
approximately 1100°K less than that of 02-BsHg and approximately 1400°K 
less than that of NF3-B5Hg- The combustion products of the latter system 
were egsentislly the seme as those of the F2-B5H9 system, except nitrogen 
gags was also present. Metaboric acid, HBOp, was the principal boron- 
vontaining specie not only in the 02-B5Ho and Hp02-BsHg systems but also 
in the others containing boron, oxygen, and hydrogen in the absence of 
fluorine. Liquid boric oxide becomes a major specie only when operating 
fuel-rich at lower temperatures. 

Nitrogen tetroxide, perchloryl fluoride, and nitronuim perchlorate 
constitute the next performance plateau with chlorine trifluoride and 
ammonuim perchlorate being the least energetic. Despite the strong affinity 
of boron for fluorine, chlorine trifluoride is not a particularly potent 
oxidizer for the boron hydrides. This constitutes still further evidence 
of the undesirability of the C1-F bond in propellant ingredients relative 
to the N-F or O-F bond. ‘The oxidizers No0),, NO2C10,, and NH\C10} possess 
gredually-sloping performance curves on the fuel-rich side of maximum 
performance. This may be attributed to the formation of crystalline 
boron nitride and liquid boric oxide under these conditions. 

When oxidizing pentaborane with Fo, N¥3, or CiF3, the maximum specific 
impulse occurs near the stoichiometry for which all the boron and fluorine 
exist as BF3. When oxygen is the only oxidizer, the stoichiometric reaction 


of all the boron and oxygen to HBOs will define the peak performance very 
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closely. In the OFo-BsHo system, maximum performance occurs near the 
point of stoichiometric reaction to HF and BOF. 

Whereas fluorine has a very strong affinity for boron, that of 
chlorine is not as great. Large concentrations of the boron chlorides 
are rare, with hydrogen chloride being formed in preference. Also, the 
borcn interhalides are not formed in large quantities relative to the 
boron fluorides. 


Te Aluminum Hydride Systems 


In general, the performance of aluminum hydride with the various 
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oxidizers ciosely ,arallels that of pentaborane. The most noticeable 
differences are OFa, Oo, HoOe, and NY\C1O,. Whereas the first two 
oxidizers are superior with pentaborane, the latter two yields higher 
performances with AlH3. Combustion temperatures are also comparable, 
with those for the AlH3 systems tending to be slightly higher. 

The principal energy source of the alvuminum-oxygen systems is con- 
densed Alp03, which has a highly exothermic heat-of-formation. Hydrogen 
gas is the main constituent of the working fluid. Using 02, Ho05, Noy, 


and NH,C10, as oxidizers for AlLH3, the maximum specific impulse occurs 





near the stoichiometry for which a?l the aluminum and oxygen react to 
Al203. The peak performance of the NO2C10),~AlH3 system is defined by 
reacting stoichiometrically to Ala03 and AlCl. 


Using fluorine as the oxidizer for aluminun, the point of maximm 





performance is more difficult to predict. With Fo, NF; » and C1F3, the 
peak performance occurs between reacting to AlF3 and AF, pius HF. The 


reaction products are usually highly dissociated, with AIF, AlFo, and AIF; 
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all existing in appreciable quantities. The first of these rredominates 
in the chamber while the last is most significant at exhaust conditions. 
Of course, HF also exists in large concentrations, particularly in the 
cu. ber. During e "pansion, the quantity of HF decreases in favor of 
AIF3. In the ALH3-C1F3 system, the aluminum chlorides and mixed halides 
are formed. The oxyhalide AlOF is a minor specie and does not possess the 
same importance of its boron counterpart, BOF. 

As in the case of boron, fluorine appears superior to oxygen as an 


oxidizer for aluminum. Liquid fluorine ranks as the best oxidizer for 





aluminum hydride, followed by OF2, Ho0o, NF3 » and Op. Closely grouped 
behind these oxidizers are NoO,, NOocC10),, ClO3F, and NH,C10,. It is 
interesting to note that ammonium perchiorate is almost as good an oxidizer 
as nitronium parctiiaeate and nitrogen tetroxide with Al1H3. This is 
probably due, in part, to the relatively small quantities of hydrogen 
available for working fluid in the latter two systems. Since A1203 is 


a relatively large molecule with a molecular weight of 102, 4 considerable 
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amonnt of iow molecular weight gas must accompany its formation to obtain 


high performance. Thus, if another source of working fluid, such as Ho 





or Noo), were added to the AlH3-No0) and A1H3-NOsC10, systems, the per- 


formance could probably be increased. The fact that hydrogen peroxide 
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excels over liquid oxygen as an oxidizer for AlH3 further emphasizes 





the importance of having large quantities of working fluid to accompany 
the formation of highly exothermic condensed phases, such as Al203. 
The spread between frozen and shifting specific impulse is con- 


siderably greater in the aluminum-fluorine systems than in the aluminum- 
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oxygen systems. Upon expanding to sub-atmospheric exhaust pressures, 
this spread becomes greater. 
8. Beryllium Hydride Systems 

Beryllium hydride was found to be an extremely energetic fuel. 
Only liquid hydrogen yielded greater performance at pc/pe = 68 with 
the oxidizers Fo, Oo, OFo, and NOeC1l0,. However, at pc/pe = 5000 
beryllium hydride was the most energetic fuel in every case. Under 
these conditions the BeHe-Fo system yielded a maximum specific impulse 
of 501 seconds, exceeding that of the Ho-Fo system by l2 seconds. 

Using fluorine as the oxidizer for beryllium hydride, the principal 
combustion products are BeFp and Ho. Because of tie high combustion 
temperatures, a significant degree of dissociation exists in the chamber. 
Hydrogen atoms, beryllium gas, and hydrogen fluoride are each present in 
fairly large amounts. Upon expansion, these shift to form BeFo and Ho, 
releasing considerable energy in the process. The maximum performance 
occurs near the point at which all beryllium anid fluorine react to form 
BeFo. Liquid fluorine, nitrogen trifluoride, and chlorine trifluoride 
each follow this pattern. 

With oxygen as the oxidizer for beryllium hydride, the principal 
combustion products are condensed beryllium oxide and hydrogen. As in 
the case of the alvminum-oxygen system, large quantities of working fluid 
must accompany the formation of the highly exothermic BeO in order to 
realize the optimum performance. This is illustrated by the fact that 
hydrogen peroxide yields greater performance than liquid oxygen with 


beryllium hydride, and ammonium perchlorate is almost as energetic as 
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nitronium perchlorate. In the beryllium-oxygen system, the maximun 
specific impulse wus found to occur on the fuel-rich side of stoichio- 
metric combustion to BeQ. This varied from 5% by weight oxidizer in the 
case of the BeHo-Ho02 system to 18% in the case of the BeHp-No0, system. 
The beryllium systems are different from the other metal hydride systems 


in that feirly significant quantities of the metallic gas exists at the 


point of maximum specific impulse. Like hydrogen gas, the beryllium gas 
serves eas additional working fluid to which the energy of formation of 
BeO may be transferred. Furthermore, when nitrogen was present in the 
oxidizer, significant quantities of solid beryllium nitride vere present 
at the point of maximum specific impulse. The affinity of beryllium 
for chlorine was also found to be quite high. In the chamber, BeCl 
replaces HCl as the predominant chloride, particularly at high tempera- 
tures. In the exhaust, BeClo is the principal chloride. ‘The gaseous 
polymers of the beryllium oxide were found to exist in relatively small 
quantities, having only a very slight, if not negligible, influence on 


performance. The principal beryilium-oxygen gas specie is beryllium | 


hydroxide, BeOH. Another feature of the beryllium systems is the relatively 
higher performance available upon expanding ic low pressures. The gain 

in performance over the other light metel systems is greater at the lower 
exhaust preesures. 


ulke the other chemical syatems studied, the spread between frozen 





and shifting specific impulse is considerably greater when carrier elements, 
such as nitrogen and chlorine, are not present. The spread varies from as 
much es 40 seconds in the BeHo-Fo system to as little as 10 seconds in the 


BeHo~NH},C10), system. F 
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9. Lithium Bydride Systems 


Lithium hydride was found to be considerably more energetic with 
fluorine than with oxygen. In faet, liquid fluorine, oxygen difluoride, 
nitrogen tetroxide, and chlorine trifluoride all yielded substantially 
higher performance with lithium hydride than did any of the all-oxygen 
oxidizers, Combustion temperatures ranged between 4000°K and 5000°K. 
The spread in performance was relatively small for liquid oxygen, hydro- 
gen peroxide, nitronium perchlorate, ammonium perchlorate, and nitrogen 
tetroxide, varying from 263 to 250 seconds at pe/pe = 68, 

fhe maximum specific impulse of the LiH-Fo and LiH-NF3 systems 
occurs near the point of stoichiometric combustion to Lif and HF. 
However, peak performance for the interhalide C1F3 ie defined by 
stoichiometric combustion to HF, LiF, and Licl. The affinity of chlorine 
for lithium appears stronger than for any of the other Jight metals. 
Lithium chloride is the principal chlorine-containing specie in the 
products of combustion. Both LiF and Hf are highly exothermic energy 
sources as well as fairly low molecular weight working fluids. The 
dimer and trimer of lithium fluoride become significant exhaust species 
at low pressures and temperatures. Although their concentrations are 
negligible in the combustion chamber, these polymers may constitute 
more than 50% of the lithium fluoride at the sub-atmospheric exhaust 
conditions. The dimer of lithium chloride was found to be considerably 
less significant. 

For the oxygen-containing oxidizers, maximum performance usually 


occurs nesr the point of stoichiometric combustion to Lig0 and Ho. 
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One exception is hydrogen peroxide, for which the peak is very close to 
stoichiometric combustion to LiOH and Hp0. In all of these systems, 
lithium hydroxide is the principal lithium-containing gas specie. 
Liquid liti‘um oxide and water are the other major products of com- 
bustion. Upon reacting fuel-rich, considerable quantities of hydro- 
gen gas are also evolved. This accounts for the gradual slope of the 
performance curves over the fuel-rich composition range. Once again 
He05 and NH,C10,, with relatively exothermic heats-of-formation, yield 
almost as much performance as their more oxygen-rich counterparts, Op 
and NOoC10,. Again this points out the necessity for a low molecular 
weight working fluid to accompany the formation of an exothermic con- 
densed phase, such as liquid lithium oxide. Although lithiw hydride 
is competitive with pentaborane, aluminum hydride and hydrazine when 
reacted with Fo, NF3, and ClF3, it is quite inferior with the oxygen- 
containing oxidizers. 
10. Magnesium Hydride Systems 

Moenesium hydride was found to be somewhat like lithium hydride 


in that superior performance igs obtained by oxidizing with fluorine, 





However, this effect is not as pronounced as with lithiwm hydride. 
Although liquid fluorine yielded the greatest performance with mag- 
nesium hydride, this was not exceptionally energetic. At pe/pe = 68, 


& maximum specific impulse of 329 seconds was obtained. This occurred | 





near the point of stoichiometric combustion to MgFo and HF. The same 


held for the NF3-MgHo system. With chlorine trifluoride, maximum per- 


| 


formance was achieved by reacting to HF, HCl, and MgFp. 
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Hydroger. peroxide yielded the best performance of the oxygen- 
containing oxidizers. Closely grouped behind were Oo, NH,C10,, NoO, 
and NOsC10) respectively. Again, the lack of low molecular weight 
working fluid limited the usefulness of Oo and NOpC10,. Condensed 
magnesium oxide, hydrogen gas, water, and magnesium hydroxide consti- 


tuted the principal combustion products. The first tyo are particularly 





predominant over the entire range of combustion stoichiometry. As in 


ee 


the case of lithium and beryllium, the predominant gas metal-oxygen 


specie was the metallic hydroxide. One outstanding characteristic of 


the magnesium-cxygen systems is the very gradual slope of the perforn- 
ance curves over relatively wide ranges of composition. 
ik. Generel Conclusions 

Of the non-metallic fuels, liquid hydroger was found to yield the 
greatest performance with each of the ten oxidizers. This was to be 
expected because of the energetic reactions of hydrogen with both oxygen 
and fluorine. In most cases, the remaining order of performance was 
hydrazine, UDMH, polyethylene, and dicyanoacetylene. From this, it may 
be Séeciwaad that the N-H bond is more desirable than the C-H bond iu a 


rocket fuel. Most important of all is that the hydrogen content of the 





fuel be high. Even the presence of a highly endothermic reactant, such 


as dicyanoacetylene, cannot overcome the value of hydrogen as in ingre- 
dient element. Of the metallic fuels, beryllium hydride was clearly 

superior, both with fluorine and oxygen as oxidizers. Pentaborane and 
aluminum hydride were found to be highly energetic with both oxygen and 


fluorine. Lithium appears to be a desirable fuel ingredient if oxidized 
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with fluorine. However, neither lithivm nor magnesium yields high 
performance if combined with -xvgen. 

Both fluorine and oxygen difluoride were found to be highly 
energetic oxidizers. Thus, it may be concluded that fluorine or 
mixtures of fluorine and oxygen yield greater performance than 
oxygen alone. An exception occurs when the carbon-to-hydrogen ratio 
of the ingredients is large. With both metallic and non-metallic 
fuels, the performance of oxygen difluoride exceeded that of nitrogen 
core iuorias which, in turn, surpassed that of chlorine trifluoride. 
These data indicate that O-F bonding, N-F bonding, and Cl-F bonding 
are preferred in that order. 

The relative ranking of the propellant systems at an exhaust 
pressure of one atmosphere was not always the same at the lower exit 
pressures. For example, although BeHo was the most energetic fuel in 
six of the ten cases at pe/pe = 68, it excelled in all ten cases at 
pe/pe = 5000. In general, expansion to sub-atmospheric pressures 
yielded relatively more performance for the metallia fuela. Any 
similer changes in the ranking of the oxidizers was not as sharply 


pronounced. 
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APPENDIX A 


SAMPLE CALCULATION USING THE METHOD OF EQUILIBRIUM CONSTANTS 


1. Statement of the Problem 


Moles/100 gm 
ight Prope! 
2 tole ty. 


| 85277 





0.8528 








| 
(Clip), 14.026 | 0 | 0.78 


Chamber Pressure = 1000 FPSIA 


15 | 1.0694 


Exhaust Presaure = 14.7 PSIA 


Density ef Propellant = 1.7715 @/cc 
Heat-of-Formation of Propellant = -36.633 Keal/100 gm 
Concentration of Ingredients per 100 gm propellant: 


C 


1.0694 gneatoms 


i] 


H = 4.2496 gm-atons 

O = 2.1107 gm-atans 

N = 0,5277 gm-atons 

Cl = 0.5277 gmeatons 

Al = 0.8528 gn-atons 
2. Mase Balances and Chepicel Equilibrium 

In solving this problem on a high-speed digital computer, all products 

of combustion which could exist would be considered until proven to exist 
in negligible quantities. However, thia rigorous epproach is not feasible 
when using a desk calculator. Therefore, the following species are selected 


for consideration in this sample calculation: H, Cl, A1(g)s CO, COo. Hp0, Bo, 
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fel, Np, Alo, 4100, Alci, ALCL a, 41203( 5 )s and 41203 (13° Since, in ell 


| 

| 

| 

probability, aluminum oxide will exist in only one condensed phase in the i 
chember and one condensed phase at the nozzle exit, this propellant system | 
] 


can be solved by 14 simultaneous eaustions. 


The meas bélances are: 


ZH = 42496 = ny + 2m, + 20,9 + BC] 
EN = 0.5277 = 2m, 





21 = 0.5277 = ng, * myc) * 30,10), * Barcr 
- a 
% a= 0.8526 =n4y/ ) * Mano * 204130 * March *Pa101,* 2241203, 2 * 2141003( 9) 


TO = 2.1107 = Ogg > 2N¢0, + 2HO0 + PAL0* PALZO* 39a1 2094) * FBAlgOg, g ) 








The equations of chemical equilibria may be written as: 


Be CO * H20 4g  605 + H, An=0 


= Ro02*FHD = Peo, * "He i *He 
7 Poo*Puso0 = "co * [HAO 


Speen eee Ee E 
a ee, 


be 1/2 HH An = 1/2 
ne me Ea 
Pig Mig ® Logs 
O. HOLE? Cl * 1/28 An = 1/2 
aes i/2 


_ Boastig’/? ney. 2 (2 


*P3 Puc. DBD) ap 
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| 
| 


ad. 1/2 41,0 


(a) * 2/2 COMA 4) # 1/2 Op An=1 


1 
a Palo + Pc ae DALOsNgog!/? =| 
*%, Poot mogt2 Lng 


6. 2Al0 + CO wg. 4120 2) * 005 An =i 


| Hp, = attr Thee = BAZ", 


5 
Pao *Foo = Fanos*®co) =P 


£5 Alo + CoO E> Aly.) # CO, An=@ 


PaisPco, = Alt*C02 
Pal0ePco Ra1o-2co 


Ge Al * 301 ao A1Gig An = +3 


Pall AAlCls ey 





“a Par Poi? Baa sto} Pr 


hs AiCl og. Al * CL An=1 


x, = AAG = “ALO. HE 


Paica Maici |= AT 


3+ Ghamber Composition iteration 


At an assumed copbustion temperature, the chamber composition is salou- 


iated. Although general methods of iteration utilising general equations of 
equiliorium and matrix techniques are preferred for use on computers, alge- 
bracio routines are faster for desk calculator solutions. Therefore, this 1i- 
lustrative exemple will be solved using trial-and-error algebracic iteration. 
An initial canbustion temperature of 3500°K ia assumed. At 3500°K, the 


equilibrium constants for the previously chowen equilibrium reactions ere: 


| | 
| 
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K, = Ky = 0.1248 

KE, = Ky = 0.5888 

Ky = Kip, = 9-3937 
K, = Eno = O.0746 
k,, = Euro = 0.1722 
Kg = Ky, = 0.01564 
K > uc, = 1.177 


Ky = Bassi = 06.134 
An initial estimate of ny, the total moles of gaseous products is 
np = 1/20LCl * FH + 220 + Fal * EN) 
Dp = 4.1483 moles 
Also, the initial estimate of 7H, is 
PH, = 0.45 (ZH = EC1l * 3Za1 - £0 * ZC) 
Dy, = 2.3576 
Thea» 
Bas hati? Boles = 0.061004 
Se 1000 RS XT ae? PSE 


1/2 1/2 
Pt) =o.zy70 {RLY Oxy, 
B 4 |? 4.0486 


Evaluating H, 


Dy = Kp (ag. 9(28)? = (0.5883) (1.5354) (0.2470) = 0.2233 





Evalueting HCl and Cl, 





€2 = 0,5277 = Bg) * Day * wei, * Dyion 
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ly ER 
eal a i ne ae) : 


i al Ai ee a an all Cot ao 





| 


Initially, assume that all chlorine exists as HOl or Cl. 





1/2 1/2 z 
Tapa (neg) =2) a 212 a5h ies 
1¢} i (= 0.3937 (4.0486) 15.7892 
« SOL . 25277 . 
ae ‘ ) 16.7692 = 90324 
Roy 
Byol = 6.4963 


Evaluating H,0, 

ZH = 4e2h96 = ny + 2ny, + 2ny 9 + Dy) 

10 21/2 (ZH - any, = Dy - Byg) 0.5926 

However, it would be physically impossible to have a negative quantity 
of water present. This meana that the original estimates of nx and My 
especially the latter, were too large. Assume a 5% smaller np and a 20% 
ameller Dy, and re~calculate, 
Let ap = (401483) (0.95) = 369409 moles 

om = (2.9576) (0.80) = 1,8861 moles 


Dp 309409 
i am = 0.057954 


Bg]? o.24074 Pap? 358 
Py| ep 
Re-evaluating H, 


2 
Ry = Ke (ng) = (0.5888) (1.373h) (0.24074) = 0.1947 





Re-evaluating HCl and a 


1/2 . 
DEOL = "Ho EW Pada (4.1538) = 14.4903 
7 OL ae 


22h. s 25277 wo, 
Rel (=: } 15.4903 0.0341 
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Din = 004936 
Re-eveluating H,0, 
0 =1/2 (2H ~ 2m = Dy = Dyg,) = -0-1055 
Again, this is a physica] impossibility. This time oo will be assumed 10% 


smaller than before while np remains 3.9409 moles. 





"H = (1.8861) (0.90) > 1.6975 ny? = 1.3029 
Re-eValuating H, 


| 1 
ty = Ky (ny, (38) 42 (0.5888) (1.3029) (0.24074) = 0.1849 


Re-evaluating HCl and Cl, 


Myo. 222029 
, Noi 0.3937 








(401538) = 13.7464 


Sd 35 


sary 


nei" SEER, = 0.0059 ty = O49 
Re-evaluating #50, 
"Ho = 1/2 ( ZH - ny, - Rp - Dysp) = ¢..9889 
Since water ia now present in realistic quantities, the iteration my he 


continued to other species. 


Evaluating CO and C05, 


2088: 
Pcp = Ky PHO - 0.1248 anata 0.0065359 
Bag "HL 1.6975 


¢ 
TaaepH chagg, 7 070065359 


2 
"00, * 0.00694 and nog = 1.0625 
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Bvaluating 410, 


ayo = Ki er pee] ‘ (.047h6) (6057954) = 40340 
ALO le EB 0.0069 . . 203403 


Evaluating 41,0, 
Barco = (Bano)? (Ks) “BO . PE = (03409)? (1722) (153610) (274255) 
"0p OF = 0.5269 
However, this is another physical impossibility since 2n4,_ 9 & Al. 
< 


The concentrations of Al0 and Al,0 are greatly affected by the ratio Peo, 
: Bg0, 
Therefore, this ratio will be introduced as a third unknown variable to 


facilitate rapid eonvergance of the aluminum balance. Next, sssume 


Xo = 756 fey 66 


Re~evaluating Al0, 


1/2 1 
Baio = K,|-22 BE] = (04746) (8.66) (057954) » 6.02362 
"0, Py 


Re-sevaluating 41,0, 


= 2 Xo , Pr 2708 a 
341,,0 K, (2,30) we oe (3722) (02382) (75-0) (17.255) 


Then, other specicsz are: 


Byy = (Myx) (K) Po = (.02982) (01564) (75.0) = 0.02794 
2 


Dale (02794) (0358) _ 
Cs ae ae = Ollhs (17.255) = 9.1507 


By1c13 7 (Mgq) (0¢3)9 (7) [zy 


Bgior, * (402794) (3658 & 1072)? (1.77) (417255 X 102)? = 0.0075 
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Re-avaluating Hp0, CO, and CO, 


a 
2 =x, 22 = (0,2048) (75) = 9.360 
"20 CO, 


If My, = 1.6975, then Bi,0 = 0.1814. 


Boo, 10691 = Ngo 
Bog = 160553 
Re-evelueting HCl and Cl, 


typ, = ZH = 2Dy,, - By - 200 = 0.3069 


=Cl-n - 3n = 9035375 

ie A1C1 ALCIig = I=0O, 

oa lg = 53 0.0228 
PHC] + t 
M1 


The other species are: 


Raye (02794) (.0228) 
Palcl * — = = STE (17255) = 0.0960 





n = (02794) (2.28 X 3072) (1.177) (37255 x 102) = 06,0026 
‘ALCL, = he 9h) « ob f ( Ve 
Pal, ~ 1/2 (41 ~ Baacg) ~ Maro ~ *Mazz0 ~ Malcn m41C14) 
Checking by oxygen balance, = 0.2252 
LO = G9 * 28co, * Bua * Maro * M120 * 3%a2203/ 5) 
£0 = 2.1107 = 2.0907 
Checking by total moles balance, 
Dp = Dy, * 5 * BHOL * Aco * Sco, * MH * ACL * Baro * "algo * Mad 


* Dalch, * DNS * M4101 
165 
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ap = 329409 S 4.0029 
Checking by chlorine balance, 

2 C1 = ng, * Oya * My3¢2 * MarcL, 

2Gl = 0.5277 # 0.4517 

None of the preceding three material balances is completely satisfied. 
However, it is highly improbable that the correct combustion temperature was 
guessed. Therefore, an enthalpy balance will be made to determine the degree 
of error in the assumed combustion temperature. Of course, in a high speed 


digital computer, new values of ny, Ry ngof/tgg,, Would be assumed until the 





metsriel balances each converged. 






































b he Epthalpy Balance at 350 
K Specie 1598 
j ce 1.6975 | 25.702 G 
t BR 0.1849 | 15.906 4522089 
KE 
¢ HO 0.1814 36.933 #57 0798 
i co 4.0559 | 26.821 26.416 
H ca, 4.605 =94.052 
k Hol 0.3069 | 26.347 22.063 
| cr 0.0228 | 16.359 428.042 
4 Ald 0.0238 | 28.400 416.9 
i 
41,0 021264. | M4dje.080 =39 0h; 
AL (g) 0.0279 115.946 *77 oh 
alcl 0.0960 | 28.414 50.0 
ac, 0.020 | 62.700 | ~137.7 
A1-G5,)| 0.2252 | 151-10 | ~373005 
N 0.2639 | 26,605 | a 
: | 136,618 | | Sib ear 
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erie Seo oS ed ~ : — _ Te Ss wate ae vas et < Se, 


i Sey ot, rat tani ee = 





“CEnjHy) = AH, = (En, ABe - (in; Ar,» 
=136.8:.8 £ ~127.681 + 36,633 
#136.818 # -91.048 
Thus, the enthalpy balance is not satisfied. Since the assumed temperature 
wes too high, & combustion temperature of 3100°K will be assumed, 
Se Ghamber Composition Iteration at 3100°K 


Ky = Ky = 0.1349 


% 


ae” 








«| Xp = Ky = 0.2221 

K, = Kina = 0.3453 

’ Ky = Kyo = 30618 X 1079 
Ks = Kyroo = 1-2u0 

z Ky = Kyy = 1-44) X 107% 


Ky = Eyo1, = 336.2 
Ky ® Kyoy = 1-hal X uo"? 


Assume ny = 3.8600 moles Then, nq/pp = 0.056764 
iy, = 1.8410 moles pp/n; = 1726166 
Dgo/tge, = 220 (np/pg)¥4? = 0.29825 


(ppdog)*4? © 4.1972 
The chember camposition then becomes: 
Dy, = 0.2639 
Xoo = 1.0646 
+ (ne 42) [ EP? 

ny * (ko) (RED p<) E = (.2121) (1.3568) (23825) = 0,0686 
we = x, 22 = (3349) (220) = 29.678 

™H20 = C02 
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nar cap enone mean ee ee ya tet Et neti = gee eS ane pH 


ees 


ee ne nt me 


130 = 0.0620 


Dp, = LCl = 2ny, - 2nHo ~ By = 0.3750 


1/2 12 (2.63568) 
Pee) a oe at Pp tithe Gle é | = 99.) 
aay K, BE .1h53) (4.1972) = 392193 


Bg, = 0.009578 


i/2 
Dao sa (x,) BE = (2448524) (003618) (056764) = 0.0030 
No. Pp 


Byp0 BS) (taro)? 3 PE = (2.241) ( 003046)" (220) (2746366) 
2 Zovolks 
nay = CK) Gano) Bee = (.01441) (003046) (220) = 0.0097 


ie 006 ‘ 8 
Daici * aoa) z = {-cestse) (ee (1726166) © 0.1147 


Eason, ~ (Ep) (man) (Mon)? | EL = 0.0156 
3 ny 
PAL 203¢ 5) = 1/2 (Fal -Dy0 7 204150 - Day - 24101 * 41014) = 0.31603 
Checking oxygen balance, 
FO = 2.1107 = 2.1147 
Checking chlorine balance, 
Tol = 0.5277 = 0.5461 (Satisfactory checks} 


Checking total moles balance, 


Dp = 3.8771 4 3.8600 
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6. Enthaipy Balance at 3100°K 








Specie 1598 
1.8410] 22.210 | 40.889 

x 0.0686 | 13.919 | 0.955 #52089 
0 0.0620) 31.530 | 1955 576798 
ca 1.0646 | 23.246 | 24.748 =26.416 
co, | @,0048 | 38.021 | 0.183, ~94,.052 
Ho1 0.375 22.756 8.534 =22,063 . 
Ci 0.0096 | 14.529 | 0.139 #28942 
Al0 0,003 | 24.726 | 0.074 #169 
41,0 0.0446 | 38.51 | 1.718 =39 01: 
AL | 0.0097 | 13.957 | 0.235 +77 oh 
Alcl 0.1147 | 24.838 | 2.849 -50.0 
aici, 0.0156 | 54.78 | 0.855 “13707 
A1o0q/,) | 93103 | 115678 | 35-927 -373009 
N, 42639 |23.046 | 6,082 0 

| 308771 125-03 


~(in,H)), = AH, = (E03 AB 298) ~ (ny Ake og )n 


#125043 £ 161.197 + 36.639 
7125.043 S -124.564 


Thus, the enthalpy balance is satisfied, 
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+3053 
30583 
28,122 
~0.451 
8.274, 
+0.278 
+9.051 
“1.757 
40.751 
5 +735 
2.148 
115-780 
Q 


-161.197 


— a ecencesenemtetsseanstssanaanineeitaaetne eat e228 : 
i 
i 


| 
| 
i 
! 
\ 
fe 


et ee ae re ce ter ee 





7e Calculation of Chamber Entropy end Heat Capacity 







Specie “4 ln = | = ln z 


=0.354 

















1.8410 | 48.761 | 89-769 8.912 | 16.407 
































H 0.0686 | 39.025 | 2.677 ~0.072 | 4.968 | 0.342 
10 bites 68.856! 4.269 20,066 =| 13.374| 0.829 
ca: 1.0646 65.662 69.904 1| -0.355 | 8918 | 9.486 
ca, 0.0048 | 80.328 | 0.386 | .0.007 /14.962| 0.072 
M1 | 0.3750 | 62.656] 23.4h6 | | 0.226 { 8.929 | 3.348 
Cl | 0.0096 | 51.788 | 0.497 | 0.002 | =6.215 | -a.012 | 5.030 0.048 
Alo 0.0030 | 72.164] 0.216 0,001 | ~6.908 | «0.007 9.165 | 0.027 
AL@ | 0.0446} 96010 | 40286 | 0,012 | 4.423 | -0.053 113-90 | 0.620 





Al e) 0.0097 | 51.005 | 0.495 0.002 | -6.215 
acl GeLIA7 | 7h.99h | 8.602 | 0,030 | -3.507 
4101, 0.0156 | 119.40} 1.863 | 0,004 | 5.521 





0.012 | 4.970 

9-086 

19.84. 

4103/9) oad ace 28,616 | ~---= | --- wae | --n--- 38015 

m | 242629 | 64.046 162902 | 0,068 | «2.688 8.868 
| 


| 3.8771 | 251.978 


9,048 
1,042 
0.310 





0,105 


11.838 


8, z 22385 - et = ln = = De.Reln P, 


S, = EnjS; - Reng (eo) in?i + In Py) 
Dp OT 


S, = 251.978 = 7.704 (-1.47h * 42220) 
S_ = 251.978 - 21.155 
8, * 230.623 
200 gn propellant _ 
Vean Molecular Weight = Eyal moles of gas ~ 79°79 
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8. Exhaust Composition Iteration 


At an assumed exhaust temperature, the exhaust composition is calculated. 
Aesume the exhaust temperature is 2150°K, Since this is below the melting 
point of 41203. the species Al0, Al20, Alig). and 1205/9) will be neglected 
in the exhaust composition. Also, the concentrations of Cl end AlCl will be 
Mall in comparison to those of HCl and AlC1,. Therefore, the following 


species will be considered: CO, COn, HoG, H5, Hol, Hy No» AlCl, and A1202( 93° 


The nsv¥ mass balances became: 
2H = 4.2496 = any, * 2ny0 * Nap) * Ay 
EC = 1.0694 = nog + 2605 
ZN = 0.5277 = 2ny,, 
ZCl = 0.5277 = ny + SaiCly 
Zal = 0.8528 = 2241003," M4161, 


LO = 2.1107 = noo * 22¢05 + °H,0 + FAL 2054) 


Necesaery equilibria are: 





@e CO + Hp0 ~—g—— C05 + H, An =0 
= Bo0p*Pig 200, 
* Poo*Pi,o = "0° "20 
be 1/2 Ho =q— H an = 1/2 


he ae 7 St By 
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Assume np = 3.6400 moles 


K, = Ky = 0.19786 
y= K, = 0.004366 
Then, np/Pp = 3.6400 


= 2,0950 moles (np/ pp) /2 = 1.9079 


"Hp 
Qco = 1.0693 
B05 = 0.0001 


Ty, = 0.2639 


1/2 
n, = (Kp) (ng?) [Z| = 0,012 


711005) = 1/3 ( £0 = Ngg = 2M¢0, ~ nyo) * 0.3467 


"AlCl, = Dal - 20412055) = 0.1594 


Checking chloriuze balance 


TC = 0.5277 = 0.5235 (Satisfactory Check) 


Checking total moles 


Rg = 306462 = 320400 
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9. Entropy Balance at 2150°K 


Constituent ny 


2.0950 






45213 | 94.721 































H 0.0121 | 37.206 | 0.450 

H,0 0.0011 64.2121 0.072 | 0,000 | ~---- 
co 1.0693 | 62.434 66,671 — 

co, 0.0001 740939 0,007 D | meme 
Hl | 0.0453 | 590446 | 2.693 

AlCl, 91594 | 112.11 | 17.870 

Al503¢g) | 063467 | 674763 | 230493 | ---n- | ------ 
Ny 0.2639 | 60.836 | 16,055 


3.6462 222.121 


Sy = Inj8j - Rony (ri In Z + in P,) 
S, = 222.121 + (1.074) (7.245) 
S, = 222.121 + 7.781 


S, = 229.902 


S, = 229.902 = 230.623 = S, (Satisfactory check) 
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10. Calculation of Exhaust Enthal Heat Capacity, and Heat-of-Formation at 
2150°K 








m4 


2.0950 | 13.596 | 29cll2 

















xy 0.0121 | 9.200 0.111 $52,089 +0630 
#0 @.0012 | 19.222 | 0.021 572798 6,064 
rote) 1.0693 | 14.866 | 15.894. -26,416 -28,247 
co, 0.0001 | 24.027 | 0.002 94.052 | 0,009 
Hol 0.0453 | 12.710 | 0.576 =222063 | 0.999 
s1C1, Qe1594 | 35093 | 5-727 -137.7 “21.949 
$1203(9)| 963467 | 55¢29 | 19.169 “400.4 -138.819 
Ny 202639 | 14-711 | 3,882 Q 0 

306462 7heh96 | 7189457 


i00 t oT 
Mean Molecular Weight = 3°¢)20 moles of gas ~ 27 h3 


ll. thalpy Ba. 
Since absolute quantities of enthalpy cannot be assigned, the cverali 
enthelpy balance must be celculated for the complete path followed in the 


rocket motors 


Ty 
Te 
Reference Temp = 298K 
Chamber Exhaust 


t v 
nH i098 + ARDS + En jug + (ZnsAze,), - (Eny Bie), =e 9 
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See A 


Sid le cae ca 


T, = 3100 


AH =? 
AH = 125.043 Te = 2150 
AH = 7heh96 
| oe 
none be) 


AH, = 28,260 


125,03 - AEC - 7heh96 * 28.260 = @ 
Ave = 78.807 Keal/10¢ ¢m 
8 


ue 12. Sypmery of Chamber Gonditions 
Pressure = 1009 PSTA 


pe 


carey eeu wt 


Temperature = 3100°K 


rc 


¢ 


Enthalpy = 125.043 Keal/100 gm propellant 
Entropy = 230.823 calories /1L00 gm = °E 
Heat Capacity = 46.756 calories /100 gm - %K 
Moles of Gas = 308771 moles/10C gm 
Total Moles = 4.187) moles/100 gn 
Mean Molecular Weight = 25.79 ga/mole 
Gas Composition, moles/100 gn propellant 
H, = 1.8410 
H = 0.0686 
HpO: = 0.0620 
CO = 1.0646 
CO, = 0.0045 
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13. 


HO1 = 0.3750 
Cl = 0.0096 
A410 = 0.0030 
41,0 = 0.0446 
Aly g) = 0.0097 
AiG] = 0.1147 
AlCl, = 0.0156 


A1,03/,, = 0.3103 


M2) 
Hp = 022639 
Summary of Bxheust Condi tiong 
Fressure = 14.7 BSIA 
Temperature = 2156°%E 
Bnthelpy = 74.496 Keal/100 gm propellant 
Entropy = 229.902 calories/100 gui - °E 
Heat Capacity = 4).357 calories/100 gm - °K 
Moles of Gas * 3.6462 moles/100 gn 
Total moles = 3.9929 moles/100 gm 
Mean Molecular Weight = 27.43 gn/mole 
Gas Composition, moles/10@ gm propellant 
H, = 2.0950 
H # 90,0121 
Hp0 = 0.0011 
CO = 1.0693 
CO, = 0,009L 


HCL = 0.0453 
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14. 


ae 


be 


d. 


Ge 


fe 


pews Wie aes aE MRE eee 


AlCL, = 0.1594 
A1,03¢4) = 963467 
¥, = 0.2639 

ation cket Perf ce Paremete 
Isp = 9.330 \/H, = Hy = 9.330 \/788.07 
Isp = (9.330) (28-071) 
Isp = 261.9 seconds 


Ivep = Petey = 1°75 & a ch PP SBbror tas (261-9 sec 


Ivep = 16,75 lb-sec/in.2 
, = Ispeg = (261.9 seconds) (32.2 ft/sec) 


u, = 6433, ft/sec 


r 
=: 
ka=2 =e * = = ane eke = 0.0868 





2. 
Cy = dk : a 1 - Te. 
keel kel z, 


Cy = \/ c23enu) (0.953) "97 (o.gi2) = 7.513) (0.373) 


Gy = 1.652 
4guk [2 a Py 1/ 
a [FI 


= 12096 22639 (payed _ 
Lor (0.953) (68) (665), (357) (46.7) 


= 11.09 


wl? I 


: 











bee Tp 14652 31 


c* = 5100 ft/sec 


o 60756 
is op =sk=e= ae = g g 
C Tp-mp-R 466756 = (328771) (1.987) 
chember 
k = 1.199 


ge Pp Ae a L o 
f= = Hels. = 
m a ) shifting TeI5Cu7 (326462) 261.9 


Pete = 25.9 


np T, k+l 


2 
jo gefce* = =e 1/2 2. 2(k-1) 


2 
1/2. a8 
2. 
g/e* = 1.074} 308771) ( sri5] | Eas5| 


@/c* = (1,074) (0.01) (0.587) = 0.00631 


ct = = 5100 ft/sec 


0.00631 


Taree _ (26 
Ke Cy = c* = 5100 


Cp = 1.655 


ls 4, 6 Be Bebe _ 2-2 1000, 
ie SG em * Hoe * Tey 9) 


Se #1112 
a* 


Me Ideal boost velocity data are evaluated sa follows: 


178 





Vy = Ispegeln | wk |- Ispegeln l1+* i ae | 
Py 


We A, 


x 


| 

| 

¢ =mass of missile hardware 
V.. = volume of propellant 
f = density of propellant 
Wy = total mass of miasile 
*) 


Wi, = masa of propellant 









on * eee ; 

b ‘m rw, Ve-¥> y 

é 0.95 20 3.00 25, 300 

: . 0.90 10 2.30 19,400 

i 0.8 6.67 1.90 16,020 

4 | 0.80 5,0 1.61 13,580 

is 

f 0.70 3633 1.20 10,210 

t 0.50 2,0 0.69 5820 

i 

\* 

I: be os Ta 1+ cr in [1 ea i Vp 

u 2! hela «45 dpa 1.690 Th.250 
: 6a 1.842 2,842 1.040 8,770 
é 

t 100 1,005 2,005 0.693 5.840 
E 

t 

€ 

‘ 

> 
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APPENDIX B 


SAMPLE CALCULATION USING THE METHOD OF MINIMIZATION OF FREE ENERGY 


1. Statement of the Problem 






AHpa98 
Keal/mole | @m/cc 


Moles/100 gm 
Weight % | Propellant 






Ingredient 





1.8162 


ave 


$12.05 1.004 


92 = Soh 1.43 


58.207 






NoH), 
N20), 41.793 0.4542 
Chamber Pressure = 500 PSIA 
Exbaust Pressure = 14.7PSIA 
Density of Propellant = 1.148 g/cc 
Heat-of-Formation of Propellant = +19. Keal/100 gm 
Concentraticn of Ingredients per 100 gm Propellant: 
H = 7.2648 gm~atoma 


N 


4.5408 gm-atons 


0 = 1.8168 gm~atoms 


2. Mass Balance 





In solving this problem on a high speed digital computer, all pro~ 
ducts of combustion which could exist would be considered until proven 
to be present in negligible quantities. However, such a rigorous ap- 
proach is beyond the scope and feasibility of this sample calculation. 
In this example, the purpose is solely to demonstrate the calculation 
method. Therefore, the products-of-combustion will be limited to the 


following species: Hy, Ho0, No, Hy N, OH, and 0. The mass balances are: 
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a a a Gee 
a ee PS 





ZH = 7.2648 = any, + 2nH.0 * Noy * By 
IN 4.5408 = ang, + Dy 
LO = 1.8168 = 70 * No, * A, 


3- Initial Conditions 
Initially, in a general solution, it is assumed that the only products of 


combustion are the elemental gas species. At 298°K, the properties of N, h, 


and 0 are: 


Specie 


H 
a 


From these data, an enthalpy balance is made with AT the variables 
~ E00) products. * C2 ny ABs 58 products ~ (any ABs 96 reastants 


( ies Ak, 298) 





#19,440 calories/100 gn = 998,500 + 68.1 AT 
AT = -14,400°K 
Thia, of course, is a hypothetical temperature corresponding to tae existence 


of species N, H, and O at the constant enthalpy of the propellant. Whenever 


reactants 


Cc PAY | 
4.97 cal/deg-mole #113.0 Keal /mole 
4.968 cal/deg-mole 452,09 Keal/mole 
52 cal/deg-mole +59.55 Keal/mple 


@) 


=(5 ny AH 98) 


hal nC A 
products (z J Pjo98™ “products 
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the hypothetice? temperature of tha mixture falls below 298K, thermodynamic 
data at 298°K are used for the various product species until the trial 


temperature raises above that value. 


Each of the various products of combustion will now be examined, one by 
one, to determine which specie will cause the greatest reduction in total 
free energy of the system if introduced as a constituent. Speciea are exemined 
and intreduced one at a time. 


4. Formation of No 





Nitrogea gas is introduced into the system by formation from nitrogen atoms: 


2N—="N, 





2 Any, ze- Any 


The formation of No wiil cause both a temperature and concentration change. At 


-AT = Eno - 2 Hy Any 
ani", 


‘ 
or -AT=k Any, 
e 


where k = -2 
"No By 


| 

| 
a given temperature, this relationship is | 
| 
| 

Xa,C 
"ORs; | 


When the temperature of the mixture is a hypothetical value below 298°K, x iz 





considered zero and the above relationship loses significance. However, when 


the temperature rises above 298°K, the above expression is used to relate 


| 

AT end An. 7 
1 

| 
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~ ] 
4t 298°K, the ideal entropy and free energy of N and Ny ares 
< 
Sy * 36.6] cal/deg-mole 
{ 


6 
Sy, ® 45.77 cal/deg-mole 
Fy = HATS" = # 113.0 - 10,90 = 4102.1 Keal/mole 
Fy, = HTS” = 0 - 13.65 = -13.65 Keal/mole 


For the general case, the approximate quantity of No which should be introduced 
in order to minimize the free energy of the system is defined by the equation | 
; 3 6 2. Ate A -2 

ArSa Any, +B (Any,)° + 6'( Aay,)? 


wuere A’ =cel/mole=k'ay (-Sy)q + ' oy, (Sh Mp # i! ny SB)p + 'ng( 85), = 20E ly 
+ yp + 2RT (1 # In ny). + RE (1 * In ny) + k'R (ny ln ny) + ie" RC ny, in ny) 


+ &'R (ny ln ny) + k'R (np ln ng) = RE (1n np #1) - RE In pp + k Reng Is pray 


B' = oel/(mole)* = URE + BT. ~2k'(-Sy)y * x’ (28 Ye - 2k'R (1 * An ny) tek R 
"NBN 7 
(1 + In ny) = k'R (1m np #1) = k'R In pp * BP (np, = np) + mg an OE 
a 
(np we 7 Bl? 
c= cal/(mole)2 = uk B *k'R + k'R (Dg ox - Xp) + ee in ot 
(type ~ Be) 
Then, An, = =p. + Vip)? ~ a'o! 
2 30° 
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However, k' will be zero at first corresponiing to the low hypothetical 
caaparature: Under this condition, the general cubic equation reduces 
to a quedratic, Furthermors, Dy, is initially zero, Same finite value 
of Np must be asgigned, however, in order to avoid infinite contri- 
butions to the scefficients A’, BY, C', For the general cease, 2/3 of 
the maximum aliowable concentration is chosen os the initial quantity 
to be introduced. This velue corresponds to the point of tangency at 
which, on the average, the y = ex? = bx curve epproximetes the 

y = x ls x curve. Thus, the scefficients are evalucted es: 

AY = -2(102,100) + (13.650) = 2(2.987) (298) (1 # In 4.5408) + (1.987) 
(298) (2 # In 1.5136) = (1.987) (298) (in 13.6224 #1) -(1.987) (298) 
(ln 34.0) 

A! = =204,200 ~ 13,650 = 2976 + 837 -2139 -2089 

A’ & «224,217 


BY = 4(1.987) (298) + (3.987) (298) + (2.987) (298) (-0.5) 
b.5408 1.5136 


B! & 616.7 
ofa 
In this sasz, 
AF Sa (any) * BY (Amy)? 
A¥ © -224,217 (Amy) + 616.7 (Amy,9? 
Differentiating and setting the derivative equal to zero, 
Gs 224,217 + 125304 Ary, 
Ay, = 181.8 
However, according to the masa balance limitation, a maximum of 2.270h 
moles of N- can be formed. 
Therefore, since 161.8 exceeds 2.2704, all nitrogen should exist 


16h, 











as Noe If all N were transformed into No, the total change in ideal free 
shergy of the system at 298°K sould be: 


AF = #494.0 Keal 


5+ Formation of H, 


Consider the introduction of hydrogen gas into the system by formation 
from hydrogen atoms: 
2H Hp 
2 Any, = ~ Any 


The relationship between temperature and concentration changea is 


«AT = 





Hn - 2Hy 
a A s Any, 
3 Ps 
&t 298°K, the ideal entropy and free energy of H ani H, are: 
Ss = 27.39 cal/deg-mole 


oO 
S = 31.21 cal/deg-mole 


re =H -TS” = 452,09 - 8.17 = +43.92 Keal/mole 


°o 
F. =H-TS = 0 = 9.30 = -9.30 Keal/mole 


He 


Proceeding ss before, 
a’ = -2Fa ey + a ~ 2RP (1 # In ng) # RE (1 * In ng) - RP (1n ny + 2) 


- RE ln pp 


A’ = ~103,783 
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a re 


| 
| 
I, 
ee 
\ 
I 
iE 
V, 


~ np 
p' + URD + RT + RT (Ot yax ap) ae 
my OR, aca Sa ara 
(nga ~ Or? 
B= 27h 
c'2@ 


Then, 548.8 Any, = 103,783 =@ 
Any = 189.2 


However, according to the mass balance limitations, a maximum of 3.6324 moles 
of H, can be formed. Therefore, #ince 189.1 exceeds 3.6324, all hydrogea 
should exiat as He in preference to existing as H. For this case, the change 
in ideal free energy of the system at 298°K would bes 

AF = 3.6324 (-9.30) = 762648 (43.92) 

A¥ = ~352,8 Keoal 
6. Formation of H;0 





Next, consider the formation of water from hydrogen and oxygen: 
2H + OL” Ho 
40 ee Ano “ 2Ang 
Tram the oxygen balence, a maximm of 1.6168 moles of water could be formed. 


From the hydrogen balance, a maximum of 3.6324 moles of water could be formed, 
Hence, the oxygen balance is the restricting one, The general relationahip 


between AT and Anu.o is 


~ At = tz0__“0 - °8 fo ~ Pha Atyog 


Znc 
ae 


186 


aot Aap oT To a 5 : 
a kt 











At 298°K, the ideal entropy and free energy of H, O, and Ho0 are: 


Sy = 27.39 cal/deg-mole 


Lo) 
oo 
it 


38.47 cal/deg-mole 


9 = 45.11 cal/deg-mcle 


we 


mo 
u 
ize 
$ 
a 
" 


443.92 Keal/mole 


ty 
° 

u 
ie) 

t 
3 
n 

© 

NN 


448.07 Keal/mole 


So 


oO BO Bae 6h Sol Bek te ae 
F.. =SH=TS = “71.21 Keal/mois 


H,0 


In this case, the respective constants to the cubic equation are: 


ad 
tt 


° ° 8 
CR dg - (Foe * (Fa odp > 2RE (1 + In ng) - RE (1 ln no) + RE 





(1 + In ny 09) - RE (in np +1) + RT in pp 


> 
j 


= 215,119 





: (aq - og) * Dpnex ln Z 
b p RE Lae mex © fax 
F a 0 "Hoa (nq - np)* 





B = Bhs? 
cl =a 
Then, 1689.4 Ang,o - 215,119 = a 
Anyoo = 127-3 
However, a maximum of ),8168 moles of water can exist. Since Ary o exceeds 
this value, all oxygen shouid scxist as water in preference to oxygen etoms. 


Furthermore, 3.6324 moles of hydrogen should also be combined. The ideal free 





energy change accanpanying this concentration change is 
AF = 1,8168( -71.21) = 1.8168( 948.07) - 3.6336(+43.92) 


AF = -37569 Keal 
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7e Introduction of Major Species 

In a general conputer solution, all other species, including OH, would 
be evaluated in the same manner. In fact, the existence of Ho, No, and HHO 
as the predaninant species could be reasoned from intuition in this rele 
tiv.iy simple example. Nevertheless, the preceding procedure was followed to 
demonstrate the fundemental approach. Since the greatest single reduction iam 
free sueigy occurs by transforming N to No, the first new trial composition 
will bes 

N, = 2.270 moles 


O = 1.8168 moles 


H = 7.2648 moles 









¥, 2.2704 | 0 

z 722648 +52 409 437705 4.968 36 el 

0 | 1.2168 | 159.55 | 100 5.2h “265 | 
#48525 61.4 


19.440 = 485,500 + 61.4 AT 
AT = =7580°K 
Since the introduction of Nz was not sufficient to raise the temperature 


of the mixture to a realistic value, water will also te introduced’ 








-105.0 


3.6336 | +52.09 *389,0 
+ B4.G 





2 


1Bek 
419,440 = #84,000 + 48.4 AT 
AT = -1330°K 
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2 270k 






0 1.8168 
H, 1.8168 


529040 
#19,440 = -105,000 + 42.9 AT 
AT =.2900°K 
This corresponds to T, = 3198°K, However, the heat capacity data used in 
obtaining this temperature were extrapolated over a wide range, An enthalpy 
balance is therefore needed to adjust to the corresponding chamber temperature. 
4t 2700°K, 
(z *j AH 1098 \seactante . 
» 12h.4y S 125.3 


and the enthalpy balance is essentially satisfied. Actually, by extrapolation 


a T 
(Z ay OPE y ercunte ie Az ny Joga products 


t = 2687°K. The free energy of the system at 2700°K is 


Specie n [ ¥ Sed In wd Bi ae Et 
} j oT | np 






| 348.4 | 1.000 | 


5920 


Fo injHy - T2ZnsS, + Rfnp bes in Zi] Ring 1a pp 





F = 420,300 ~ (2700) (348.4) + (1.987) (£700) (5.9040) (3.526 - 1.093) 
F = 420,300 ~ 940,080 * 77,065 


¥F = ~842,715 calories/100 gm propellant 
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8. Introduction of Minor Species 
At 2700°K, same dissoci-“ion will be present. Therefore, the species 
N, H, 0, and OH will be re-introduced. First. atomic hydrogen will be 
considered. Its concentration may be approximated fran the decamposition of 
hydrogen gas by the reaction 
Hou 2H 
2 Any, = Any 
Any = + 1/2Ang 
At constant enthalpy, the formation of E will cause both q temperature and a 


concentration change. The enthalpy balance becomes 


AH =Q= Any (Hy; - 1/2 Hy) + Enjep, At 


- At = hy - 1/2 By Ang 


Lajey 


At 2700°K, na ole02 = 9230 


59.0 Ay 
At = -0,928 Any 
Si, = 47.23 cal/deg-mole 


S_ = 38.34 cal/deg-mole 


m0 


fo) 


F =H-TS = 18.59 = (2700)(,94723) = -109.6 Keal/mole 


° 


F,th- TS’ = 64.02 ~ (2700)( 403834) = -39.4 Keel /mole 


Thus, at thia point, the species Ho, Nos HHO, and EK are present az products 


of combustion. The constenta for the cubic equation become 
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Als ka, (-Sho de + k'ny (Shp + k'ny,, (“Shy + king (Sino) ~1/2 
(Fi dp + (Fgdp © 1/2 Rt (1 4 In ng.) + RE (1 + In ny) + k'R (my, In Dgp) 
* k'R (ng In nq) + k'R (ny, In nya) + k'R (npg9 In nyo) ~ 1/2 RT 
(in np #1) #1/2 RE In pp * k'R ap in Ee 


BY = BT + i - = 1/2 (“Si Jap + (“Sg )y > 2/2 KR (1 + In ngp) + KIR 
hoy, 


nice tannery omer et mee 
A (one, ° ~ np)? 


areas 


LN EIEN ERI SIRES LEN IE ARIS TION 


* 1/2 k'R ln pp 
Cc! = #kIR ( )+ n 
a ELE Pypax * By) * Oppay 2 eS 





Mi ny)? 
Sama finite value of 5, mst be initially assigned. Although 2/3 
of the maximum allowable concentration would be initially assumed for 
the general case, the first estimate for this sample problem wes made 


from equilibrium data: | 


1/2 
mg * Hy (age?) 2 72 (0405669) (14340) 52308” 5 9,056 





Eveluating the constants, 

At # (0.928) (1.8168) (47.55) + (-0.928) (0.0326) (-38.34) * (-C.928) 
(2.2704) (62.84) # (0.928) (1.8263) (-66.92) - 1/2 (-10926) + (+3904) 
= 1/2 (2-987) (2700) (2.599) + (1.987) (2700) (-2.45) + (-0.928) 
(2.967) (1.081) + (-0.928) (1.987) (.109) + (-0.928) (1.987) (2.861) + 
(~0,928) (1.987) (1.082) = 2/2 (2.987) (2700) (2.775) + 1/2 (1.987) 
(2700) (3.526) + (-0.928) (1.987) (5.9040) (1.750) 
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ee re ne 


a’ © 80.17 + 1.12 + 132.40 + 112.83 + 54,800 - 39,400 = 4289.2 = 13514306 


= 1.99 0.2 - 3043 - 1.99 = 7443.8 * 9.458.3 ". 19,05 
a! = # 281.56 


} 8 2700 
pt = Ga EBD + MO ost ) = 1/2 (0.928) (47.55) # (9.928) 


(+38.34) = 1/2 (-0.928) (1.987) (1.599) # (0.928) (1.987) (2.45) 


(1326224) = 5.9040) + 13.6224 in oe 


* (1.987) (2700) ) 15.6224 + 1/2 (0.928 
(13.6224 = 5.9040) pines 


(3.967) (30526) = 2/2 (=0.928) (1.987) (2.775) 


B = 738.23 + 169477303 - 21.61 * 35258 * 1eh7 * A.52 =2 056 82.30 = 3245 
BY © 4170, 44304 


CG" © 20,2537 = 58.3525 * 0.0283 
O° = -58,578 
The change ir. concentration of hydrogen atoms is eveluated as: 

— 3 + Vos? - 31's! 

x 
ng = 27O+4beok #129,051,003,083.23, 
175 073A 

An, = -0.0114 moles 

Thus, the new concentration of H is 


By = 0,9316 - 0.0114 = 0.0202 moles 
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By continuing the iteration until all species have been re-introduced and the 


free energy of the system has been minimized, the gas composition 1s found te 


he< 
BN, = 2.270) moles 


"HO = 1.8112 


hy, = 1.8052 
He 

ny = 0.0266 

Qo = 0,0G00 
Ty = 0.0000 


9. Epthalpy B. e at 2650°K 


Specie Ry Awe 
2.2704 | 19.08 






43.319 





HO | 1.8112 45-914 
H, 1.8052 

H 6.0266 0.311 
OH 0.303 






2.0056 

5.9189 122.419 
“Easy, = (Za; ee Po ~ (Eny Abe nope 
122.419 = -103-246 - 19044 

-122.419 -122.686 


Thus, the enthalpy balance is satisfied at 2650°K, 
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-104.684 
o 
+1.386 


=e 





we i ae, Shag te, Sa 


o) 
3 
& 
@ 
ee 
aL 
Qe. 
= 
LP 









2.2704 142.31 | 0.384 ~0.367 19.957 





































0 1.8112 120.77 | 0.306 -0.362 22.966 
H, 1.8051 85.36 | 0.305 -0.362 152704 
H 0.0266 | 1.02 | 0.004 -0.022 0.132 
OH 0.0056 9 3% | 0,002 | 20,097 8.64 | 0,0 
559189 349 8G | 1.060 1.120 58.807 


S. = InyjSj - np.R. zed In mo =npeReln By 


S, = 349.80 - (5.9189) (1.987) (+1.120) = (5.9189) (1.987) (3.5264) 
8S, = 349.80 * 33.172 = 41.473 
S, = 321.50 61/(100 gm) (°K) 


— 100 t. = 16.8 
Mean Molecular Weight = See 95, 


Frey -TS 

F = (103,246 * 122,419) - (2650) (321.50) 
¥ = 419,173 - 851,975 

¥ = -832,802 calories/10C gm propellant 


#ithough this free energy is numerically greater than thet shown for the 


ayatem at 2700°K, the enthalpy balance was not completely satisfied st 
either temperatuse, Free onergy chunges became very smmll as equilibrium 
4s approached, If the exact temperatures were found, the free energy of 


the dissociated system containing H and OH would be slightly less. 
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ll. Calculation of Exhaust Conditions 
During expansion through the nozzle, the minor species H and OH will be 


converted to Hp and HO so that these, along with No, constitute the exhaust 
gas composition. This is known for this relatively simple example. However, 
for solving 8 general case on a high-speed digital computer, the same method 
of minimization of free energy previously employed would be utilized at 
constant entropy. At an exhaust temperature of 1250°K, the entropy of the 


system iss 










Specie (ad int | Bian Bi 









105.175 | 0.308 -0,363 


—L52h07 | 92308 
308.083 | 1.000 






Ss, = ins, - al la = * in P| 


S_ = 308.083 * (1.987) (5.9040) (1.093) 

S, = 308.083 + 12,822 

S, = 320.905 

S, = 320.905 a 321.50 =S, (Satisfactory Check) 


The enthalpy, heat capacity, and free energy are evaluated as: 







18.413 










#0. =105 007 19.149 
= | 8 Q 434590 
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F=H-TS 


= (4h.h53 = 105.007) ~ (1250) (320.905) 
F = 60,554 = 401,132 


F = -461,685 calories/105 gm propellant 


ropeliant — _ 16.95 


Mean Molecular Weight = 5, cnt moles of gas 


12. Calcuiation of Rocket Performance Parameters 


(a) Isp = 9.330 VHj-H, = 9.330 797.27 


(b) 


(2), 


(a) 


ve) 


(f) 


(g) 


(h) 


Isp = (9.330) (28.236) 
Isp = 263.0 seconds 


Ivap = f *Isp = ( ee ) {28340} ) = 10.90 1b-sec/in2 


u, = Ispeg = (263.0) (32.2) = 8465 ft/sec 


k cp 8 280 


Sy wR iy 7 584807 ~ (1.967) (59189) 18°58 





Pea 2 Eo = (529040) (1250) = 1, 
“t| Ten ae (1.15667) (263) a 


g/c* = olf a 


BT + 
450 
_ 2.250 | | _ 2 
li casoras| eA as [> 35| 


gafe* = (1.974) (0.00892) (0.591) = 6.00567 
= & = 
c* 0.00567 5680 ft/sec 


= Janeg = = 1.491 
ee ee 


te. & Po, Bele = (0.00567) (34.0) (24.25) = he69 


ae ce Pe 
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APPENDIX C 


Theoretical Performance Curves for the 


an Se Fes 


Various Propellant Systems 
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SPECIFIC IMPULSE, lb-sec/1b 


TING Isp 
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LIQUID HYDROGEN - LIQUSD OSYCEN PRRPUtu iG CURVES » TCURE 1 
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LIQUID HYDROGEN - LIQUID OXYGEN PERFORMANCE CURVES FIGURE 2 
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Chamber Temperature, OX 


c# £t/sec 


Chanber 


Molecular Weisht 


LIQUID HYDROGEN ~ LIQUID OXYGEN PERFORMANCE CURVES FIGURE 3 
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SHIFTING Isp 


SPRCIFIC IMPULSE, lb-sec/1b 
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LIQUID HYDROGEN ~ biQUS0 FIUURTAD CRG ts 0 CYR FP DGURE 4 
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LIST OF NOMENCLATURE 


4 = cross-sections] area of the nozzle at any goneral position 
4 = a chemical reactant in the general equation of equilibrium 
A! = frontal area of an individual particle in multiphase flow. 
4* = area of the nozzle at the throat, i.e. et Mech i 
| & = area of the nozzles at ite exit position 
a j = ateas of the i th element present in the j th exhaust specie 
A j=8 numerical coefficient which relates an incremental chenge in 
concentration to an incremental change in free energy 
& = the number of moles cf reactant 4 in the general equation of equilibrium 
@ = an ompirical coefficient relating enthalpy, entropy, and heat capacity 
to temperature 
a = acceleration of a solid particl* during Zlow through the nozzle 
®), = mass fraction of the ges phase in miltiphase flow. 
> = mass fraction of the condensed phase in muitiphase flow 
. a; j = a coefficient which relates ti2 amount cf an exbaust specie present to 
the maximum amount which coul. thesretically exist, according to the 
mass baience 


B = a chemical reactant in the general equation .f equilibrium 

B, = the maxima number of moles of a gives specie which can theoretically 

ah exist for a given propellant compesition 

ay Bo = the maximum total moles of gas which can theoretically exist for a 

given propellant canposition 

B j74 numerical coefficient which .. lates at incremental m<cry. 4. 

| concentration to one in free ene gy 

| b = the mmber of moles of reavtent Rin the ¢eneral equation of equilibriuy 

! b = an arbitrary parameter proportional to boost velocity 

b = an empirical coefficient relating enthalpy, entropy, and heat capaci.7 
to temperature 
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| Gea chenfeal product of combustion in ths generel equation of equilibr’ 
4 Q*, = @ numerical coefficient which relates an incremental change ia 
| :  gonceiitratioa to one in free energy 

Cp = drag coefficient 

c D.. the number of moles of product C in the general equation of equiiihviun 

@ © en empirical coefficient relating enthalpy, entropy, and heat COPEL - ty 

to temperature 
c * velocity of presaure wave propagetion 
1 c* = charecteristic exhaust velocity, or the velocity of sound et the ;. zzle 
a throat 
= | «Gy = thrust coefficient 
= CG, = specific heat. capacity at constant pressure 
| Gy * specific heat capacity at constant volume 


ea 


yoy ; Cy = discharge coefficient 
.8V¥g = constent colume specific heat capacity of the gas phase in nai :° hase 
flow 
“Dg = conesant pressure specific heat capacity of the gas phase in! - ti:nase 
flow 
{ 





ee * constant pressure specific heat capacity of the condensed pha:- in 

t multiphase flow 

} ep, - constant pressure specific heat capacity of a newly-formed cr... sused 
phat, otther aolid or Vicaid_ 
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= a chemical product of canbustion tn the general ecuation of equilibriun 

= substantial derivative 

= the number of moles of product Din the general equation of equilibrium 

= an empirical coefficient relating enthalpy, entropy, and heat capacity 
to temperature 

= standard derivative 


wn eauoys 


= internal energy of the system 

= a chemical product of canbustion in the general equation of equilibrium 

= energy emitted fran an individual condensed particle hy rudietion 

= total grameatoms of the i th clement present in a given propellant 

formiletion 

® = the number of moles of product E in the general equation of equilibrium 

e = an empirical coefficient relating enthalpy, entropy, and heat capacity. 
to t¢ nperature 

© = base of netural logarithme 


bj 1 ted 
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F = Gibbs foe energy of tie systex 

F =a chemical product of combustion in the general equation of equilibrium 
FO = ideal free energy at i atm presenre 

Ft = thrust 

Fy, = field forces 

¥, = force acting on an individual condensed particle in mltiphase flow 
al =the number of moles of product ¥ in the general equation of aquilibriun 
f = mass ratio of condenged species to gae species in multiphase flow y 


&@ = asceleration of gravity (32.17% 24/sec“) 
&’ = Local eeceleration of gravity 
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E = enthalpy of the system 

K' = reference enthalpy at a refesance temperature 

H, = stegnation enthalpy of tie ayztem 

H; = enthalpy of the combustion produci# in %°- chamber 

Re = enthalpy of the combustion produces gs tho nozzle exit 

Hy = enthalpy of the jth specie 

AH, = moler latent heat of phase chenge 

Aue, = heat-of-formation of the j th camponent at a stenderd reference 
J temperature 

& = avorege heat transfer coefficient 

Ryo = height of a missile at time of burnout 
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I * total impulse, lb-sec 

Isp = specific impulse, lb-sec/1b 

Ivsp = volumetric specific impulse 

i = designation for any elerent ip a propellant formulation 


J = mechenical equivalent of heat 
rues j = designation for any product of combustica 
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K = thermodynamic equilibrium constant 
= empirical constant 
2 = thermodynamic equilibrium pressure constant 
= thermodynamic equilibrium concentration constant 
= mean ratio of specific heat capacities, C p/Gys at any general time 
= mean specific heat ratio during nozzle espana on 


ln = natural logarithm to the base ¢ = 2.718 


M = Mach number 


Mj = moles of 4 given apecis j at any general time 
Mg = total moles of gas existing at any general time 
MW = 


mean molecular weight of the canbustion products 

MW* = mean molecular weight of the cambustion products at the throat, 
i.e. at Mach 3% 

m = masa rate of flow 

m = designation for any gaseous element in the products of cambustion 

m’ = mass of an individual condensed particle in multiphase flow 


Nas Nps Ng.+--M; = molar equilibrium concentrations of the respective 
species in the general equation of equilibrium 
Np + total moles of gas 
n, = .4e8 of gas in multiphase flow 
= designation for any guseous molecule in the products of combustion 
a, = moles of the j th specie 
maximum possible moles of gas which can theoretivally te formed 
from a given propellant composition 
An = total moles of products, minws tetal moles of resctants 


"Tnax = 


.,a°tat of nozzle cross-section at any general position 
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x = préasi.tea 

i 8 Gohre: oe es <a Yor say condeasad ora in the oroducta of canbustica 
Bg * est reL rer ere of is aurrouacina rinis medium 

Py * JFeaeae : * $a wembuttion sxecdibe.t 

Er. * prewsui< < ? the Senta yeseuets av the nozzle exit 

& < ebaganvd Se acne we the Fludad (po) = ty) 

oF x plese ss uhe vheocet of the nugale, 1.6, at Mach 1 
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q = energy absorbed by She 23% 10 Atom bis s sondusdinge 
q' = total heat transferred between the gas ae sowdensed phases in milti- 


phase flow — s é 
Q; = reversible energy change at const ant temperature 
R = universal gas constant e e 
Re = Reynolds’ number a oo 


Ps R/M 
= radius of an individual condensed partic “in mul ci plage flow 


109, 
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& = absolute entropy of the syz<em 

S* = ideal entropy at 1 atm pressure 

S' = reference entropy at a reference temperature an 1 atm pressure 

S, = entropy of the j th specie 

Ss. = horizontal diatence of a missile from ites launching position at time of 
burnaut 

® = empirical constant dependent upon missile design conditions 


T = temperature of the system 
tT, = stagnation temperature of the system 
Ts = combustion or chamber temperature (T, = t,) 
T. = temperature of the products of ecmbustion at the nozzle exit 
T* = temperature of the system at ths nozzle throat, i.e. et Mach 1 
a8 = température of the gas phase in multiphase flow 
2 = temperature of the condensed phase in 11]tiphase flow 
= time 
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U = unit tensor 


“V_= volume of the fluid system at any general tine 
vis syecific voluws of fluid system 
v*, 3 = apecific volurs of the gas phase in mal these flow 
V,°= volume of propeiiant ia a missile 
o = velocity vector of the fluid syaten 
v = fluid velocity at «ny general position 
v' = effective exhaust relocity 
Ve = velocity of the products Of cazbustion at the nezzle exit 
_v® = fluid Welocity st the ‘nozzle throat, ice. ak Mech 2 
“aex = maximum attainable 6riwust velocity i 
Vx * one-dimensional contribution to fluid v."ocity (vy = =v for 
_ whidirectional flow) 
Yg = velocity of the gas phese in multinhese flow. ae 
Vg = velocity of the condensed phase in multiphase flow eh 
‘go = velcecity of the gas phase at the instant a condensed particle enters 
the gas streen 
Vp = incremental doort velocity imp: ted to a missile 


We = mase of missile hardware, or the tote] mess of a missile at burnout 
Wy = total initial mass of the missile hardware and its propellant 
vn = total mass of the missile (W, = Wy) 
= total mass of propellant consumed duxing combustion in time t. 
w = non-expansion work done by the system 
w = work done by the system on its surroundings 


x = linear distance along the nozzle axis 

x = general concentration of a specie in the products of caabustion 

xj = the ratio of the moles of specie j in the cambustion products to the 
total possible moles of specie j which could theoretically be formed 

Xq = the ratio of the total moles of gas in the combustion products to the 
total possible moles of gas which could theoretically be formed 
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Y © general uathematical representation of the free energy of mixing eid 


pressure as a function of concentration 


y ‘empirical cosfficient dezendent upon physical conditions of the system 


end surrounding, ¢.g- missiie surface area, density of surrounding 


medium, ste, 
% * potential energy of positioa 
o © half-engle. of nozzle divergence 


(f° = density of the fluid syetem or of the propellant 


f, = stagnation density of the fluid system 


x 
f = density of the fluid system at the nozzle throst, i.e. at Mach I 


Pg = density of the ges phase in multiphase flow 
fa = density of the condensed phase in multiphase flow 


AL = fluid viscosity 


Ag = viscesity of the gas phase in multiphase flow 
$@ = meen ratio of the specific heat capacities, C)/Cy, for multiphase flow 


V = vector operator 
A = incremental cueage in a property 
@ = Stefen-Boltanegn constant 


= anissivity coustant 


™ engic of missile trajectory with the horizontal 


€ 

p 

6 = partial derivative 

Y = ideal cycle efficiency — 
T 
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